File: TABLE2.FT of Tape: Various/ETH/eth11-1

(Source file text) 

---

C	TABLE2.FTN - MATHEMATICAL TABLES USING:
C  BESJ		- J BESSEL FUNCTION
C  BESY		- Y BESSEL FUNCTION
C  BESI		- I BESSEL FUNCTION
C  BESK		- K BESSEL FUNCTION
C
C	RANGE  N=0(1)10  X=0(1)20
C
	DIMENSION BJ(11),BY(11),BI(11),BK(11),IM1(11)
	DIMENSION IJ(11),IY(11),II(11),IK(11)
C
C	INITIALIZE ARRAY FOR OUTPUT
	DO 1 I=1,11
1	IM1(I)=I-1
C
C	OUTPUT CHANNEL = IOUT
	IOUT=2
	X=-1.
	DX=1.
	WRITE(IOUT,100)
100	FORMAT(/' BESSEL FUNCTIONS N=0(1)10  X=0(1)20'//
	1'   X   N'4X'BESJ'7X'IER'3X'BESY'7X'IER'3X'BESI'7X'IER'3X'BESK'
	27X'IER'//)
	DO 2 I=1,21
	X=X+DX
	DO 3 J=1,11
	JM1=J-1
	CALL BESJ(X,JM1,BJ(J),1.E-5,IJ(J))
	CALL BESY(X,JM1,BY(J),IY(J))
	CALL BESI(X,JM1,BI(J),II(J))
3	CALL BESK(X,JM1,BK(J),IK(J))
2	WRITE(IOUT,101) X,(IM1(K),BJ(K),IJ(K),BY(K),IY(K),BI(K),II(K),
	1BK(K),IK(K),K=1,11)
101	FORMAT(F5.1/(5X,I3,1X,E12.5,I4,1X,E12.5,I4,1X,E12.5,I4,1X,
	1E12.5,I4))
	STOP
	END

---