y[i] *= (1.0+SPREAD*ndev.dev());
sig[i]=SPREAD*y[i];
}
FGauss fgauss = new FGauss();
Fitmrq myfit=new Fitmrq(x,y,sig,guess,fgauss);
myfit.fit();
// System.out.printf(setw(18) << "chi-squared:%f\n", setw(13) << myfit.chisq);
// System.out.printf(fixed << setprecision(6);
// for (i=0;i<MA;i++) System.out.printf(setw(9) << myfit.a[i];
// System.out.printf(endl;