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