8spl10.c

第一种三次样条插值函数 SPL1为函数文件

#include "8spl1.c"
#include "stdio.h"
#define MAX 100

int  main()
  { int i;
    double s;
/*    static double dy[12],ddy[12],z[8],dz[8],ddz[8];
    static double x[12]={0.52,8.0,17.95,28.65,50.65,104.6,
                       156.6,260.7,364.4,468.0,507.0,520.0};
    static double y[12]={5.28794,13.84,20.2,24.9,31.1,36.5,
                         36.6,31.0,20.9,7.8,1.5,0.2};
    static double t[8]={4.0,14.0,30.0,60.0,130.0,230.0,
                        450.0,515.0};*/
      static double x[MAX];
      static double y[MAX];
      static double t[MAX];
      static double dy[MAX],ddy[MAX],z[MAX],dz[MAX],ddz[MAX];
      int n, m;
      printf("请输入节点数：\n");
      scanf("%d", &n);
      for(i=0; i<n; i++) {
          printf("请输入第%d个节点：", i+1);
          scanf("%lf%lf", &x[i], &y[i]);
      }    
      printf("请输入插值节点数：\n");
      scanf("%d", &m);
      for(i=0; i<m; i++) {
          printf("请输入第%d个节点的横坐标：", i+1);
          scanf("%lf", &t[i]);
      }    
      printf("请输入y1'与yn'");
      scanf("%lf%lf", &dy[0], &dy[n-1]);
    printf("\n");
    s=spl1(n, m, x,y,dy,ddy,t,z,dz,ddz);
    printf(
"      x(i)          y(i)          dy(i)         ddy(i)\n");
    for (i=0;i<n;i++)
   printf("%14.6e%14.6e%14.6e%14.6e\n",x[i],y[i],dy[i],ddy[i]);
    printf("\n");
    printf("s=%e\n",s);
    printf("\n");
    printf(
"      t(i)          z(i)          dz(i)         ddz(i)\n");
    for (i=0;i<m;i++)
   printf("%14.6e%14.6e%14.6e%14.6e\n",t[i],z[i],dz[i],ddz[i]);
   printf("\n");
   while(1);
  }