void dcopy(long,double[],long,double[],long);
void dcov3(double*,long,long,double[],double,double[],long*);
double ddot(long,double[],long,double[],long);
void dscal(long,double,double[],long);
void ierm1(STRING,long,long,STRING,STRING,long,byte);