「微分方程式をMathematica で解く」
 
最初にお断りしておくが、以下の式・図は
 
「Mathematica 微分方程式(小林 道正著)」
 
を参考にさせていただいた。アニメーションはオリジナルである。
 [Graphics:math04gr2.gif][Graphics:math04gr1.gif]
 
[Graphics:math04gr3.gif]
 
 
[Graphics:math04gr2.gif][Graphics:math04gr4.gif]
[Graphics:math04gr2.gif][Graphics:math04gr5.gif]
[Graphics:math04gr2.gif][Graphics:math04gr6.gif]
[Graphics:math04gr2.gif][Graphics:math04gr7.gif]
[Graphics:math04gr2.gif][Graphics:math04gr8.gif]
 
 
[Graphics:math04gr2.gif][Graphics:math04gr9.gif]
 
 
[Graphics:math04gr2.gif][Graphics:math04gr10.gif]
[Graphics:math04gr2.gif]
 

  [Graphics:bibungr2.gif][Graphics:bibungr1.gif] [Graphics:bibungr2.gif][Graphics:bibungr3.gif] [Graphics:bibungr2.gif][Graphics:bibungr4.gif]

[Graphics:bibungr2.gif](初期値の違いによるグラフの例)

初期値c を、1から12まで変化させると次のようになる。

 

「変数分離形の微分方程式」の例

 y'=ax(b-y) の形 を解くと

[Graphics:bibun2gr2.gif][Graphics:bibun2gr1.gif]
[Graphics:bibun2gr2.gif][Graphics:bibun2gr3.gif]
 a=b=1,初期値をy(0)=k としたとき、k を変化させたときのグラフは次のようになる。
 
 
[Graphics:bibun2gr2.gif][Graphics:bibun2gr4.gif]
[Graphics:bibun2gr2.gif][Graphics:bibun2gr5.gif]

[Graphics:bibun2gr2.gif][Graphics:bibun2gr6.gif]

[Graphics:bibun2gr2.gif][Graphics:bibun2gr7.gif]
次に初期値を y(0)=2、と固定して b の値を変えたときの解の変化を式とアニメーションで表してみると、
 
[Graphics:bibun2gr2.gif][Graphics:bibun2gr8.gif]
[Graphics:bibun2gr2.gif][Graphics:bibun2gr9.gif]

[Graphics:bibun2gr2.gif]

 

    Mathematicaについて   高校数学とMathematica   数学あれこれ   数学のページへ トップページへ