Google
 

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