ode_discrete 常微分方程式ソルバ, 離散時間シミュレーション

## 呼び出し手順

    y = ode("discrete",y0,k0,kvect,f)

## 引数

y0

  実数のベクトルまたは行列 (初期条件).

t0

  実数スカラー (初期時間).

f

  外部 (関数または文字列またはリスト).

k0

  整数 (初期時間).

kvect

  整数ベクトル.

## 説明

(最初の引数を`"discrete"`とする ) この構文により,`ode`は 初期状態量 `y(k0)` から再帰的に `y(k+1)=f(k,y(k))`を計算し, `kvect`の中にある`k`について `y(k)`を返します. `kvect(1)` は`k0`以上である 必要があります.

他の引数およびその他のオプションは `ode`と同じです. odeのヘルプを参照ください.

## 例

    y1 = [1;2;3];
    deff("yp = a_function(k,y)", "yp = A*y + B*u(k)")
    A = diag([0.2,0.5,0.9]);  B = [1;1;1];  u = 1:10;  n = 5;
    y = ode("discrete", y1, 1, 1:n, a_function);
    y(:,2) - (A*y1+B*u(1))

    // ここで, [y3,y5,y7,y9] で y を評価します
    y = ode("discrete", y1, 1, 3:2:9, a_function)

## 参照

- ode