Google
 

Trailing-Edge - PDP-10 Archives - decuslib20-02 - decus/20-0026/expon.smp
There are 2 other files named expon.smp in the archive. Click here to see a list.
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