void daxpy(long,double,double[],long,double[],long);
void dblse(double*,long,long,long,long,double[],double[][2],
	double,long,double,long*,double[],double*,long*,
	double[],double[],double[],double[],double[],double[],long[],
	long[],double[],double[],double[]);
void dcopy(long,double[],long,double[],long);
double ddot(long,double[],long,double[],long);
void derv1(STRING,double,byte);
void dhtcc(long,long,long,long,double[],double*,double[],
	long,long);
double dnrm2(long,double[],long);
double dranu(void);
void drotg(double*,double*,double*,double*);
void dswap(long,double[],long,double[],long);
void ermor(STRING,byte);
void ermsg(STRING,long,long,STRING,byte);
void ierm1(STRING,long,long,STRING,STRING,long,byte);
void ierv1(STRING,long,byte);