Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-02 - 43,50145/pint.ssp
There are 2 other files named pint.ssp in the archive. Click here to see a list.
C                                                                       PINT  10
C     ..................................................................PINT  20
C                                                                       PINT  30
C        SUBROUTINE PINT                                                PINT  40
C                                                                       PINT  50
C        PURPOSE                                                        PINT  60
C           FIND INTEGRAL OF A POLYNOMIAL WITH CONSTANT OF INTEGRATION  PINT  70
C           EQUAL TO ZERO                                               PINT  80
C                                                                       PINT  90
C        USAGE                                                          PINT 100
C           CALL PINT(Y,IDIMY,X,IDIMX)                                  PINT 110
C                                                                       PINT 120
C        DESCRIPTION OF PARAMETERS                                      PINT 130
C           Y     - VECTOR OF COEFFICIENTS FOR INTEGRAL, ORDERED FROM   PINT 140
C                   SMALLEST TO LARGEST POWER                           PINT 150
C           IDIMY - DIMENSION OF Y (EQUAL TO IDIMX+1)                   PINT 160
C           X     - VECTOR OF COEFFICIENTS FOR ORIGINAL POLYNOMIAL,     PINT 170
C                   ORDERED FROM SMALLEST TO LARGEST POWER              PINT 180
C           IDIMX - DIMENSION OF X                                      PINT 190
C                                                                       PINT 200
C        REMARKS                                                        PINT 210
C           NONE                                                        PINT 220
C                                                                       PINT 230
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  PINT 240
C           NONE                                                        PINT 250
C                                                                       PINT 260
C        METHOD                                                         PINT 270
C           DIMENSION OF Y IS SET AT DIMENSION OF X PLUS ONE, AND THE   PINT 280
C           CONSTANT TERM IS SET TO ZERO. INTEGRAL IS THEN CALCULATED   PINT 290
C           BY DIVIDING COEFFICIENTS BY THEIR RESPECTIVE EXPONENTS.     PINT 300
C                                                                       PINT 310
C     ..................................................................PINT 320
C                                                                       PINT 330
      SUBROUTINE PINT(Y,IDIMY,X,IDIMX)                                  PINT 340
      DIMENSION X(1),Y(1)                                               PINT 350
C                                                                       PINT 360
      IDIMY=IDIMX+1                                                     PINT 370
      Y(1)=0.                                                           PINT 380
      IF(IDIMX)1,1,2                                                    PINT 390
    1 RETURN                                                            PINT 400
    2 EXPT=1.                                                           PINT 410
      DO 3 I=2,IDIMY                                                    PINT 420
      Y(I)=X(I-1)/EXPT                                                  PINT 430
    3 EXPT=EXPT+1.                                                      PINT 440
      GO TO 1                                                           PINT 450
      END                                                               PINT 460