lagrange.cpp

来自「这是拉格朗日插值的源码」· C++ 代码 · 共 52 行

/*该程序完成拉格朗日插值法通用子程序*/
                         #include"stdafx.h"
                        #include <iostream.h>
class Lagrange
{
public:
	void flag(void);
private:
	float *p,*q; //用指针动态分配内存空间
	int n;
};
void Lagrange::flag(void)
{    
	int i,j;
	float x,y,t;
	  cout<<"请输入插值阶数:N=";
    	cin>>n;
    	cout<<"请一共输入"<<n+1<<"个插值结点:"<<endl;
	n++;
	p=new float[n];  //据输入n值指针定位
	q=new float[n];
	y=0;
	for (i=0;i<n;i++)
	{   
		cout<<"第"<<i+1<<"个结点 ";
		cin>>p[i]>>*(q+i); //指针的两种表示形式
	}
	  cout<<"请输入X=";
	  cin>>x;
            for (i=0;i<n;i++)    /*据拉格朗日公式*/
			{
	    	t=1;
		       for (j=0;j<n;j++)
			   {
			    if (i!=j)      
		        t*=((x-p[j])/(p[i]-p[j]));
			   }
		      y+=t*q[i];
			}
	cout<<"y="<<y<<endl;
}
void main()
{
    Lagrange a;
	a.flag();
}