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