Mathematica在微分方程上的应用.PDF

Mathematica 在微分方程上的應用  驗證微分方程式的函數解 利用Mathematica 可以輕鬆地驗證任一函數是否為某微分方程式的解。例如 若要驗證y (x) x ex 是方程式y 2y y 0 的解必須先定義y (x) 函數。在 Mathematica 中定義一個函數的寫法跟傳統數學函數的寫法相當相似，不過需將 小括弧改成中刮號，刮號中的自變數名稱的後面也需加上一橫底線： y x_  x x x  x 定義完所要驗證的函數後，就可以馬上套用，執行下面的指令： y x 2 y x y x Simplify 0 上式代表y (x) x ex 果然符合微分方程y 2y y 0 ，注意到上述的指令不能 [x] 漏掉 項。 磱껳 Mathematica 的內建指令Plot y (x) ，下圖繪出 函數及其一階與二階導 Hue (Hue 0) (Hue 數的函數曲線，並利用不同的 值設定分別加以紅色 值為 、綠色 值為0.3) 與藍色 (Hue 值為0.6) 著色區分之。 Plot y x , y x , y x , x, 1, 1, PlotStyle  Hue 0, Hue.3, Hue.6 . 1 . 4 2 -1 -0.5 0.5 1 -2 -4 Graphics  [1, 1] 3 4 y (x) 最下面一條曲線在 定義域範圍內皆位於第 及第 向限中，隱含 函數 在該區間內為一嚴格凹函數 (strickly concave function)。 除了顯函數之外，Mathematica 也可以驗證微分方程式的隱函數解。下面我 2 dy 2 y 4 x 5 們將以隱函數2 x y 2 x y 5 x 0 為例，驗證其為 的解。首 dx 2 (y x) 先仍先需定義出所需驗證的隱函數： 2 2 eq 2 x y 2 x y 5 x 0 2 2 5 x 2 x 2 x y y 0 eq == == 上述指令將隱函數命名為 ，並利用符號 連結該式之左式與右式。符號 一 方面可以用以判斷左式與右式是否相等，例如： 1  0 False 所謂的相等是指抽象概念上相等，而不是說左式與右式一定要一模一樣，例如下 == 面兩個例子分別顯示 可用於跨越整數與實數間、甚至符號算式間是否相等的判 斷：