int bertero(double g,double p,double q,double x,double y,
            double *Hx,double *Hy,int CMD);