牛顿插值公式求结点的值C语言实现

数值计算基础的实验之一 插值方法 牛顿插值公式求结点的值C语言实现 代码摘要: double Newton(double x) { int i, j; // 求均差 for (i =0; i<n; i++) A[i][0] = Y[i]; for (i=1; i<n; i++) for (j =1; j<=i; j++) A[i][j] = (A[i][j-1] - A[i-1][j-1]) / (X[i] - X[i-j]); // 求结点 double result = A[0][0]; for (i=1; i<n; i++) { double tmp = 1.0; for (int j=0; j<i; j++) tmp *= (x - X[j]); result += tmp * A[i][i]; } return result; }