牛顿插值.cpp

来自「计算方法的4个小程序,牛顿插值,lagrange插值,牛顿迭代,改进欧拉公式」· C++ 代码 · 共 35 行

#include<iostream.h>
void main(void)
{
	int n,j;
	double x[100],y[100],a[100],X;
	cout<<"输入插值次数 n"<<endl;
	cin>>n;
	cout<<"输入x"<<endl;
	cin>>X;
	for(int i=0;i<n+1;i++)
	{	cout<<"输入 x["<<i<<"]"<<endl;
		cin>>x[i];
		cout<<"输入 y["<<i<<"]"<<endl;
		cin>>y[i];
		a[i]=y[i];
	}
	for(int k=0;k<n+1;k++)
	{
		for(j=n;j>k;j--)
		{
			a[j]=(a[j]-a[j-1])/(x[j]-x[j-k-1]);
		}
	}
	double sum=a[0],num=X-x[0];
	for(int m=1;m<n;m++)
	{
		num*=X-x[m];
		sum+=a[m]*num;
	}
	cout<<sum;
}