Lagrange の補間公式を用いて関数値を計算する手順
計算機で,いろいろなaに対して, (2) を計算するには, "チm; (i=0,1, , n) のとき
Pn(z)
(- 2o)( - 1) (日 1 ) ×
1 f(r)
三
2- C; di
4 = (T- o).… (g- m-1)(H- 41)… (エ-En)
da
となることから16, 次の手順でするのがよいでしょう.このようにするのが速いとされています. ただし, a
におけるf(z) の推定値 faを計算するものとします。
1. Co, I1,…
Tn をそれぞれ配列の要素 a[0], z[1], , a[n] に代入する。
2. f(zo),f(z1),…, f(エn) をそれぞれ配列の要素 f[O], f[1]),……, f\n]に代入する
3. a を与える。
4.i= 0,1,2, ..,n に対して, 順に,
di] = (z] - [0)) * * (z[i]-2i-1]) * (z]-2i+1)) * ** (z]-alnl)
を計算する。
5.c=1, s=0 と初期化し,