计算方法_牛顿插值_C++实现

#include<iostream>

#include<ctime>

using namespace std;

 

float * f(int n,float *x,float *y)

{

    //该函数的返回值是一个长度为n＋1的动态数组，存放各阶差商

    float * f = new float[n+1];

        

    f[0] = y[0];

    //cout<<f[0]<<endl;

    for (int i = 1; i <= n; i++) {

        float t = 1;

        float sum = 0;

        for (int k = 0; k <= i; k++) {

            for (int j = 0; j <= i; j++) {

                if (j != k) {

                    t=t*(x[k] - x[j]);

                }            

            }

            sum = sum + y[k] / t;

            t=1;    

        }    

        f[i] = sum;

    

    }

        

    return f;

}

int main(){

    

    float x[5]={0,2,3,5,6},y[5]={0,8,27,125,216};

    int n=4;

    

    float * answer=f(4,x,y);

    for (int i = 0; i <= n; i++) {

        cout<<answer[i]<<endl;

    }

    

    //2阶差乘

    cout<<y[0]/((x[0]-x[1])*(x[0]-x[2]))+y[1]/((x[1]-x[0])*(x[1]-x[2]))+y[2]/((x[2]-x[0])*(x[2]-x[1]))<<endl;

    

        float t = 1;

        float sum = 0;

        for (int k = 0; k <= 2; k++) {

            for (int j = 0; j <= 2; j++) {

                if (j != k) {

                    t=t*(x[k] - x[j]);

                }

                        

            }

            

            sum = sum + (y[k] / t);    

            t=1;

        }    

        

    cout<<sum;

    

    

    

    

    return 0;

}