void snqfj(long,float[],float[],float*,long*);
void sckder(long*,long,long,float[],float[],float*,long,float*,long*,long*,float*);
void snqsol(void(*)(long,float[],float[],float*,long*),long,float[],float[],
      float,long[],float[],long);
float snrm2(long,float[],long);