TEST OF THE CONTINUED-FRACTION COULOMB & BESSEL PROGRAM - COUL90 WHEN LAMBDA IS AN INTEGER (L-VALUE) F IS REGULAR AT THE ORIGIN ( X = 0 ) WHILE G IS IRREGULAR ( => -INFINITY AT X = 0 ) L F(ETA,X,L) G(ETA,X,L) D/DX (F) D/DX (G) ETA = -.500 X = 20.000 XLMIN = .00 KFN = 0 (Coulomb) 0 -1.0237230D-01 -9.8257145D-01 -1.0068693D+00 1.0431459D-01 50 2.0210828D-15 1.0720740D+14 4.7241181D-15 -2.4419563D+14 ETA = -.500 X = 200.000 XLMIN = .00 KFN = 0 (Coulomb) 0 8.1994516D-01 -5.7026249D-01 -5.7168125D-01 -8.2199599D-01 50 -9.2809974D-02 -1.0110000D+00 -9.8083865D-01 9.0206999D-02 ETA = -.500 X = 200.000 XLMIN = 50.00 KFN = 0 (Coulomb) 50 -9.2809974D-02 -1.0110000D+00 -9.8083865D-01 9.0206999D-02 ETA = .000 X = 20.000 XLMIN = .00 KFN = 0 (Coulomb) The 3 sets of results are COUL90(KFN), SBESJY & (1/X) RICBES 0 9.1294525D-01 4.0808206D-01 4.0808206D-01 -9.1294525D-01 0 4.5647263D-02 -2.0404103D-02 1.8121740D-02 4.6667468D-02 0 4.5647263D-02 -2.0404103D-02 2.0404103D-02 4.5647263D-02 50 1.1300162D-15 1.9085083D+14 2.6533576D-15 -4.3681189D+14 50 5.6500808D-17 -9.5425417D+12 1.2984284D-16 2.2317722D+13 50 5.6500808D-17 -9.5425417D+12 1.3266788D-16 2.1840594D+13 ETA = .000 X = 200.000 XLMIN = .00 KFN = 0 (Coulomb) The 3 sets of results are COUL90(KFN), SBESJY & (1/X) RICBES 0 -8.7329730D-01 4.8718768D-01 4.8718768D-01 8.7329730D-01 0 -4.3664865D-03 -2.4359384D-03 2.4577708D-03 -4.3543068D-03 0 -4.3664865D-03 -2.4359384D-03 2.4359384D-03 -4.3664865D-03 50 8.1836187D-01 -6.0313111D-01 -5.8372987D-01 -7.9174553D-01 50 4.0918094D-03 3.0156556D-03 -2.9391084D-03 3.9436494D-03 50 4.0918094D-03 3.0156556D-03 -2.9186494D-03 3.9587277D-03 ETA = .500 X = 20.000 XLMIN = .00 KFN = 0 (Coulomb) 0 -8.1319612D-01 6.0386591D-01 5.8912930D-01 7.9223803D-01 50 6.2842419D-16 3.4159736D+14 1.4822479D-15 -7.8556494D+14 ETA = .500 X = 200.000 XLMIN = .00 KFN = 0 (Coulomb) 0 9.1775470D-01 -4.0029450D-01 -3.9929828D-01 -9.1545496D-01 50 9.4121859D-01 3.8776478D-01 3.7403121D-01 -9.0835846D-01 ETA = -5.200 X = 1.000 XLMIN = .00 KFN = 0 (Coulomb) 0 -2.4254092D-01 4.8973701D-01 1.5854720D+00 9.2164110D-01 50 8.8135718D-78 1.1258996D+75 4.4850689D-76 -5.6166364D+76 ETA = -5.200 X = 30.000 XLMIN = .00 KFN = 0 (Coulomb) 0 -3.9680743D-01 -8.3923161D-01 -9.7470166D-01 4.5865966D-01 50 7.8285385D-06 5.2407113D+04 9.7746374D-06 -6.2302749D+04 ETA = -5.200 X = 1000.000 XLMIN = .00 KFN = 0 (Coulomb) 0 -9.6858590D-01 2.3807890D-01 2.3931122D-01 9.7361013D-01 50 -5.7254138D-01 8.1749270D-01 8.2069434D-01 5.7478529D-01 j,y Bessels X = 1.000 XLMIN = .00 KFN = 1 (SphBess) The 3 sets of results are COUL90(KFN), SBESJY & (1/X) RICBES 0 8.4147098D-01 -5.4030231D-01 -3.0116868D-01 1.3817733D+00 0 8.4147098D-01 -5.4030231D-01 -3.0116868D-01 1.3817733D+00 0 8.4147098D-01 -5.4030231D-01 5.4030231D-01 8.4147098D-01 100 7.4447277-190 -6.6830795+186 7.4443610-188 6.7495744+188 100 7.4447277-190 -6.6830795+186 7.4443610-188 6.7495744+188 100 7.4447277-190 -6.6830795+186 7.5188083-188 6.6827436+188 j,y Bessels X = 100.000 XLMIN = .00 KFN = 1 (SphBess) The 3 sets of results are COUL90(KFN), SBESJY & (1/X) RICBES 0 -5.0636564D-03 -8.6231887D-03 8.6738253D-03 -4.9774245D-03 0 -5.0636564D-03 -8.6231887D-03 8.6738253D-03 -4.9774245D-03 0 -5.0636564D-03 -8.6231887D-03 8.6231887D-03 -5.0636564D-03 100 1.0880477D-02 -2.2983850D-02 2.2873004D-03 4.3590946D-03 100 1.0880477D-02 -2.2983850D-02 2.2873004D-03 4.3590946D-03 100 1.0880477D-02 -2.2983850D-02 2.3961052D-03 4.1292561D-03 J,Y Bessels X = 1.000 XLMIN = .00 KFN = 2 (CylBess) 0 7.6519769D-01 8.8256964D-02 -4.4005059D-01 7.8121282D-01 100 8.4318288-189 -3.7752878+185 8.4314114-187 3.7750971+187 J,Y Bessels X = 100.000 XLMIN = .00 KFN = 2 (CylBess) 0 1.9985850D-02 -7.7244313D-02 7.7145352D-02 2.0372312D-02 50 -3.8698340D-02 7.6505264D-02 -6.6001492D-02 -3.4025651D-02 COUL90 ERROR! IFAIL -1 SBESJY ERROR! IFAIL -1 RICBES ERROR! IFAIL -1 ETA = .000 X = -1.000 XLMIN = .00 KFN = 0 (Coulomb)