double dpquad(long,long,double[],double*,double,double);
void dsfind(double[],long,long,double,long*,long*);