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

(Source file text) 

---

C	TABLE1.FOR - PROGRAM TO GENERATE TABLES OF:
C  GAMMA	- GAMMA FUNCTION
C  LEP		- LEGENDRE POLYNOMIALS
C  CS		- FRESNEL INTEGRALS
	DIMENSION A(10),POLY(11),B(10)
	DATA PI/3.1415926/
C
C	OUTPUT CHANNEL = IOUT
	IOUT=2
C	GAMMA FUNCTION
	X=.99
	H=.01
	XOLD=.9
	WRITE(IOUT,100)
100	FORMAT(/' GAMMA FUNCTION X=1(.01)2, 2(1)34'//'     X'/)
	DO 1 J=1,10
	XOLD=XOLD+.1
	DO 2 K=1,10
	X=X+H
2	CALL GAMMA(X,A(K),IER)
1	WRITE(IOUT,101) XOLD,A
101	FORMAT(F8.3/(5X,5E15.8))
	X=1.
	H=1.
	XOLD=-1.
	DO 3 J=1,11
	XOLD=XOLD+3.
	DO 4 I=1,3
	X=X+H
4	CALL GAMMA(X,A(I),IER)
3	WRITE(IOUT,101) XOLD,(A(K),K=1,3)
C
C	LEGENDRE POLYNOMIALS
	WRITE(IOUT,102)
102	FORMAT(//' LEGENDRE POLYNOMIALS X=0(0.05)1  N=0(1)10'//'     X'/)
	X=-0.05
	H=0.05
	DO 5 I=1,21
	X=X+H
	CALL LEP(POLY,X,10)
5	WRITE(IOUT,101) X,POLY
C
C	FRESNEL INTEGRALS
	WRITE(IOUT,103)
103	FORMAT(//' FRESNEL INTEGRALS U=(PI*X**2)/2  X=0.1(.1)5.0'/)
	PIBY2=PI/2.
	X=0.
	H=0.1
	XOLD=-.9
	DO 6 I=1,5
	XOLD=XOLD+1.
	DO 7 J=1,10
	X=X+H
	U=PIBY2*X**2
7	CALL CS(A(J),B(J),U)
	WRITE(IOUT,104) XOLD,A
104	FORMAT(F8.3,'   C(X)'/(5X5E15.8))
6	WRITE(IOUT,105) XOLD,B
105	FORMAT(F8.3,'   S(X)'/(5X5E15.8))
	STOP
	END

---