for (j=0;j<NP;j++) y[i] += amp[j]*f[j];
y[i] += ndev.dev();
sig[i]=stdev;
}
Fitsvd myfit = new Fitsvd(x,y,sig,this,TOL);
myfit.fit();
for (j=0;j<NP;j++) {
// System.out.printf(myfit.a[j] << " %f\n", amp[j] << " " ;
// System.out.printf(sqrt(myfit.covar[j][j]));