计算方法_拉格朗日插值_C++实现_方法1

#include <iostream> 
#include<cmath>
using namespace std;
void La(float (* xy)[8],float x,int n)
{
   float fx(0),t;
   int i,j;
   for(i=0;i<=n;i++)                        
   {
       t=1;
       for(j=0;j<=n;j++)
       {
           if(j==i)
               continue;
           else
               t=t*(x-xy[0][j])/(xy[0][i]-xy[0][j]);//求插值基函数
       }
       fx=fx+t*xy[1][i];
   }
   cout<<"插值点"<<x<<"的值为："<<fx<<endl;
}


int main()
{
    cout<<"次程序最多可以用于计算8个节点的插值"<<endl;
    float xy[8][8];
    int i,j;
    float m;
    cout<<"请输入节点个数：";
    cin>>j;

    cout<<"请 输 入 x的 值：";   //输入各个节点的x 的值
    for(i=0;i<j;i++)
        cin>>xy[0][i];
    cout<<"请输入对应的y值：";    //输入对应于x 的y的值
    for(i=0;i<j;i++)
        cin>>xy[1][i];
    cout<<"请输入需要求的插值点m：";     //输入对分的次数
    cin>>m;
    cout<<"用拉格朗日插值法求解为："<<endl;
    La(xy,m,j);                              //用拉格朗日法求解
   
    return 0;
}