Google
 

Trailing-Edge - PDP-10 Archives - decuslib20-02 - decus/20-0026/qa9.ssp
There are 2 other files named qa9.ssp in the archive. Click here to see a list.
C                                                                       QA9   10
C     ..................................................................QA9   20
C                                                                       QA9   30
C        SUBROUTINE QA9                                                 QA9   40
C                                                                       QA9   50
C        PURPOSE                                                        QA9   60
C           TO COMPUTE INTEGRAL(EXP(-X)*FCT(X)/SQRT(X), SUMMED OVER X   QA9   70
C                               FROM 0 TO INFINITY).                    QA9   80
C                                                                       QA9   90
C        USAGE                                                          QA9  100
C           CALL QA9 (FCT,Y)                                            QA9  110
C           PARAMETER FCT REQUIRES AN EXTERNAL STATEMENT                QA9  120
C                                                                       QA9  130
C        DESCRIPTION OF PARAMETERS                                      QA9  140
C           FCT    - THE NAME OF AN EXTERNAL FUNCTION SUBPROGRAM USED.  QA9  150
C           Y      - THE RESULTING INTEGRAL VALUE.                      QA9  160
C                                                                       QA9  170
C        REMARKS                                                        QA9  180
C           NONE                                                        QA9  190
C                                                                       QA9  200
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  QA9  210
C           THE EXTERNAL FUNCTION SUBPROGRAM FCT(X) MUST BE FURNISHED   QA9  220
C           BY THE USER.                                                QA9  230
C                                                                       QA9  240
C        METHOD                                                         QA9  250
C           EVALUATION IS DONE BY MEANS OF 9-POINT GENERALIZED GAUSSIAN-QA9  260
C           LAGUERRE QUADRATURE FORMULA, WHICH INTEGRATES EXACTLY,      QA9  270
C           WHENEVER FCT(X) IS A POLYNOMIAL UP TO DEGREE 17.            QA9  280
C           FOR REFERENCE, SEE                                          QA9  290
C           CONCUS/CASSATT/JAEHNIG/MELBY, TABLES FOR THE EVALUATION OF  QA9  300
C           INTEGRAL(X**BETA*EXP(-X)*F(X), SUMMED OVER X FROM 0 TO      QA9  310
C           INFINITY) BY GAUSS-LAGUERRE QUADRATURE, MTAC, VOL.17,       QA9  320
C           ISS.83 (1963), PP.245-256.                                  QA9  330
C                                                                       QA9  340
C     ..................................................................QA9  350
C                                                                       QA9  360
      SUBROUTINE QA9(FCT,Y)                                             QA9  370
C                                                                       QA9  380
C                                                                       QA9  390
      X=25.48598                                                        QA9  400
      Y=.1565640E-10*FCT(X)                                             QA9  410
      X=18.04651                                                        QA9  420
      Y=Y+.2093441E-7*FCT(X)                                            QA9  430
      X=12.77183                                                        QA9  440
      Y=Y+.3621309E-5*FCT(X)                                            QA9  450
      X=8.769757                                                        QA9  460
      Y=Y+.0001836225*FCT(X)                                            QA9  470
      X=5.694423                                                        QA9  480
      Y=Y+.003777045*FCT(X)                                             QA9  490
      X=3.369176                                                        QA9  500
      Y=Y+.03728008*FCT(X)                                              QA9  510
      X=1.692395                                                        QA9  520
      Y=Y+.1946035*FCT(X)                                               QA9  530
      X=.6032364                                                        QA9  540
      Y=Y+.5696146*FCT(X)                                               QA9  550
      X=.06670223                                                       QA9  560
      Y=Y+.9669914*FCT(X)                                               QA9  570
      RETURN                                                            QA9  580
      END                                                               QA9  590