« Maximaで微分方程式 | トップページ | Scilabで2階の微分方程式 »

Scilabで微分方程式

微分方程式のはなし9章ではコンピュータを使って微分方程式を解いてます。というわけで私もやってみました。ただしルンゲ・クッタ法を自分でインプリメントしたわけでなく、Scilabを使って解いてみただけです。

deff("wdot=df(t, w)" , "wdot=-0.139*w")
t0=0 ; w0=1000 ; t=0:1:30;
y=ode(w0 , t0 , t , df );
plot(t,y)

とすると下のようなグラフが得られます。かんたん、かんたん(^^)

|

« Maximaで微分方程式 | トップページ | Scilabで2階の微分方程式 »

コメント

コメントを書く



(ウェブ上には掲載しません)


コメントは記事投稿者が公開するまで表示されません。



トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/44043/11180576

この記事へのトラックバック一覧です: Scilabで微分方程式:

« Maximaで微分方程式 | トップページ | Scilabで2階の微分方程式 »