Google
 

Trailing-Edge - PDP-10 Archives - decus_20tap2_198111 - decus/20-0026/ql7.ssp
There are 2 other files named ql7.ssp in the archive. Click here to see a list.
C                                                                       QL7   10
C     ..................................................................QL7   20
C                                                                       QL7   30
C        SUBROUTINE QL7                                                 QL7   40
C                                                                       QL7   50
C        PURPOSE                                                        QL7   60
C           TO COMPUTE INTEGRAL(EXP(-X)*FCT(X), SUMMED OVER X FROM 0    QL7   70
C                               TO INFINITY).                           QL7   80
C                                                                       QL7   90
C        USAGE                                                          QL7  100
C           CALL QL7 (FCT,Y)                                            QL7  110
C           PARAMETER FCT REQUIRES AN EXTERNAL STATEMENT                QL7  120
C                                                                       QL7  130
C        DESCRIPTION OF PARAMETERS                                      QL7  140
C           FCT    - THE NAME OF AN EXTERNAL FUNCTION SUBPROGRAM USED.  QL7  150
C           Y      - THE RESULTING INTEGRAL VALUE.                      QL7  160
C                                                                       QL7  170
C        REMARKS                                                        QL7  180
C           NONE                                                        QL7  190
C                                                                       QL7  200
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  QL7  210
C           THE EXTERNAL FUNCTION SUBPROGRAM FCT(X) MUST BE FURNISHED   QL7  220
C           BY THE USER.                                                QL7  230
C                                                                       QL7  240
C        METHOD                                                         QL7  250
C           EVALUATION IS DONE BY MEANS OF 7-POINT GAUSSIAN-LAGUERRE    QL7  260
C           QUADRATURE FORMULA, WHICH INTEGRATES EXACTLY WHENEVER       QL7  270
C           FCT(X) IS A POLYNOMIAL UP TO DEGREE 13.                     QL7  280
C           FOR REFERENCE, SEE                                          QL7  290
C           V.I.KRYLOV, APPROXIMATE CALCULATION OF INTEGRALS,           QL7  300
C           MACMILLAN, NEW YORK/LONDON, 1962, PP.130-132 AND 347-352.   QL7  310
C                                                                       QL7  320
C     ..................................................................QL7  330
C                                                                       QL7  340
      SUBROUTINE QL7(FCT,Y)                                             QL7  350
C                                                                       QL7  360
C                                                                       QL7  370
      X=19.39573                                                        QL7  380
      Y=.3170315E-7*FCT(X)                                              QL7  390
      X=12.73418                                                        QL7  400
      Y=Y+.1586546E-4*FCT(X)                                            QL7  410
      X=8.182153                                                        QL7  420
      Y=Y+.1074010E-2*FCT(X)                                            QL7  430
      X=4.900353                                                        QL7  440
      Y=Y+.02063351*FCT(X)                                              QL7  450
      X=2.567877                                                        QL7  460
      Y=Y+.1471263*FCT(X)                                               QL7  470
      X=1.026665                                                        QL7  480
      Y=Y+.4218313*FCT(X)                                               QL7  490
      X=.1930437                                                        QL7  500
      Y=Y+.4093190*FCT(X)                                               QL7  510
      RETURN                                                            QL7  520
      END                                                               QL7  530