Duffing = {x'[t] == y[t], y'[t] == -0.25 y[t] + x[t] - x[t]^3 + 0.3 Cos[t]}
//定義 Duffing 這個微分方程組系統。並設定一些基本參數。
sol = NDSolve[{Duffing, x[0] == 1, y[0] == 0.1}, {x, y}, {t, 0, 100},
MaxSteps -> 10000]
//利用 numerical solve 微分方程組，並設定最大運作數
ParametricPlot[Evaluate[{x[t], y[t]} /. sol], {t, 0, 100}]
//利用參數式 plot 相圖。
//然後在設定其他參數就可以觀察是否出現非週期解的現象了!!
--
"因為,我們只是生命臻向完美的途中所出現的物種而已..."
==== 張草,明日滅亡 ====
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 220.139.209.50