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