double dcpval(double[],long,double);
void dpfit(long,double[],double[],double[],long,LOGICAL32,LOGICAL32,LOGICAL32,
      double[],long*,double*,double*);