void dcomqr(long,long,long,long,double*,double*,double[],
	long*);
double dzabs(double[]);
void zquo(double[],double[],double[]);
void zsqrtx(double[],double[]);