LISTING OF PROGRAM KEPLERTA
1 PROGRAM
2 L2=(3/2)**2
3 FOR I=1,18
4 E=-0.27 + 0.01*I
5 R=INTS(450,12000)/500
6 P2=E+2/R-L2/R**2
7 RPEE=R
8 CONSTRAIN(P2,R:P2.GE.0)
9 DR=R(2)-R(1)
10 P=SQRT(P2)
11 PI=ACOS(-1)
12 LEFT=(DR/PI)*SUM(P)
13 PER=MIN(R)
14 AP=MAX(R)
15 EEE(I)=E
16 PEREE(I)=PER
17 APEE(I)=AP
18 LEFTSIDE(I)=LEFT
19 ENDLOOP I
20 TABULATE EEE PEREE APEE
LEFTSIDE
21 END
EEE PEREE APEE
LEFTSIDE
**** ***** ******
********
-.26 1.37 6.322
.461159
-.25 1.356 6.644
.499996 LHS(0+1/2) gives E=-1/4
-.24 1.342 6.992
.541242
-.23 1.328 7.366
.585145
-.22 1.316 7.774
.632009
-.21 1.304 8.220
.682182
-.2 1.292 8.708
.736067
-.19 1.282 9.244
.794158
-.18 1.272 9.840
.857019
-.17 1.260 10.504
.925355
-.16 1.250 11.248
.999993
-.15 1.242 12.092
1.08199
-.14 1.232 13.054
1.17261
-.13 1.224 14.162
1.27350
-.12 1.214 15.452
1.38675
-.11 1.206 16.976
1.51511 LHS(1+1/2) gives E=-1/9
-.1 1.198 18.802
1.66228
-.09 1.190 21.032
1.83333
n = n_r + ell + 1