FTN,B,L C C .................................................................. C C SUBROUTINE QTFG C C PURPOSE C TO COMPUTE THE VECTOR OF INTEGRAL VALUES FOR A GIVEN C GENERAL TABLE OF ARGUMENT AND FUNCTION VALUES. C C USAGE C CALL QTFG (X,Y,Z,NDIM) C C DESCRIPTION OF PARAMETERS C X - THE INPUT VECTOR OF ARGUMENT VALUES. C Y - THE INPUT VECTOR OF FUNCTION VALUES. C Z - THE RESULTING VECTOR OF INTEGRAL VALUES. Z MAY BE C IDENTICAL WITH X OR Y. C NDIM - THE DIMENSION OF VECTORS X,Y,Z. C C REMARKS C NO ACTION IN CASE NDIM LESS THAN 1. C C SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED C NONE C C METHOD C BEGINNING WITH Z(1)=0, EVALUATION OF VECTOR Z IS DONE BY C MEANS OF TRAPEZOIDAL RULE (SECOND ORDER FORMULA). C FOR REFERENCE, SEE C F.B.HILDEBRAND, INTRODUCTION TO NUMERICAL ANALYSIS, C MCGRAW-HILL, NEW YORK/TORONTO/LONDON, 1956, PP.75. C C .................................................................. C SUBROUTINE QTFG(X,Y,Z,NDIM) C C DIMENSION X(1),Y(1),Z(1) C SUM2=0. IF(NDIM-1)4,3,1 C C INTEGRATION LOOP 1 DO 2 I=2,NDIM SUM1=SUM2 SUM2=SUM2+.5*(X(I)-X(I-1))*(Y(I)+Y(I-1)) 2 Z(I-1)=SUM1 3 Z(NDIM)=SUM2 4 RETURN END END$ FTN,B,L C CALL--QTFG,QTFE,QSF,QHFG,QHFE DIMENSION X(316),Y(316),DRY(316),Z(316) 4 WRITE(2,5) 5 FORMAT("WRITE DIM,H") READ(1,*) DIM,H NDIM=DIM-1.0 NK=0 X(1)=0 Y(1)=0 DRY(1)=1 DO 7 J=1,NDIM N=J+1 X(N)=H*FLOAT(J) Y(N)=SIN(X(N)) 7 DRY(N)=COS(X(N)) NDIM=NDIM+1 CALL QTFG(X,Y,Z,NDIM) GO TO 30 30 CONTINUE WRITE(2,9)(Z(M),M=1,5) WRITE(2,9) Z(NDIM) PAUSE 1 16 GO TO 4 9 FORMAT(5E14.7) END END$