double dpquad(long,long,double[],double*,double,double);
void dprpl(double,byte,byte[],long,double,double,LOGICAL32);
double dpval(long,long,double[],double*,double,long);
void dsfit(double[],double[],double[],long,long,long,double[],double[],double*,
      long*,long,double*);
double dsquad(long,long,double[],double[],double,double);
void dstop(long,long,double[],double[],double*,long*,double[],double*);
double dsval(long,long,double[],double[],double,long);