LISTING OF PROGRAM SCREEN
1 PROGRAM
2 RY=109737.3155-60.2/AM
3 AL=1/137.03602
4 C0=RY*(AL**2)/128
5 ALPI=AL/ACOS(-1)
6 GM2=0.0023192
7 PRINT " "
8 DATE
9 TIME
10 PRINT " "
11 PRINT " DATA
BASE"
12 TABULATE Z AM P32 P12 E2
DE2
13 PRINT " HIGHEST
DIRAC " TERM,LAMB,GREEN
14 PRINT " "
15 I=INTS(1,7)
16 WHERE (I.GT.TERM) I=0
17
DIRAC=(I/INTS(1,7))*COEFFDIR
18 CORR=1
19 FOR I=1,25
20 ZS=(E2/(C0*CORR))**.25
21 X=(AL*ZS)**2
22 EUCLID=1/(1-DIRAC(7)*X/COEFFDIR(6))
23 CORR=1+DIRAC(1)*X+DIRAC(2)*X*X+DIRAC(3)*X**3+DIRAC(4)*X**4
24 CORR=CORR+DIRAC(5)*X**5+DIRAC(6)*X**6+DIRAC(7)*EUCLID*X**7
25
CORR=CORR+LAMB*(GM2+(5/4)*ALPI*(X**2)*LOG(X**2)*GREEN)
26 NEXT I
27 WZ=DE2*ZS/(4*E2)
28 WZ=1/WZ**2
29 SZ=Z-ZS
30 XZ=1/ZS
31
AZ=LSQPOL(XZ,SZ,1,RR,WZ,CHI2)
32 ZZMIN=AZ(1)+2*SQRT(AZ(2))
33 PRINT " FITTING
PARAMETERS"
34 PRINT " "
35 PRINT "
AZ(1)="AZ(1),"
AZ(2)="AZ(2),ZZMIN
36 PRINT " "
37 ZFIT=INTS(31,56)
38 XZFIT=1/ZFIT
39 FREE CORR
40 FOR I=1,25
41 SZFIT=POLYVAL(AZ,XZFIT)
42 XZFIT=1/(ZFIT-SZFIT)
43 NEXT I
44 DEGFREE=NOELS(E2)-2
45
PROBAB=CHIPROB(CHI2,DEGFREE)
46 PRINT CHI2, PROBAB
47 ZSCR=ZFIT-SZFIT
48 X=(AL*ZSCR)**2
49 EUCLID=1/(1-DIRAC(7)*X/COEFFDIR(6))
50 CORR=1+DIRAC(1)*X+DIRAC(2)*X*X+DIRAC(3)*X**3+DIRAC(4)*X**4
51 CORR=CORR+DIRAC(5)*X**5+DIRAC(6)*X**6+DIRAC(7)*EUCLID*X**7
52
CORR=CORR+LAMB*(GM2+(5/4)*ALPI*(X**2)*LOG(X**2)*GREEN)
53
RYFIT=109737.412-60.2/(2.25*ZFIT)
54
E2FIT=(RYFIT*(AL**2)/128)*CORR*ZSCR**4
55 TABULATE ZFIT SZ SZFIT XZ
XZFIT
56 TABULATE ZFIT E2 E2FIT
57 END
EXECUTE SCREEN
October 8, 2007 9:25 PM
DATA BASE
Z AM P32
P12 E2 DE2
** ** ******
****** ******* *******
31 70 66887
65170 1717.8 8.589
32 73 84103
81315 2788 13.94
33 75 101245
97135 4110 20.55
34 79 118469
112771 5697.9 28.49
35 80 135854
128274 7580 37.9
36 84 153475
143697 9778 48.89
37 85 171409
159075 12334 61.67
38 88 189703
174437 15266 76.33
39 89 208433
189811 18622 93.11
40 91 227627
205202 22425 112.13
41 93 247344
220624 26720 133.6
42 96 267632
236085 31547 157.74
HIGHEST DIRAC TERM = 7
LAMB = 1 GREEN =
1
FITTING PARAMETERS
AZ(1)= AZ(1) = 9.8749 AZ(2)= AZ(2) = 100.78 ZZMIN = 29.952
CHI2 = 4.7013
PROBAB = .91022
ZFIT SZ SZFIT
XZ XZFIT
**** ****** ******
******* *******
31 17.103 17.153
.071957 .072215
32 16.321 16.29
.06378 .063652
33 15.731 15.7
.057908 .057804
34 15.27 15.25
.05339 .053332
35 14.894 14.885
.049736 .049714
36 14.582 14.58
.04669 .046684
37 14.313 14.318
.044077 .044088
38 14.082 14.09
.041809 .041823
39 13.876 13.888
.039803 .039822
40 13.695 13.708
.038015 .038034
41 13.531 13.546
.036405 .036424
42 13.383 13.398
.034944 .034963
43 13.264 .033629
44 13.141 .032405
45 13.027 .031276
46 12.921 .030231
47 12.824 .02926
48 12.732 .028355
49 12.647 .027508
50 12.567 .026714
51 12.492 .025969
52 12.421 .025266
53 12.354 .024603
54 12.291 .023976
55 12.231 .023382
56 12.174 .022818
ZFIT E2 E2FIT ZFIT E2 E2FIT
**** ******* ********
**** ******* ********
31 1717.8 1693.3
44 42852
32 2788 2810.7 45 49501
33 4110 4139.8 46 56849
34 5697.9 5723
47 64946
35 7580 7593.5 48 73839
36 9778 9782.6 49 83581
37 12334 12322 50 94226
38 15266 15245 51 105827
39 18622 18586 52 118444
40 22425 22380 53 132136
41 26720 26664 54 146966
42 31547 31477 55 162998
43 36859 56 180299