pdens function

Auxiliary functions for pitman