float sasum(long,float[],long);
void saxpy(long,float,float[],long,float[],long);
float sdot(long,float[],long,float[],long);
void sgeco(float*,long,long,long[],float*,float[]);
void sgefa(float*,long,long,long[],long*);
void sscal(long,float,float[],long);