void dmpdrv(double[],long,double[],long*);
void dmpint(double[],long,double[],long*);
double dmpval(double[],long,double);