void dsdif(long,long,double[],double[],long,double*);