void dsfind(double[],long,long,double,long*,long*);
double dsquad(long,long,double[],double[],double,double);
double dsval(long,long,double[],double[],double,long);