Google
 

Trailing-Edge - PDP-10 Archives - decus_20tap2_198111 - decus/20-0026/ql10.ssp
There are 2 other files named ql10.ssp in the archive. Click here to see a list.
C                                                                       QL10  10
C     ..................................................................QL10  20
C                                                                       QL10  30
C        SUBROUTINE QL10                                                QL10  40
C                                                                       QL10  50
C        PURPOSE                                                        QL10  60
C           TO COMPUTE INTEGRAL(EXP(-X)*FCT(X), SUMMED OVER X FROM 0    QL10  70
C                               TO INFINITY).                           QL10  80
C                                                                       QL10  90
C        USAGE                                                          QL10 100
C           CALL QL10(FCT,Y)                                            QL10 110
C           PARAMETER FCT REQUIRES AN EXTERNAL STATEMENT                QL10 120
C                                                                       QL10 130
C        DESCRIPTION OF PARAMETERS                                      QL10 140
C           FCT    - THE NAME OF AN EXTERNAL FUNCTION SUBPROGRAM USED.  QL10 150
C           Y      - THE RESULTING INTEGRAL VALUE.                      QL10 160
C                                                                       QL10 170
C        REMARKS                                                        QL10 180
C           NONE                                                        QL10 190
C                                                                       QL10 200
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  QL10 210
C           THE EXTERNAL FUNCTION SUBPROGRAM FCT(X) MUST BE FURNISHED   QL10 220
C           BY THE USER.                                                QL10 230
C                                                                       QL10 240
C        METHOD                                                         QL10 250
C           EVALUATION IS DONE BY MEANS OF 10-POINT GAUSSIAN-LAGUERRE   QL10 260
C           QUADRATURE FORMULA, WHICH INTEGRATES EXACTLY WHENEVER       QL10 270
C           FCT(X) IS A POLYNOMIAL UP TO DEGREE 19.                     QL10 280
C           FOR REFERENCE, SEE                                          QL10 290
C           V.I.KRYLOV, APPROXIMATE CALCULATION OF INTEGRALS,           QL10 300
C           MACMILLAN, NEW YORK/LONDON, 1962, PP.130-132 AND 347-352.   QL10 310
C                                                                       QL10 320
C     ..................................................................QL10 330
C                                                                       QL10 340
      SUBROUTINE QL10(FCT,Y)                                            QL10 350
C                                                                       QL10 360
C                                                                       QL10 370
      X=29.92070                                                        QL10 380
      Y=.9911827E-12*FCT(X)                                             QL10 390
      X=21.99659                                                        QL10 400
      Y=Y+.1839565E-8*FCT(X)                                            QL10 410
      X=16.27926                                                        QL10 420
      Y=Y+.4249314E-6*FCT(X)                                            QL10 430
      X=11.84379                                                        QL10 440
      Y=Y+.2825923E-4*FCT(X)                                            QL10 450
      X=8.330153                                                        QL10 460
      Y=Y+.7530084E-3*FCT(X)                                            QL10 470
      X=5.552496                                                        QL10 480
      Y=Y+.009501517*FCT(X)                                             QL10 490
      X=3.401434                                                        QL10 500
      Y=Y+.06208746*FCT(X)                                              QL10 510
      X=1.808343                                                        QL10 520
      Y=Y+.2180683*FCT(X)                                               QL10 530
      X=.7294545                                                        QL10 540
      Y=Y+.4011199*FCT(X)                                               QL10 550
      X=.1377935                                                        QL10 560
      Y=Y+.3084411*FCT(X)                                               QL10 570
      RETURN                                                            QL10 580
      END                                                               QL10 590