PDP-10 Archive: decus/20-0026/dqa24.doc from decuslib20-02

Trailing-Edge - PDP-10 Archives - decuslib20-02 - decus/20-0026/dqa24.doc

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

SUBROUTINE DQA24

PURPOSE
   TO COMPUTE INTEGRAL(EXP(-X)*FCT(X)/SQRT(X), SUMMED OVER X
		       FROM 0 TO INFINITY).

USAGE
   CALL DQA24 (FCT,Y)
   PARAMETER FCT REQUIRES AN EXTERNAL STATEMENT

DESCRIPTION OF PARAMETERS
   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 24-POINT GENERALIZED GAUSS-
   LAGUERRE QUADRATURE FORMULA, WHICH INTEGRATES EXACTLY
   WHENEVER FCT(X) IS A POLYNOMIAL UP TO DEGREE 47.
   FOR REFERENCE, SEE
   SHAO/CHEN/FRANK, TABLES OF ZEROS AND GAUSSIAN WEIGHTS OF
   CERTAIN ASSOCIATED LAGUERRE POLYNOMIALS AND THE RELATED
   GENERALIZED HERMITE POLYNOMIALS, IBM TECHNICAL REPORT
   TR00.1100 (MARCH 1964), PP.15-16.