lagelangri.cpp

拉格朗日算法插值源程序--对于学习数值计算的是一个很好资料

//拉格朗日算法插值源程序:
#include<iostream.h>
#define c 10
int main()
{
	int n,a=1,b=1;
	double j,k=1,result=0,x[c],y[c];

	cout<<"输入插值结点数目:   ";
	cin>>n;
	cout<<"输入结点的坐标:   x y"<<endl;
	while(a<=n)
	{
		cout<<"第 "<<a<<" 个点的坐标为  ";
		cin>>x[a]>>y[a];
		a++;
	}
	cout<<"输入所求插值结点x坐标:   ";
	cin>>j;

	for(a=1;a<=n;a++)
	{
		for(b=1,k=1;b<=n;b++)
			if(b!=a)k=k*(j-x[b])/(x[a]-x[b]);
		result=result+k*y[a];
	}
	cout<<"所求y坐标为:  "<<result<<endl;
}