void dcdiv(double,double,double,double,double*,double*);
void devbh(double*,long,long,long*,long*,long[],double[]);
void devvun(double*,long,long,double[],double[],double*,long[],
	double[]);
double dnrm2(long,double[],long);
void dscal(long,double,double[],long);
void ermsg(STRING,long,long,STRING,byte);