home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c129 / 1.ddi / SPLINEDE.C < prev    next >
Encoding:
C/C++ Source or Header  |  1988-02-02  |  1019 b   |  47 lines
  1. # include <math.h>
  2. # include "miscio.h"
  3. # include "curvefit.h"
  4.  
  5.  
  6. float xv[11];
  7. float yv[11];
  8. float coef[13][4];
  9. int i;
  10. int j;
  11. int n;
  12. float x;
  13. float YActual;
  14. float YEst;
  15. float r;
  16. float error;
  17. char c;
  18.  
  19. void main()
  20. {
  21.    n = 10;
  22.    for ( i = 0; i <= 10; ++i ) {
  23.       xv[i] = M_PI * i / 5.0;
  24.       yv[i] = 3.0 * cos(xv[i]);
  25.    }
  26.    CubicSplines(xv,yv,n,&coef[0][0]);
  27. /*   ClrScr();
  28.   */  printf( "Spline equation coefficients\n");
  29.    for ( i = 0; i <= n - 1; ++i ) {
  30.       printf("%6d %8.3f %8.3f %8.3f %8.3f\n", i, coef[i][0],
  31.              coef[i][1], coef[i][2], coef[i][3]);
  32.    }
  33.    printf("Press return to continue");
  34.    c = getch();
  35.    printf("\n");
  36.    printf("\n");
  37.    printf("Y-actual Y-estimate Error in estimate\n");
  38.    printf("\n");
  39.    for ( i = 0; i <= 20; ++i ) {
  40.       r = i * M_PI / 250.0;
  41.       YActual = 3.0 * cos(r);
  42.       CalcSpline(xv, &coef[0][0],n,r,&YEst);
  43.       error = YActual - YEst;
  44.       printf("%8.3f %8.3f %8.3f\n", YActual,  YEst, error);
  45.    }
  46. }
  47.