牛顿插值多项式算法

NewtonDifference.h

#pragma once

class NewtonDifference
{
public:
	NewtonDifference(void);
	virtual ~NewtonDifference(void);

	// 三阶插值函数
	float f3(float x);

	// 牛顿插值多项式算法
	void NewInt(double* x, double* y
		, int n, double& xi, double* yint, double* ea);
};

NewtonDifference.cpp

#include "NewtonDifference.h"
#include <string.h>
#include <stdio.h>

NewtonDifference::NewtonDifference(void)
{
    //f3(2);

    // 开始牛顿插值
    double *x = new double[4];
    x[0] = 1.0f;
    x[1] = 4.0f;
    x[2] = 6.0f;
    x[3] = 5.0f;

    double *y = new double[4];
    y[0] = 0.0f;
    y[1] = 1.386294f;
    y[2] = 1.791759f;
    y[3] = 1.609438f;

    int n = 3;
    double * yint = new double[n];
    double * ea =