long isamax(long,float[],long);
void saxpy(long,float,float[],long,float[],long);
void sgefa(float*,long,long,long[],long*);
void sscal(long,float,float[],long);