Google
 

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