Google
 

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