int nmmm(double Hc,double Mr,double Hg,double d,double s,double D,
    double *a1,double *a2,double *a3);