/*qintr integrates Q for any continuum between two limits, used for normalization */
double qintr(float *qenlo, 
  float *qenhi);