void drotg(double*,double*,double*,double*);