void drotmg(double*,double*,double*,double,double[]);