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