double ddot(long,double[],long,double[],long);
void dherql(double*,double*,long,long,double[],double*,double*,double[],long*);
void dmatp(double*,long,long,long,STRING);
void dvecp(double[],long,STRING);