// --------------------------------------------------------------------------
//  関数RungeKuttaを用いて例題3 (p.141)を解くスクリプト
// --------------------------------------------------------------------------
clear;                                                   // 定義済みの変数の値の消去
exec('RungeKutta.sci');                                  // 関数の読み込み
deff('dy=f(x,y)','dy=y-y^2');                            // 例題3の微分方程式の関数の定義
a = 0;                                                   // 区間の始まり
b = 20;                                                  // 区間の終わり
y0 = 0.1;                                                // 初期条件
N =  8; [x,y] = RungeKutta(f,a,b,y0,N); plot2d(x,y,1);   // N = 8  の場合
N = 10; [x,y] = RungeKutta(f,a,b,y0,N); plot2d(x,y,2);   // N = 10 の場合
N = 20; [x,y] = RungeKutta(f,a,b,y0,N); plot2d(x,y,3);   // N = 20 の場合