lag.c

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

#include <stdio.h>
main()
{float s,p,x,X[10],Y[10];
 int i,j,n;
 s=0;
 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=0;i<n;i++)
 {   p=1;
	 for(j=0;j<i;j++)
		 p=p*(x-X[j])/(X[i]-X[j]);
	 for(j=i+1;j<n;j++)
		 p=p*(x-X[j])/(X[i]-X[j]);
     s=s+p*Y[i];
 }
 printf("此题的结果为：");
 printf("f(%f)=%f",x,s);
 printf("\n");
}