N次元連立微分方程式のMATLABでの解法コード

Question

正典小林 2020-9-20

回答： michio 2020-9-24

n次元連立微分方程式のMATLABでの解法コードを教えてください

具体的な例としては　x、y、z　の3方向のばねで支えられた球の3方向から変動力を受けた時の3次元運動（変位と速度）

dy1/dt=f1(t,y1,y2,・・・,yn)

dy2/dt=f2(t,y1,y2,・・・,yn)

・

dyn/dt=fn(t,y1,y2,・・・,yn)

Answer 1

michio 2020-9-24

例えば ode45 関数のヘルプページ（https://jp.mathworks.com/help/matlab/ref/ode45.html）の

"ノンスティッフ方程式の求解" や "時間依存の項をもつ ODE" のサンプルコードは 2 次元連立微分方程式を解いています。まずはこちらを参考にして頂くとよいかと思います。