rhopot function

Calculates potential density