double dasum(long,double[],long);
void daxpy(long,double,double[],long,double[],long);
double ddot(long,double[],long,double[],long);
void dgeco(double*,long,long,long[],double*,double[]);
void dgefa(double*,long,long,long[],long*);
void dscal(long,double,double[],long);