/*pintr integrates L for any continuum between two limits, used for normalization */
double pintr(double penlo, 
  double penhi);