home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l292 / 1.ddi / SMOOTHWE.FOR < prev    next >
Encoding:
Text File  |  1989-10-10  |  1.3 KB  |  55 lines
  1. !!      PROGRAM Smoothwe
  2.  
  3.       INCLUDE 'STDHDR.FOR'
  4.       REAL xdata(0:maxv), ydata(0:maxv)
  5.       REAL smoothdata(0:maxv), weights(0:maxc)
  6.       INTEGER smoothnum,n,i
  7.       REAL wdivisor, rr
  8.       CHARACTER * 80 pathname
  9.  
  10.       PRINT 10, 'Smoothnum =   '
  11.       READ *, smoothnum
  12.  
  13.       DO i = 1, smoothnum
  14.         PRINT 20, 'Weight',i,':'
  15.         READ *, weights(i)
  16.       END DO
  17.  
  18.       N = 100
  19.  
  20.       PRINT 10, 'Weight divisor   '
  21.       READ *, wdivisor
  22. 10    FORMAT (1X, A15 \)
  23. 20    FORMAT (1X, A7, I3, A2\)
  24.  
  25.       DO i = 0,  n
  26.         CALL Random(rr)
  27.         xdata(i) = i
  28.         ydata(i) = COS(i/10.0) + rr
  29.       end do
  30.  
  31.       pathname = 'c:\for\lib\*.fon'
  32.       CALL InitSEGraphics(-1,pathname)
  33.       CALL SetCurrentWindow(3)
  34.       CALL ClearWindow
  35.       CALL SetAxesType(0, 0)
  36.       CALL SelectColor(2)
  37.       CALL AutoAxes(xdata, ydata, n, 0)
  38.       CALL LinePlotData(xdata, ydata, n, 1, 0)
  39.       CALL TitleWindow('original data')
  40.  
  41.       CALL SetCurrentWindow(4)
  42.       CALL SelectColor(2)
  43.       CALL SetAxesType(0, 0)
  44.       CALL DataSmoothWeights(ydata, n, smoothnum, weights,
  45.      +   wdivisor, smoothdata)
  46.       CALL AutoAxes(xdata, smoothdata, n, 0)
  47.       CALL LinePlotData(xdata, smoothdata, n, 15, 0)
  48.       CALL TitleWindow('Smoothed Data')
  49.       READ (*,*)
  50.       CALL CloseSEGraphics
  51.       END
  52.  
  53.  
  54.  
  55.