void dgefsc(double*,long,long,double*,long,long,long[],double*,double[]);
void dmatp(double*,long,long,long,STRING);