double dcdnml(double,double,double);
double dppnml(double,double,double);