牛顿前插值.c

来自「计算方法的牛顿前插值算法vc实现」· C语言 代码 · 共 860 行 · 第 1/3 页

#include"stdio.h"
#include"math.h"

main()
{
 int n,i,j;
 float x[10],y[10],z,h,t,s;
 printf("请输入自变量X个数:\n");
   scanf("%d",&n);
 printf("请输入自变量X值:\n");
   for(i=0;i<n;i++)
     scanf("%f",&x[i]);
 printf("请输入因变量Y值:\n");
   for(i=0;i<n;i++)
     scanf("%f",&y[i]);
 printf("请输入要求的自变量X的值:\n");
   scanf("%f",&z);
   h=x[1]-x[0];
   t=(z-x[0])/h;
   for(i=1;i<n;i++)
   {
	   for(j=n-1;j>=i;j--)
	   {
		   y[j]=(y[j]-y[j-1])*h/(x[j]-x[j-i]);
	   }
	}
   for(i=1;i<n;i++)
   {
	   for(j=i;j<n;j++)
	   {
		   y[j]=y[j]*t;
	   }
		   t--;
   }
 		   s=0;
		   for(i=0;i<n;i++)
		   {
			   s=s+y[i];
		   }
		   printf("所求自变量X的函数值为:");
		   printf("%f",s);
}

/*