int szczech(double p,double q,double g,double x,double y,
        double usp,double usq,double *Hx,double *Hy,int CMD);