// --------------------------------------------------------------------------------
//  テキストp.24 例題2
//  4つの補間点から3次のラグランジュ補間多項式を求め，3√20の推定値を求めるスクリプト
// --------------------------------------------------------------------------------
clear;                                    // 定義済みの変数の値を全て消去
exec('Lagrange.sci');                     // Lagrange関数の読込み

x = [19.9466, 19.9907, 20.0349, 20.0793]; // 補間点のx軸の値の系列
y = [  2.712,   2.714,   2.716,   2.718]; // 補間点のy軸の値の系列
f = Lagrange(x, y)                        // Lagrange関数でラグランジュ補間多項式fを求める
yy = horner(f, 20);                       // 補間多項式fに20を代入
printf("f(20) = %f\n", yy);               // 3√20の推定値f(20)の値の表示