newtonchazhi.c

来自「两种曲线拟合的算法」· C语言 代码 · 共 36 行

#include <stdio.h>
main()
{float s,p,x,X[10],Y[10];
 int i,j,n;
 
 printf("请输入结点个数(0<n<11)：");
 scanf("%d",&n);
 printf("\n");
 printf("请输入结点：\n");                  
 for(i=0;i<n;i++)
 {
	 scanf("%f,%f",&X[i],&Y[i]);
   	 
 }
 printf("\n");
 printf("请输入x值：");
 scanf("%f",&x);
 printf("\n");
 
 for(i=1;i<n;i++)
	 for(j=n-1;j>i;j--)
		 Y[j]=(Y[j]-Y[j-1])/(X[j]-X[j-i]);
 p=1;s=Y[0];
 for(i=1;i<n;i++)
 {
	 p=p*(x-X[i-1]);
	 s=s+Y[i]*p;
 }
 
 printf("此题的结果为：");
 printf("N%d(%f)=%f",n-1,x,s);
 printf("\n");
}