拉格朗日.cpp

来自「牛顿 拉格朗日 分段线性」· C++ 代码 · 共 30 行

#include<iostream.h>
void main()
{
	const int N=3;
    double a[N],b[N];
	cout<<"Input x[N]:";
	for(int i=0;i<N;i++)
		cin>>a[i];
    cout<<"Input y[N]:";
	for(i=0;i<N;i++)
		cin>>b[i];
	double x,y=0,t;
	cout<<"x=";
	cin>>x;

	for(int k=0;k<N;k++)
	{
		t=1;
		for(int j=0;j<N;j++)
		{
			if(j==k) continue;
			t*=(x-a[j])/(a[k]-a[j]);
		}
		y+=t*b[k];
	}
	cout<<"y="<<y<<endl;
}