雖然向前差除和向後差除在計算順序上似乎有些不同，但對程式設計而言，事實上前差法和後差法是一樣的東西，只是表格輸出不同而已。說明圖如下所示。
/*
Newton Forward and Backward Interpolation Algorithm
Assume function is -4x^4 + 2x^3-3x^2+4x-1
- 11月 28 週日 201006:01
[C] 牛頓 向前差除/向後差除 內差法
雖然向前差除和向後差除在計算順序上似乎有些不同，但對程式設計而言，事實上前差法和後差法是一樣的東西，只是表格輸出不同而已。說明圖如下所示。
/*
Newton Forward and Backward Interpolation Algorithm
Assume function is -4x^4 + 2x^3-3x^2+4x-1
- 11月 28 週日 201003:26
[C] Lagrange 內插法
/*
Largrange Interpolation Algorithm
(1) give n, give x[0..n], give y[0..n]
(2) Ln,k(x) = (x-x0)*(x-x1)*...*(x-xn) / (xk-x0)*(xk-x1)-...*(xk-xn)
(3) P(x) = Ln1(x)f(x1)+Ln2(x)f(x2)+...+Lnn(x)f(xn)
assume function = -4x^4 + 2x^3-3x^2+4x-1
-by EdisonX
*/
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>
#define N 4
Largrange Interpolation Algorithm
(1) give n, give x[0..n], give y[0..n]
(2) Ln,k(x) = (x-x0)*(x-x1)*...*(x-xn) / (xk-x0)*(xk-x1)-...*(xk-xn)
(3) P(x) = Ln1(x)f(x1)+Ln2(x)f(x2)+...+Lnn(x)f(xn)
assume function = -4x^4 + 2x^3-3x^2+4x-1
-by EdisonX
*/
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>
#define N 4
1
開發手札 (2)