File: SPHBES.FC of Tape: Various/ETH/fc1

(.FC,.FO Focal format converted to 'WRITE' listing) 

---

C-SPHBES.FC OCT-75

15.10 A ?L,X?
15.20 D 16
15.30 T %,AJ,!
15.40 G 15.1

16.01 C SPHERICAL BESSEL FUNC: J(L,X). L=-8 TO 8, X ANY POSITIVE NUMBER, RESULT = AJ
16.10 S X=FABS(X);S L=FITR(L)
16.20 S BJ(1)=1.04*X+6.5
16.22 S BJ(2)=1.8284*FLOG(X)+3
16.25 S BJ(3)=2
16.30 I (BJ(2),>,2)16.35;S BJ(3)=BJ(2)
16.35 S NS=L+BJ(3)
16.40 I (BJ(1),>,L)16.45;S NS=BJ(1)+BJ(3)
16.45 S M=FITR(NS+1)
16.50 S BJ(3)=0;S BJ(2)=1E-9
16.60 S T=3;S NS=FSIN(X);I (NS,<,.707107)16.65;S NS=FCOS(X);S T=2
16.65 F M=M,-1,1;D 16.95;S BJ(3)=BJ(2);S BJ(2)=BJ(1)
16.70 S NS=NS/BJ(T)
16.75 F M=0,-1,L+2;D 16.95;S BJ(3)=BJ(2);S BJ(2)=BJ(1)
16.80 S AJ=NS*1.0000002725*BL/X
16.85 R
16.95 S BJ(1)=(2*M-1)*BJ(2)/X-BJ(3);I (FITR(M),=,FITR(L+2))16.99;S BL=BJ(1)
16.99 C

---