double ddot(long,double[],long,double[],long);
void dgeco(double*,long,long,long[],double*,double[]);
void dgei(double*,long,long,long[],double[]);
void durev(double*,double*,long,double*,long[],double*);
void ermsg(STRING,long,long,STRING,byte);