Google
 

Trailing-Edge - PDP-10 Archives - decuslib20-02 - decus/20-0026/expon.cdk
There are 2 other files named expon.cdk in the archive. Click here to see a list.
$JOB EXPON[30,30]
$FORTRAN EXPON
C                                                                       EXPO  10
C     ..................................................................EXPO  20
C                                                                       EXPO  30
C        SAMPLE MAIN PROGRAM FOR TRIPLE EXPONENTIAL SMOOTHING - EXPON   EXPO  40
C                                                                       EXPO  50
C        PURPOSE                                                        EXPO  60
C           (1) READ THE PROBLEM PARAMETER CARD AND A TIME SERIES,      EXPO  70
C           (2) CALL THE SUBROUTINE EXSMO TO SMOOTH THE TIME SERIES,    EXPO  80
C           AND (3) PRINT THE RESULT.                                   EXPO  90
C                                                                       EXPO 100
C        REMARKS                                                        EXPO 110
C           A SMOOTHING CONSTANT SPECIFIED IN THE PROBLEM PARAMETER     EXPO 120
C           CARD MUST BE GREATER THAN ZERO BUT LESS THAN ONE IN ORDER   EXPO 130
C           TO OBTAIN REASONABLE RESULTS.                               EXPO 140
C                                                                       EXPO 150
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  EXPO 160
C           EXSMO                                                       EXPO 170
C                                                                       EXPO 180
C        METHOD                                                         EXPO 190
C           REFER TO R. G. BROWN, 'SMOOTHING, FORECASTING AND PREDICTIONEXPO 200
C           OF DISCRETE TIME SERIES', PRENTICE-HALL, N.J., 1963,        EXPO 210
C           PP. 140 TO 144.                                             EXPO 220
C                                                                       EXPO 230
C     ..................................................................EXPO 240
C                                                                       EXPO 250
C     THE FOLLOWING DIMENSION MUST BE GREATER THAN OR EQUAL TO THE      EXPO 260
C     NUMBER OF DATA POINTS IN A GIVEN TIME SERIES..                    EXPO 270
C                                                                       EXPO 280
         DIMENSION X(1000),S(1000)                                      EXPO 290
C                                                                       EXPO 300
C     ..................................................................EXPO 310
C                                                                       EXPO 320
    1 FORMAT(A4,A2,I4,F5.0,3F10.0)                                      EXPO 330
    2 FORMAT(12F6.0)                                                    EXPO 340
    3 FORMAT(34H1TRIPLE EXPONENTIAL SMOOTHING.....,A4,A2//22H NUMBER OF EXPO 350
     1DATA POINTS,I6/19H SMOOTHING CONSTANT,F9.3/)                      EXPO 360
    4 FORMAT(13H0COEFFICIENTS,9X,1HA,14X,1HB,14X,1HC)                   EXPO 370
    5 FORMAT(9H0ORIGINAL,F19.5,2F15.5)                                  EXPO 380
    6 FORMAT(8H0UPDATED,F20.5,2F15.5/)                                  EXPO 390
    7 FORMAT(1H0,27X,13HSMOOTHED DATA/7X,10HINPUT DATA,12X,10H(FORECAST)EXPO 400
     1)                                                                 EXPO 410
    8 FORMAT(F17.5,8X,F15.5)                                            EXPO 420
C                                                                       EXPO 430
C     ..................................................................EXPO 440
C                                                                       EXPO 450
C     READ PROBLEM PARAMETER CARD                                       EXPO 460
C                                                                       EXPO 470
  100 READ (5,1,END=999) PR,PR1,NX,AL,A,B,C                             EXPO 480
C        PR......PROBLEM NUMBER (MAY BE ALPHAMERIC)                     EXPO 490
C        PR1.....PROBLEM NUMBER (CONTINUED)                             EXPO 500
C        NX......NUMBER OF DATA POINTS IN TIME SERIES                   EXPO 510
C        AL......SMOOTHING CONSTANT                                     EXPO 520
C        A,B,C...COEFFICIENTS OF THE PREDICTION EQUATION                EXPO 530
C                                                                       EXPO 540
      WRITE (6,3) PR,PR1,NX,AL                                          EXPO 550
C                                                                       EXPO 560
C     PRINT ORIGINAL COEFFICIENTS                                       EXPO 570
C                                                                       EXPO 580
      WRITE (6,4)                                                       EXPO 590
      WRITE (6,5) A,B,C                                                 EXPO 600
C                                                                       EXPO 610
C     READ TIME SERIES DATA                                             EXPO 620
C                                                                       EXPO 630
      READ (5,2) (X(I),I=1,NX)                                          EXPO 640
C                                                                       EXPO 650
      CALL EXSMO (X,NX,AL,A,B,C,S)                                      EXPO 660
C                                                                       EXPO 670
C     PRINT UPDATED COEFFICIENTS                                        EXPO 680
C                                                                       EXPO 690
      WRITE (6,6) A,B,C                                                 EXPO 700
C                                                                       EXPO 710
C     PRINT INPUT AND SMOOTHED DATA                                     EXPO 720
C                                                                       EXPO 730
      WRITE (6,7)                                                       EXPO 740
      DO 200 I=1,NX                                                     EXPO 750
  200 WRITE (6,8) X(I),S(I)                                             EXPO 760
      GO TO 100                                                         EXPO 770
999	STOP
      END                                                               EXPO 780
$DECK EXP.CDR
SAMPLE  38  0.1       0.0       0.0       0.0                                 20
   430   426   422   419   414   413   412   409   411   417   422   430      30
   438   441   447   455   461   453   448   449   454   463   470   472      40
   476   481   483   487   491   492   485   486   482   479   479   476      50
   472   470                                                                  60
$EOD
.ASSIGN CDR 5
.ASSIGN LPT 6
.SET CDR EXP
.EXECUTE/REL EXPON,WES:SSP/LIB
%FIN::
.DELETE EXP.CDR