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

(Source file text) 

---

C	TABLE3.FTN - MATHEMATICAL TABLES USING:
C  CEL1		- ELIPTIC INTEGRALS OF FIRST KIND
C  CEL2		- ELIPTIC INTEGRALS OF SECOND KIND
C  EXPI		- EXPONENTIAL INTEGRAL
C  SICI		- SINE AND COSINE INTEGRALS
C
	DIMENSION V(10),VK(9),W(10)
C	OUTPUT CHANNEL = IOUT
	IOUT=2
C
C  CEL1
	X=0.00
	H=0.01
	XOLD=-0.08
	WRITE(IOUT,100)
100	FORMAT(/' ELLIPTIC INTEGRAL OF FIRST KIND:'/
	1'  K=SQRT(M)  M=.01(.01).99'//'    X'/)
	DO 1 I=1,11
	XOLD=XOLD+.09
	DO 2 J=1,9
	X=X+H
	XK=SQRT(X)
2	CALL CEL1(VK(J),XK,IER)
1	WRITE(IOUT,101) XOLD,VK
101	FORMAT(F6.2/(5X5E15.7))
C
C  CEL2
	WRITE(IOUT,102)
102	FORMAT(//' ELLIPTIC INTEGRAL OF SECOND KIND:'/
	1' CASE 1:  A = B = 1, K=SQRT(M)  M=.01(.01).99'//'    X'/)
	X=0.0
	H=0.01
	XOLD=-0.08
	DO 3 I=1,11
	XOLD=XOLD+.09
	DO 4 J=1,9
	X=X+H
	XK=SQRT(X)
4	CALL CEL2(VK(J),XK,1.,1.,IER)
3	WRITE(IOUT,101) XOLD,VK
	WRITE(IOUT,103)
103	FORMAT(/' CASE 2:  A=1, B=1-K**2, K=SQRT(M)  M=.01(.01).99'
	1//'    X'/)
	X=0.0
	H=0.01
	XOLD=-0.08
	DO 5 I=1,11
	XOLD=XOLD+.09
	DO 6 J=1,9
	X=X+H
	XK=SQRT(X)
	B=1.-X
6	CALL CEL2(VK(J),XK,1.,B,IER)
5	WRITE(IOUT,101) XOLD,VK
C
C  SICI
	WRITE(IOUT,104)
104	FORMAT(//' SINE AND COSINE INTEGRALS:  X=0.1(0.1)10.0'//'    X'/)
	X=0.0
	H=0.1
	XOLD=-.9
	DO 7 I=1,10
	XOLD=XOLD+1.
	DO 8 J=1,10
	X=X+H
8	CALL SICI(V(J),W(J),X)
	WRITE( IOUT,105) XOLD,V
105	FORMAT(/F8.3,5X,'SI(X)-PI/2'/(5X5E15.7))
7	WRITE(IOUT,106) XOLD,W
106	FORMAT(/F8.3,5X,'CI(X)'/(5X5E15.7))
C
C  EXPI
	WRITE(IOUT,107)
107	FORMAT(//' EXPONENTIAL INTEGRAL:  X=-10.(.2)6.'//'    X'/)
	X=-10.2
	H=0.2
	XOLD=-11.
	DO 9 I=1,16
	XOLD=XOLD+1.
	DO 10 J=1,5
	X=X+H
	IF (ABS(X).LT.1.E-4)  X=0.E0
10	CALL EXPI(X,V(J),AUX)
9	WRITE(IOUT,101) XOLD,(V(K),K=1,5)
	STOP
	END

---