1. 算法:
Newton插值_子函数形式。
2. 算法组织:
参见《计算方法教程(第2版)》,凌永祥、陈明奎著,西安交通大学出版社,2005年4月第二版,第91~92页。
3. 程序代码(执行环境为VS2005):
//Newton插值_子函数形式
#include "stdafx.h"
#include
#include //exp(const complex& _ComplexNum)函数
#include //使用三角函数及log函数
void DD(double x[], double y[], const int n);//牛顿插值算法。x[],y[]存储插值数据点(xi, yi),算法DD计算牛顿插值多项式的系数
//(即各阶差商),并将系数存在y[]中。n为插值数据点数组的最大下标值,即有n+1个插值数据点。
void NI(double x[], double y[], const int n, double x_point, double& y_point);//对于给定点x_point,本算法用Newton插值公式计算值
//N(x)并存放到y_point中。
int _tmain(int argc, _TCHAR* argv[])
{
using namespace std;
//using namespace Math;
const int N = 3;//数组下标必须为确定值,不能为变量
double px = 2.0;//px为插值点