Original Mathematica File

第一章 微分方程式

§1 微分方程式と曲線群

問1

"DFChap1ver6_1.gif"の曲線群の微分方程式を求めよう。

In[3]:=

"DFChap1ver6_2.gif"

Out[3]=

"DFChap1ver6_3.gif"

実際、パラメターcを変化させて曲線群をプロットしてみる。

In[63]:=

"DFChap1ver6_4.gif"

Out[63]=

"DFChap1ver6_5.gif"

Manipulateを使ってインタラクティブに操作してみよう。

In[68]:=

"DFChap1ver6_6.gif"

Out[68]=

"DFChap1ver6_7.gif"

次に "DFChap1ver6_8.gif" の曲線群を求めてみよう。

In[2]:=

"DFChap1ver6_9.gif"

In[3]:=

"DFChap1ver6_10.gif"

Out[3]=

"DFChap1ver6_11.gif"

アウトプットをTraditional form にすると

Out[3]=

"DFChap1ver6_12.gif"

In[11]:=

"DFChap1ver6_13.gif"

Out[11]=

"DFChap1ver6_14.gif"

In[16]:=

"DFChap1ver6_15.gif"

Out[16]=

"DFChap1ver6_16.gif"

問2

 "DFChap1ver6_17.gif" の曲線群を求めてみよう。

In[17]:=

"DFChap1ver6_18.gif"

In[18]:=

"DFChap1ver6_19.gif"

Out[18]=

"DFChap1ver6_20.gif"

または、

Out[18]=

"DFChap1ver6_21.gif"

次に y-(b cos(2 x)+a sin(2 x))==0 の曲線群を求めてみよう。

In[19]:=

"DFChap1ver6_22.gif"

Out[19]=

"DFChap1ver6_23.gif"

または、

Out[19]=

"DFChap1ver6_24.gif"

問題

1.(1)

 "DFChap1ver6_25.gif" の曲線群を求める。

In[20]:=

"DFChap1ver6_26.gif"

Out[20]=

"DFChap1ver6_27.gif"

または、

Out[20]=

"DFChap1ver6_28.gif"

In[25]:=

"DFChap1ver6_29.gif"

Out[25]=

"DFChap1ver6_30.gif"

In[31]:=

"DFChap1ver6_31.gif"

Out[31]=

"DFChap1ver6_32.gif"

1.(2)

"DFChap1ver6_33.gif" の曲線群を求める。

In[32]:=

"DFChap1ver6_34.gif"

Out[32]=

"DFChap1ver6_35.gif"

または、

Out[32]=

"DFChap1ver6_36.gif"

Mathematicaで逆に微分方程式を解いてみよう。以下の解が元の
式と等しいことに注意せよ。("DFChap1ver6_37.gif"がcとなる。)

In[33]:=

"DFChap1ver6_38.gif"

Out[33]=

"DFChap1ver6_39.gif"

In[38]:=

"DFChap1ver6_40.gif"

Out[38]=

"DFChap1ver6_41.gif"

1.(3)

"DFChap1ver6_42.gif" の曲線群を求める。

In[3]:=

"DFChap1ver6_43.gif"

Out[3]=

"DFChap1ver6_44.gif"

Out[3]=

"DFChap1ver6_45.gif"

Mathematicaで逆に微分方程式を解いてみよう。以下の解が元の
式と等しいことに注意せよ。

In[4]:=

"DFChap1ver6_46.gif"

Out[4]=

"DFChap1ver6_47.gif"

In[5]:=

"DFChap1ver6_48.gif"

Out[5]=

"DFChap1ver6_49.gif"

2,(1)

x 軸上に中心を持つ半径2の円の群の微分方程式は

In[6]:=

"DFChap1ver6_50.gif"

Out[6]=

"DFChap1ver6_51.gif"

または、

Out[6]=

"DFChap1ver6_52.gif"

In[15]:=

"DFChap1ver6_53.gif"

Out[15]=

"DFChap1ver6_54.gif"

2,(2)

原点を焦点とし、軸がx軸の放物線の群 "DFChap1ver6_55.gif"

In[16]:=

"DFChap1ver6_56.gif"

Out[16]=

"DFChap1ver6_57.gif"

または、

Out[16]=

"DFChap1ver6_58.gif"

In[19]:=

"DFChap1ver6_59.gif"

Out[19]=

"DFChap1ver6_60.gif"

In[24]:=

"DFChap1ver6_61.gif"

Out[24]=

"DFChap1ver6_62.gif"

In[25]:=

"DFChap1ver6_63.gif"

Out[25]=

"DFChap1ver6_64.gif"

2,(3)

原点が中心、2直線 y=+-2xを漸近線とする双曲線の群

In[20]:=

"DFChap1ver6_65.gif"

Out[20]=

"DFChap1ver6_66.gif"

または、

Out[20]=

"DFChap1ver6_67.gif"

In[32]:=

"DFChap1ver6_68.gif"

Out[32]=

"DFChap1ver6_69.gif"

§2 微分方程式の解

例題1

同心円群の微分方程式を解いてみよう。

In[33]:=

"DFChap1ver6_70.gif"

Out[33]=

"DFChap1ver6_71.gif"

これは"DFChap1ver6_72.gif" と同等であるから、曲線群を描くと

In[35]:=

"DFChap1ver6_73.gif"

Out[35]=

"DFChap1ver6_74.gif"

例題2

"DFChap1ver6_75.gif" を解いてみよう。

In[31]:=

"DFChap1ver6_76.gif"

Out[31]=

"DFChap1ver6_77.gif"

これは y = (x - c)^2 と同等であるので、その曲線群は

In[36]:=

"DFChap1ver6_78.gif"

Out[36]=

"DFChap1ver6_79.gif"

演習問題 1−1 [A]

1 (1)

In[38]:=

"DFChap1ver6_80.gif"

In[39]:=

"DFChap1ver6_81.gif"

"DFChap1ver6_82.gif"

Out[39]=

"DFChap1ver6_83.gif"

1 (2)

In[40]:=

"DFChap1ver6_84.gif"

"DFChap1ver6_85.gif"

Out[40]=

"DFChap1ver6_86.gif"

1 (3)

In[41]:=

"DFChap1ver6_87.gif"

Out[41]=

"DFChap1ver6_88.gif"

または、

Out[41]=

"DFChap1ver6_89.gif"

1 (4)

In[42]:=

"DFChap1ver6_90.gif"

Out[42]=

"DFChap1ver6_91.gif"

Out[42]=

"DFChap1ver6_92.gif"

2 (1)

まず、"DFChap1ver6_93.gif"をプロットする。

In[43]:=

"DFChap1ver6_94.gif"

Out[43]=

"DFChap1ver6_95.gif"

接線は x = c のとき 傾きがcで yは c^2/2 の値を取るので、曲線群は "DFChap1ver6_96.gif".

In[50]:=

"DFChap1ver6_97.gif"

Out[50]=

"DFChap1ver6_98.gif"

"DFChap1ver6_99.gif"

曲線群の微分方程式は以下のように求めることができる。

In[41]:=

"DFChap1ver6_100.gif"

In[42]:=

"DFChap1ver6_101.gif"

Out[42]=

"DFChap1ver6_102.gif"

Out[42]=

"DFChap1ver6_103.gif"

Spikey Created with Wolfram Mathematica 6