main()
{int i,n;
int k;
int x[5]={0,1,2,3,4};
float y[5]={1,0.95,0.82,0.71,0.43};
float p,s;
float b;
s=0.0;
printf("input b,n:");
scanf("%f,%d",&b,&n);

for(k=0;k<n;k++)
{p=1.0;
for(i=0;i<n;i++)
if(k!=i)
p*=(b-x[i])/(x[k]-x[i]);
s+=p*y[k];
}

printf("%5.4f",s);
}
看下有错么

运行精确不够

N阶拉格朗日
是什么样子的呀？

楼上的不是要我的命吗

您的程序是没错的,我昨天去看了下那本书.
您所说的运行精确不够,那是因为您的测试数据过少，自然精确不够。
您可以多加几组数据，把n直接赋成整个数组成员总数。这样就可以精确些.
不过,对于您的这个拉格朗日的插值在现实中是几乎不用的.
运用较广的是牛顿均差插值法.您想了解的话可以去查相关资料看下。

[此贴子已经被作者于2006-5-17 16:35:55编辑过]