Trailing-Edge
-
PDP-10 Archives
-
decuslib10-02
-
43,50145/ql6.doc
There are 2 other files named ql6.doc in the archive. Click here to see a list.
SUBROUTINE QL6
PURPOSE
TO COMPUTE INTEGRAL(EXP(-X)*FCT(X), SUMMED OVER X FROM 0
TO INFINITY).
USAGE
CALL QL6 (FCT,Y)
PARAMETER FCT REQUIRES AN EXTERNAL STATEMENT
DESCRIPTION OF PARAMETERS
FCT - THE NAME OF AN EXTERNAL FUNCTION SUBPROGRAM USED.
Y - THE RESULTING INTEGRAL VALUE.
REMARKS
NONE
SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
THE EXTERNAL FUNCTION SUBPROGRAM FCT(X) MUST BE FURNISHED
BY THE USER.
METHOD
EVALUATION IS DONE BY MEANS OF 6-POINT GAUSSIAN-LAGUERRE
QUADRATURE FORMULA, WHICH INTEGRATES EXACTLY WHENEVER
FCT(X) IS A POLYNOMIAL UP TO DEGREE 11.
FOR REFERENCE, SEE
V.I.KRYLOV, APPROXIMATE CALCULATION OF INTEGRALS,
MACMILLAN, NEW YORK/LONDON, 1962, PP.130-132 AND 347-352.