int main(int n, char **args)
{
double alphaSMZ,mbp,mbmb,mtp;
double Q;
int nf;
double McMc=1.4;
sscanf(args[1],"%lf",&alphaSMZ);
sscanf(args[2],"%lf",&McMc);
sscanf(args[3],"%lf",&mbmb);
sscanf(args[4],"%lf",&mtp);
sscanf(args[5],"%lf",&Q);
initQCD(alphaSMZ,McMc,mbmb,mtp);
printf("qMass : %E %E %E %E\n", qMass[3], qMass[4], qMass[5], qMass[6]);
printf("lambda: %E %E %E %E\n", lambda[3],lambda[4],lambda[5],lambda[6]);
printf("qMin=%E\n",qMin);
printf("MbPole=%f\n", MbPole);
printf("MtMt=%f\n", qMass[6]);
printf("qmass[6]=%f\n",qMass[6]);
printf("alphaS(%f)=%f\n",Q,alphaQCD(Q));
printf("MbRun=%f MbEff=%f \n",MbRun(Q),MbEff(Q) );
printf("MtRun=%f MtEff=%f \n",MtRun(Q),MtEff(Q) );
}
double MbEff(double Q) { double m=MbRun(Q)*sqrt(1+DeltaQCD(Q)); if(m>MbPole) return MbPole; else return m;}
double mbrun_(double *q){ return MbRun(*q);}