double dasum(long,double[],long);
void daxpy(long,double,double[],long,double[],long);
void dcopy(long,double[],long,double[],long);
void ddasf(double,double[],double[],double[],double[],long,double,
	long,double[],long[]);
void ddasls(void(*)(),long,double*,double[],double[],long[],double,
	double,double*,long*,long,long*,double[],long,
	long[],long);
void dgbfa(double*,long,long,long,long,long[],long*);
void dgbsl(double*,long,long,long,long,long[],double[],
	long);
void dgefa(double*,long,long,long[],long*);
void dgemv(byte,long,long,double,double*,long,double[],
	long,double,double[],long);
void dgesl(double*,long,long,long[],double[],long);
void dscal(long,double,double[],long);
void dswap(long,double[],long,double[],long);