# 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