window -t wks Templates\index temp;
open -w C:\Documents and Settings\shawn\Desktop\academic\beams\033006\roomtemp\lowerpeak\index.txt;
files=wks.maxRows;


loop(i,1,files){

window -n wks temp;
open -w C:\Documents and Settings\shawn\Desktop\academic\beams\033006\roomtemp\lowerpeak\$(index_filenumber[i]);
col=wks.nCols;
row=wks.maxRows;
window -t wks Templates\mine4 a$(index_filenumber[i])fit;
loop(ii,1,row){
a$(i)fit_frequency[ii]=a$(i)_C[ii];
a$(i)fit_real[ii]=a$(i)_A[ii];
a$(i)fit_imag[ii]=a$(i)_B[ii];
a$(i)fit_error[ii]=2E-8;
};
window -c a$(i);
nlsf.init();
nlsf.keepparam=0;
nlsf.msgPrompt=0;
nlsf.constr$="";

sum(a$(i)fit_frequency);
resfreq=sum.mean;
resmin=sum.min;
resmax=sum.max;
sum(a$(i)fit_real);
xoff=sum.mean;
sum(a$(i)fit_imag);
yoff=sum.mean;
amp=.1;
gamma=30000;
xslp=0;
yslp=0;
theta=0;

nlsf.func$="resonancefirst";
nlsf.numDepend = 2;
nlsf.numIndep = 1;
nlsf.cleanupfitdata();
nlsf.fitData1$=a$(i)fit_real;
nlsf.fitData2$=a$(i)fit_imag;
nlsf.x$=a$(i)fit_frequency;





nlsf.derivStep = 0.0001;
nlsf.chiSqrErr=1;
nlsf.wType=3;
nlsf.w$=a$(i)fit_error;
nlsf.tolerance=1E-9;
nlsf.v5=0;
nlsf.v6=0;
nlsf.v3=0;
nlsf.v4=0;
nlsf.v1=1;
nlsf.v2=1;
nlsf.v7=1;
nlsf.v8=1;
nlsf.iterate(5);

nlsf.v5=0;
nlsf.v6=0;
nlsf.v3=0;
nlsf.v4=1;
nlsf.v1=1;
nlsf.v2=1;
nlsf.v7=1;
nlsf.v8=1;
nlsf.iterate(50);


nlsf.v5=1;
nlsf.v6=1;
nlsf.v3=1;
nlsf.v4=1;
nlsf.v1=1;
nlsf.v2=1;
nlsf.v7=1;
nlsf.v8=1;
nlsf.constraints=1;
nlsf.constr$="gamma>20000;resfreq>resmin;resfreq<resmax";
nlsf.fit(50);

window -c Graph1;
worksheet -n 5 realfit;
worksheet -n 6 imagfit;

worksheet -s 2 0;
worksheet -p 201 mygraph a$(i)graph;
window -a a$(i)graph;
layer -i a$(i)fit_realfit;
layer -a;
legend;
layer -s 1;
layer -i a$(i)fit_imag;
layer -i a$(i)fit_imagfit;
layer -a;
legend;

nlsf.end(1);


window -t wks templates\results results;
results_filenumber[i]=index_filenumber[i];
results_dcvoltage[i]=index_dcvolts[i];
results_acvoltage[i]=index_acvolts[i];
results_amplitude[i]=amp;
results_amperr[i]=nlsf.e4;
results_resonance[i]=resfreq;
results_resonanceerr[i]=nlsf.e5;
results_gamma[i]=gamma;
results_gammaerr[i]=nlsf.e6;
results_theta[i]=theta;
results_thetaerr[i]=nlsf.e3;
results_chisquared[i]=nlsf.chisqr;
};