void dnqfj(long,double[],double[],double*,long*);
void dckder(long*,long,long,double[],double[],double*,long,double*,long*,long*,
      double*);
void dnqsol(void(*)(long,double[],double[],double*,long*),long,double[],
      double[],double,long[],double[],long);
double dnrm2(long,double[],long);