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

(Source file text) 

---

C  QDINT.FTN - SAMPLE PROGRAM FOR INTEGRATION OF A TABULATED FUNCTION BY
C	NUMERICAL QUADRATURE USING THE SUBROUTINE  QSF
C	THE FOLLOWING DIMENSION MUST BE AS LARGE AS THE MAXIMUM NUMBER
C	OF TABULATED VALUES TO BE INTEGRATED
	DIMENSION Z(500)
10	FORMAT(2I5,F10.0)
20	FORMAT(////' INTEGRATION OF TABULATED VALUES FOR DY/DX USING SUB
	1ROUTINE QSF'///11H FUNCTION  ,I5,3X,I5,17H TABULATED VALUES,
	25X,10HINTERVAL =,E15.8/)
22	FORMAT(/18H ILLEGAL CONDITION/)
23	FORMAT(/46H NUMBER OF TABULATED VALUES IS LESS THAN THREE)
24	FORMAT(/27H WPECIFIED INTERVAL IS ZERO)
30	FORMAT(/7X,'RESULTANT VALUE OF INTEGRAL AT EACH STEP IS',/
	1(1H ,5E15.8))
32	FORMAT(7F10.0)
C
C	OUTPUT CHANNEL = IOUT, INPUT CHANNEL = IN
	IOUT = 2
	IN = 1
C
35	READ(IN,10) ICOD,NUMBR,SPACE
	IF(ICOD+NUMBR) 70,70,38
70	STOP
38	WRITE(IOUT,20) ICOD,NUMBR,SPACE
	IF(NUMBR-3) 100,50,50
50	READ(IN,32)(Z(I),I=1,NUMBR)
	CALL QSF(SPACE,Z,Z,NUMBR)
	IF(SPACE) 60,200,60
60	WRITE(IOUT,30) (Z(I),I=1,NUMBR)
	GO TO 35
100	WRITE(IOUT,22)
	WRITE(IOUT,23)
	READ(IN,32) (Z(I),I=1,NUMBR)
	GO TO 35
200	WRITE(IOUT,22)
	WRITE(IOUT,24)
	GO TO 35
	END

---