Google
 

Trailing-Edge - PDP-10 Archives - decuslib20-02 - decus/20-0026/qa8.ssp
There are 2 other files named qa8.ssp in the archive. Click here to see a list.
C                                                                       QA8   10
C     ..................................................................QA8   20
C                                                                       QA8   30
C        SUBROUTINE QA8                                                 QA8   40
C                                                                       QA8   50
C        PURPOSE                                                        QA8   60
C           TO COMPUTE INTEGRAL(EXP(-X)*FCT(X)/SQRT(X), SUMMED OVER X   QA8   70
C                               FROM 0 TO INFINITY).                    QA8   80
C                                                                       QA8   90
C        USAGE                                                          QA8  100
C           CALL QA8 (FCT,Y)                                            QA8  110
C           PARAMETER FCT REQUIRES AN EXTERNAL STATEMENT                QA8  120
C                                                                       QA8  130
C        DESCRIPTION OF PARAMETERS                                      QA8  140
C           FCT    - THE NAME OF AN EXTERNAL FUNCTION SUBPROGRAM USED.  QA8  150
C           Y      - THE RESULTING INTEGRAL VALUE.                      QA8  160
C                                                                       QA8  170
C        REMARKS                                                        QA8  180
C           NONE                                                        QA8  190
C                                                                       QA8  200
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  QA8  210
C           THE EXTERNAL FUNCTION SUBPROGRAM FCT(X) MUST BE FURNISHED   QA8  220
C           BY THE USER.                                                QA8  230
C                                                                       QA8  240
C        METHOD                                                         QA8  250
C           EVALUATION IS DONE BY MEANS OF 8-POINT GENERALIZED GAUSSIAN-QA8  260
C           LAGUERRE QUADRATURE FORMULA, WHICH INTEGRATES EXACTLY,      QA8  270
C           WHENEVER FCT(X) IS A POLYNOMIAL UP TO DEGREE 15.            QA8  280
C           FOR REFERENCE, SEE                                          QA8  290
C           CONCUS/CASSATT/JAEHNIG/MELBY, TABLES FOR THE EVALUATION OF  QA8  300
C           INTEGRAL(X**BETA*EXP(-X)*F(X), SUMMED OVER X FROM 0 TO      QA8  310
C           INFINITY) BY GAUSS-LAGUERRE QUADRATURE, MTAC, VOL.17,       QA8  320
C           ISS.83 (1963), PP.245-256.                                  QA8  330
C                                                                       QA8  340
C     ..................................................................QA8  350
C                                                                       QA8  360
      SUBROUTINE QA8(FCT,Y)                                             QA8  370
C                                                                       QA8  380
C                                                                       QA8  390
      X=21.98427                                                        QA8  400
      Y=.5309615E-9*FCT(X)                                              QA8  410
      X=14.97262                                                        QA8  420
      Y=Y+.4641962E-6*FCT(X)                                            QA8  430
      X=10.09332                                                        QA8  440
      Y=Y+.5423720E-4*FCT(X)                                            QA8  450
      X=6.483145                                                        QA8  460
      Y=Y+.001864568*FCT(X)                                             QA8  470
      X=3.809476                                                        QA8  480
      Y=Y+.02576062*FCT(X)                                              QA8  490
      X=1.905114                                                        QA8  500
      Y=Y+.1676201*FCT(X)                                               QA8  510
      X=.6772491                                                        QA8  520
      Y=Y+.5612949*FCT(X)                                               QA8  530
      X=.07479188                                                       QA8  540
      Y=Y+1.015859*FCT(X)                                               QA8  550
      RETURN                                                            QA8  560
      END                                                               QA8  570