double dcdnml(double,double,double);
double derfc(double);