« Scilabで微分方程式 | トップページ | 「CPUの創りかた」買いました。 »

Scilabで2階の微分方程式

さらに頑張って2階の微分方程式までやってみました。

function xdot=myfunc(t,x)
xdot=zeros(2,1);
xdot(1)=x(2);
xdot(2)=-4.91*sin(x(1));
endfunction

tt=0:0.1:10;
xx=ode([%pi/2 ; 0], 0, tt, myfunc);
xbasc();
plot2d(tt,xx(1,:),1);
plot2d(tt,xx(2,:),2);

とすると図のようになります。連立方程式を配列を使って表現するんですね。
MATLAB つかいませんかを参考にしています。

|

« Scilabで微分方程式 | トップページ | 「CPUの創りかた」買いました。 »

コメント

コメントを書く



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


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



トラックバック

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

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

« Scilabで微分方程式 | トップページ | 「CPUの創りかた」買いました。 »