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