Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-02 - 43,50145/qa10.doc
There are 2 other files named qa10.doc in the archive. Click here to see a list.
SUBROUTINE QA10

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

USAGE
   CALL QA10 (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 10-POINT GENERALIZED GAUSS-
   LAGUERRE QUADRATURE FORMULA, WHICH INTEGRATES EXACTLY,
   WHENEVER FCT(X) IS A POLYNOMIAL UP TO DEGREE 19.
   FOR REFERENCE, SEE
   CONCUS/CASSATT/JAEHNIG/MELBY, TABLES FOR THE EVALUATION OF
   INTEGRAL(X**BETA*EXP(-X)*F(X), SUMMED OVER X FROM 0 TO
   INFINITY) BY GAUSS-LAGUERRE QUADRATURE, MTAC, VOL.17,
   ISS.83 (1963), PP.245-256.