Trailing-Edge
-
PDP-10 Archives
-
decuslib20-02
-
decus/20-0026/qa10.ssp
There are 2 other files named qa10.ssp in the archive. Click here to see a list.
C QA10 10
C ..................................................................QA10 20
C QA10 30
C SUBROUTINE QA10 QA10 40
C QA10 50
C PURPOSE QA10 60
C TO COMPUTE INTEGRAL(EXP(-X)*FCT(X)/SQRT(X), SUMMED OVER X QA10 70
C FROM 0 TO INFINITY). QA10 80
C QA10 90
C USAGE QA10 100
C CALL QA10 (FCT,Y) QA10 110
C PARAMETER FCT REQUIRES AN EXTERNAL STATEMENT QA10 120
C QA10 130
C DESCRIPTION OF PARAMETERS QA10 140
C FCT - THE NAME OF AN EXTERNAL FUNCTION SUBPROGRAM USED. QA10 150
C Y - THE RESULTING INTEGRAL VALUE. QA10 160
C QA10 170
C REMARKS QA10 180
C NONE QA10 190
C QA10 200
C SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED QA10 210
C THE EXTERNAL FUNCTION SUBPROGRAM FCT(X) MUST BE FURNISHED QA10 220
C BY THE USER. QA10 230
C QA10 240
C METHOD QA10 250
C EVALUATION IS DONE BY MEANS OF 10-POINT GENERALIZED GAUSS- QA10 260
C LAGUERRE QUADRATURE FORMULA, WHICH INTEGRATES EXACTLY, QA10 270
C WHENEVER FCT(X) IS A POLYNOMIAL UP TO DEGREE 19. QA10 280
C FOR REFERENCE, SEE QA10 290
C CONCUS/CASSATT/JAEHNIG/MELBY, TABLES FOR THE EVALUATION OF QA10 300
C INTEGRAL(X**BETA*EXP(-X)*F(X), SUMMED OVER X FROM 0 TO QA10 310
C INFINITY) BY GAUSS-LAGUERRE QUADRATURE, MTAC, VOL.17, QA10 320
C ISS.83 (1963), PP.245-256. QA10 330
C QA10 340
C ..................................................................QA10 350
C QA10 360
SUBROUTINE QA10(FCT,Y) QA10 370
C QA10 380
C QA10 390
X=29.02495 QA10 400
Y=.4458787E-12*FCT(X) QA10 410
X=21.19389 QA10 420
Y=Y+.8798682E-9*FCT(X) QA10 430
X=15.56116 QA10 440
Y=Y+.2172139E-6*FCT(X) QA10 450
X=11.20813 QA10 460
Y=Y+.1560511E-4*FCT(X) QA10 470
X=7.777439 QA10 480
Y=Y+.0004566773*FCT(X) QA10 490
X=5.084908 QA10 500
Y=Y+.006487547*FCT(X) QA10 510
X=3.022513 QA10 520
Y=Y+.04962104*FCT(X) QA10 530
X=1.522944 QA10 540
Y=Y+.2180344*FCT(X) QA10 550
X=.5438675 QA10 560
Y=Y+.5733510*FCT(X) QA10 570
X=.06019206 QA10 580
Y=Y+.9244873*FCT(X) QA10 590
RETURN QA10 600
END QA10 610