Trailing-Edge
-
PDP-10 Archives
-
decuslib20-02
-
decus/20-0026/dqa8.doc
There are 2 other files named dqa8.doc in the archive. Click here to see a list.
SUBROUTINE DQA8
PURPOSE
TO COMPUTE INTEGRAL(EXP(-X)*FCT(X)/SQRT(X), SUMMED OVER X
FROM 0 TO INFINITY).
USAGE
CALL DQA8 (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 8-POINT GENERALIZED GAUSS-
LAGUERRE QUADRATURE FORMULA, WHICH INTEGRATES EXACTLY
WHENEVER FCT(X) IS A POLYNOMIAL UP TO DEGREE 15.
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.