/*opfun  routine used to get continuum pumping of lines */
double opfun(double x);