void drot(long,double[],long,double[],long,double,double);