LISTING OF PROGRAM IPOT

1 PROGRAM

2 SIGNIFIC 8

3 PRINT "RITZ PARAMETER DETERMINATION OF THE IP OF SODIUM"

4 PRINT POW RZ

5 PRINT " "

6 PRINT "ITERATION TO SEARCH-FIT IP"

7 PRINT " "

8 IP1=MAX(E)+RZ/(MAX(N)**2)

9 FOR I=1,26

10 IP=IP1

11 X=(IP-E)/RZ

12 Y=N-1/SQRT(X)

13 DY=DE/(2*(IP-E)*SQRT(X))

14 W=1/DY**2

15 A=LSQPOL(X,Y,POW,R,W,CHI2,E1)

16 IPARR(I)=IP

17 CHIARR(I)=CHI2

18 X1=X

19 FOR J=1,40

20 YF=POLYVAL(A,X1)

21 X1=1/(N-YF)**2

22 ENDLOOP J

23 PRINT IP CHI2

24 IP1=AVERAGE(E+RZ/(N-YF)**2)

25 ENDLOOP I

26 END

EXECUTE IPOT

RITZ PARAMETER DETERMINATION OF THE IP OF SODIUM

POW = 2 RZ = 109734.72

ITERATION TO SEARCH-FIT IP

IP = 41328.077 CHI2 = 5534110.1

IP = 41368.194 CHI2 = 2551411.3

IP = 41394.849 CHI2 = 1174073.7

IP = 41412.681 CHI2 = 539686.89

IP = 41424.663 CHI2 = 247873.58

IP = 41432.738 CHI2 = 113749.14

IP = 41438.191 CHI2 = 52143.359

IP = 41441.878 CHI2 = 23867.459

IP = 41444.373 CHI2 = 10901.778

IP = 41446.063 CHI2 = 4964.5128

IP = 41447.208 CHI2 = 2251.0773

IP = 41447.984 CHI2 = 1014.6154

IP = 41448.51 CHI2 = 453.64919

IP = 41448.866 CHI2 = 200.83152

IP = 41449.108 CHI2 = 88.048373

IP = 41449.272 CHI2 = 38.535645

IP = 41449.383 CHI2 = 17.358384

IP = 41449.458 CHI2 = 8.6976701

IP = 41449.509 CHI2 = 5.4451649

IP = 41449.544 CHI2 = 4.4442938

IP = 41449.567 CHI2 = 4.319414

IP = 41449.583 CHI2 = 4.4893286

IP = 41449.594 CHI2 = 4.7215856

IP = 41449.601 CHI2 = 4.9328789

IP = 41449.606 CHI2 = 5.1008855

IP = 41449.609 CHI2 = 5.2261761

LISTING OF PROGRAM FIT

1 PROGRAM

2 NFIT=INTS(3,20)

3 XFIT=1/NFIT**2

4 FOR I=1,20

5 YFIT=POLYVAL(A,XFIT)

6 XFIT=1/(NFIT-YFIT)**2

7 ENDLOOP I

8 PRINT " "

9 EFIT=IP-RZ/(NFIT-YFIT)**2

10 PRINT "FITTED VALUES",IP,A

11 TABULATE NFIT E EFIT DE Y YFIT DY

12 END

EXECUTION FIT

FITTED VALUES IP = 41449.567

A (A 3 Component Array)

.85477043 .10962661 .068169309

NFIT E EFIT DE Y YFIT DY

**** ********* ********* *** ********* ********* ************

3 16973.368 16973.368 .03 .88261399 .88261399 1.2976194E-6

4 30272.58 30272.587 .05 .86664467 .86664361 7.0084971E-6

5 35042.85 35042.816 .05 .86139247 .86140324 1.6149485E-5

6 37297.611 37297.623 .05 .85902338 .85901587 3.0955145E-5

7 38540.93 38540.958 .05 .85775348 .85772407 5.2793163E-5

8 39298.84 39298.886 .05 .85702178 .85694518 8.3029788E-5

9 39795 39794.882 .1 .85614864 .85643898 2.4610218E-4

10 40137.231 40137.125 .1 .8557213 .85609133 3.4839692E-4

11 40383.16 40383.187 .1 .85596905 .8558422 4.756173E-4

12 40566.031 40566.009 .1 .85551953 .85565754 6.3067473E-4

13 40705.545 .85551685

14 40814.456 .8554072

15 40901.091 .85532007

16 40971.136 .85524968

17 41028.571 .85519201

18 41076.251 .85514417

19 41116.267 .85510403

20 41150.178 .85507003