PDP-10 Archive: decus/20-0026/dqg8.doc from decus_20tap2

Trailing-Edge - PDP-10 Archives - decus_20tap2_198111 - decus/20-0026/dqg8.doc

There are 2 other files named dqg8.doc in the archive. Click here to see a list.

SUBROUTINE DQG8

PURPOSE
   TO COMPUTE INTEGRAL(FCT(X), SUMMED OVER X FROM XL TO XU)

USAGE
   CALL DQG8 (XL,XU,FCT,Y)
   PARAMETER FCT REQUIRES AN EXTERNAL STATEMENT

DESCRIPTION OF PARAMETERS
   XL	  - DOUBLE PRECISION LOWER BOUND OF THE INTERVAL.
   XU	  - DOUBLE PRECISION UPPER BOUND OF THE INTERVAL.
   FCT	  - THE NAME OF AN EXTERNAL DOUBLE PRECISION FUNCTION
	    SUBPROGRAM USED.
   Y	  - THE RESULTING DOUBLE PRECISION INTEGRAL VALUE.

REMARKS
   NONE

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
   THE EXTERNAL DOUBLE PRECISION FUNCTION SUBPROGRAM FCT(X)
   MUST BE FURNISHED BY THE USER.

METHOD
   EVALUATION IS DONE BY MEANS OF 8-POINT GAUSS QUADRATURE
   FORMULA, WHICH INTEGRATES POLYNOMIALS UP TO DEGREE 15
   EXACTLY. FOR REFERENCE, SEE
   V.I.KRYLOV, APPROXIMATE CALCULATION OF INTEGRALS,
   MACMILLAN, NEW YORK/LONDON, 1962, PP.100-111 AND 337-340.