拉格朗日插值.cpp

来自「数值分析程序。。。。。好多程序下载 喜欢的朋友可以下 随时欢迎」· C++ 代码 · 共 24 行

#include <iostream.h>
void main()
{
	double x1,y1,t;
    int n,i,j;
	cout<<"请输入n,x1的值："<<endl;
	cin>>n>>x1;
    double x[]={1,2,3,4,5,6,7};
    double y[]={1,1.414214,1.732051,2,2.236068,2.449490,2.645751};
    y1=0.0;
	for(i=0;i<=n;i++)
	{ 
		t=1.0;
	   for(j=0;j<=n;j++)
	   {
	   if(j!=i)
	    t=t*(x1-x[j])/(x[i]-x[j]);
	   }
       y1=y1+t*y[i];
	}
      cout<<y1<<endl;
 
}