LISTING OF PROGRAM CALLI

1 PROGRAM

2 FOR IE=1,31

3 E=3.75+0.25*IE

4 EXECUTE TRYI

5 EI(IE)=E

6 AM(IE)=R1

7 AP(IE)=R2

8 LEFT(IE)=EBK

9 ENDLOOP IE

10 TABULATE EI, AM, AP, LEFT

11 END

LISTING OF PROGRAM TRYI

1 PROGRAM

2 R0=8*INTS(1,NUM)/NUM

3 P20=E-R0**2/2-(49/8)/R0**2

4 NULLS=ROOTS(P20,R0)

5 FREE R0, P20

6 R1=NULLS(1)

7 R2=NULLS(2)

8 DELTA=(R2-R1)/NUM

9 FOR I=1,NUM-1

10 R(I)=R1+I*DELTA

11 ENDLOOP I

12 P2=E-R**2/2-(49/8)/R**2

13 P=SQRT(P2)

14 EBK=DELTA*SUM(P)*SQRT(2)/PI

15 END

EXECUTE CALLI

EI AM AP LEFT EI AM AP LEFT

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

4 1.4365 2.4365 .25 8 .89792 3.8979 2.25

4.25 1.3561 2.5809 .375 8.25 .88274 3.9649 2.375

4.5 1.2929 2.7071 .5 8.5 .86835 4.0306 2.5

4.75 1.2404 2.8216 .625 8.75 .85469 4.0951 2.625

5 1.1955 2.9276 .75 9 .84169 4.1583 2.75

5.25 1.1562 3.0271 .875 9.25 .82929 4.2205 2.875

5.5 1.1213 3.1213 1 9.5 .81746 4.2816 3

5.75 1.0899 3.2112 1.125 9.75 .80614 4.3417 3.125

6 1.0614 3.2975 1.25 10 .7953 4.4009 3.25

6.25 1.0353 3.3805 1.375 10.25 .7849 4.4591 3.375

6.5 1.0113 3.4608 1.5 10.5 .77492 4.5166 3.5

6.75 .98909 3.5386 1.625 10.75 .76533 4.5732 3.625

7 .96841 3.6142 1.75 11 .75609 4.6291 3.75

7.25 .9491 3.6877 1.875 11.25 .74719 4.6842 3.875

7.5 .93099 3.7594 2 11.5 .73861 4.7386 4

7.75 .91397 3.8294 2.125

LISTING OF PROGRAM TRY3

1 PROGRAM

2 FREE P2,R,P20,R0

3 PRINT NUM

4 R0=8*INTS(1,NUM)/NUM

5 P20=9/2-R0**2/2-(49/8)/R0**2

6 NULLS=ROOTS(P20,R0)

7 R1=NULLS(1)

8 R2=NULLS(2)

9 PRINT R1 R2

10 DELTA=(R2-R1)/NUM

11 FOR I=1,NUM-1

12 R(I)=R1+I*DELTA

13 ENDLOOP I

14 P2=9/2-R**2/2-(49/8)/R**2

15 F=DELTA/SQRT(P2)

16 INTSUM=SUM(F)

17 CTHE=27*PI

18 PRINT INTSUM, CTHE

19 FF=F/INTSUM

20 AVR(1)=SUM(FF*R**2)

21 AVR(2)=SUM(FF*R**4)

22 AVR(3)=SUM(FF/R**2)

23 AVR(4)=SUM(FF/R**4)

24 EXECUTE EXPECT

25 TABULATE AVR TAVR

26 END

LISTING OF PROGRAM EXPECT

1 PROGRAM

2 NR=0

3 L=3

4 GUM=SQRT(2*(NR+L+1))

5 BUM=SQRT(2*NR+1)

6 AA=(GUM+BUM)/SQRT(2)

7 BB=(GUM-BUM)/SQRT(2)

8 RAT=AA/BB

9 PRINT AA BB RAT

10 DUM=AA**2+BB**2

11 DUM2=AA*BB

12 TAVR(1)=DUM/2

13 TAVR(2)=(3*DUM**2-4*DUM2**2)/8

14 TAVR(3)=1/DUM2

15 TAVR(4)=DUM/(2*DUM2**3)

16 END

EXECUTE TRY3

NUM = 10000

R1 = 1.2929 R2 = 2.7071

INTSUM = 2.2011 CTHE = 84.823

AA = 2.7071 BB = 1.2929 RAT = 2.0938

AVR TAVR

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

4.4952 4.5

24.17 24.25

.28535 .28571

.10447 .10496