void dsortp(double[],long,long,long[]);
void dsortq(double[],long,long,long[]);