void daxpy(long,double,double[],long,double[],long);
void dgefa(double*,long,long,long[],long*);
void dscal(long,double,double[],long);
long idamax(long,double[],long);