attach rs1.new
xrs1
call public $ez_readfile
fit function "a-b*exp(-(x-xo)/t)" to curve 1 of graph gpo218
 set printer to 'lpt80'
printout table foo noheader on file 'foo'
 set printer to 'lpt132'
printout table longfoo on file 'lfoo'
 % enscript -r -Pprinter lfoo
 # set plotter to 't4010'
 # plot graph foo bottomkey on file 'foo.tek'

          % ps4014 foo.tek | lpr -Pprinter
 /mit/rs1.new/rsfix


fit function "sqrt(c1*(x-x0)**2 +c2)" to curve 1 of gop2f1
plot graph gop2f1 bottomkey on file 'gop2f1.tek


4/4/90

columns:
0  Iz(mA)
1  sigma_Iz(mA)
2  f_low (kHz)
3 sig_f_low (kHz)
4 gmu_85(MeV/G)
5 sigma_gmu_85 (MeV)
6 1/(5)^2
7 (4)/(5)^2
8 f_hi (kHz)
9 sig_f_hi (kHz)
10 gmu_87
11 sig_gmu_87
12 1/(11)^2
13 (10)/(11)^2

!set col 4 of gmu to 8.944e-16*col 2/abs(col 0)

!set col 5 of gmu to c*sqrt((1.994e3*col 3/col 0)**2
			+ (col 2*5555*0.05/col 0)**2
			+ (1.944e3*col 2*col 1/(col 0)**2)**2)

bad:
!do n = 1 to 15;
sumtop = row n of col 4/ (row n of col 5)**2 + sumtop;
end;

!set col 6 of gmu to 1/((col 5)**2)
!set col 7 of gmu to (col 4/((col 5)**2))


!set col 10 of gmu to 8.944e-16*col 8/abs(col 0)

!set col 11 of gmu to c*sqrt((1.994e3*col 9/col 0)**2
			+ (col 8*5555*0.05/col 0)**2
			+ (1.944e3*col 8*col 1/(col 0)**2)**2)
!set col 12 of gmu to 1/((col 11)**2)
!set col 13 of gmu to col 10 * col 12

fit function "c1*abs(x-x0)" to curve 1 of ggmu
plot graph gop2f1 bottomkey on file 'gop2f1.tex

Result of (WEIGHTED) FIT is
    2.330959*ABS(X-24.292367)

FIT converged.  This does NOT necessarily imply a good fit.
Final SUM OF SQUARES is 33.082933.
0    1 PARAM.  2 FINAL    3 STD     4 T VALUE  5 SIG.  6 DEP.
                 VALUES     ERRORS               LEV.    VALUES
---------------------------------------------------------------
1    C1         2.330959  0.007068      329.8  0.0001   0.09623
2    X0        24.292367  0.185875      130.7  0.0001   0.09623

Result of (WEIGHTED) FIT is
    3.502268*ABS(X-24.169579)

FIT converged.  This does NOT necessarily imply a good fit.
Final SUM OF SQUARES is 85.496399.

                        PARAMETER TABLE

0    1 PARAM.  2 FINAL    3 STD     4 T VALUE  5 SIG.  6 DEP.
                 VALUES     ERRORS               LEV.    VALUES
---------------------------------------------------------------
1    C1         3.502268  0.011496      304.6  0.0001    0.1076
2    X0        24.169579  0.200141      120.8  0.0001    0.1076


fit function "5.3 + c2*exp(-(x-x0)/t0)" to curve 1 of gbrev

fit function "1.47 + c2*exp(-(x-x0)/t0)" to curve 2 of gbrev

fit function ".81 + c2*exp(-(x-x0)/t0)" to curve 3 of gbrev
