void dqrbd(long*,double[],double[],long,double*,long,long,
	double*,long,long);
void drot(long,double[],long,double[],long,double,double);
void drotg(double*,double*,double*,double*);
void dswap(long,double[],long,double[],long);
void ermsg(STRING,long,long,STRING,byte);