home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l292 / 1.ddi / SPLINEDE.FOR < prev    next >
Encoding:
Text File  |  1989-10-13  |  1018 b   |  32 lines
  1.       INCLUDE 'stdhdr.for'
  2.       REAL xv(0: maxr), yv(0: maxr), coef(0: maxr, 0: maxc)
  3.       INTEGER n, i
  4.       REAL r, YActual, Errr
  5.  
  6.       n = 10
  7.       DO i = 0, 10   !!one cycle of a cosine sampled at 10 intervals
  8.         xv(i) = pi * REAL(i) / 5.0       !x value
  9.         yv(i) = 3 * COS(xv(i))           !y value
  10.       END DO
  11.       !!! calculate spline equation coefficients
  12.       CALL CubicSplines(xv, yv, n, coef)
  13.  
  14.       PRINT *, 'Spline equation coefficients'
  15.       DO i = 0, n - 1
  16.         PRINT 10, i, coef(i, 0), coef(i, 1), coef(i, 2), coef(i, 3)
  17. 10      FORMAT (I2, 4(F10.2))
  18.       END DO
  19.       PRINT *, 'Press any key to continue'
  20.       READ (*,*)
  21.       PRINT *, 'Y-actual     Y-estimate   Error in estimate'
  22.       PRINT *
  23.       DO i = 0, 20
  24.         r = REAL(i) * pi / 250.0
  25.         YActual = 3.0 * COS(r)
  26.         CALL CalcSpline(xv, coef, n, r, yest)
  27.         Errr = YActual - yest
  28.         PRINT 20,  YActual, yest, Errr
  29. 20      FORMAT (3(F10.5, 5X))
  30.       END DO
  31.  
  32.       END