动态调整微分方程的参数
使用一个 Manipulate 的相位平面图改变一个有 NDSolve 解出的微分方程的初始条件.
Manipulate[
Module[{sol = NDSolve[{y''[x] + Sin[y[x]] == 0, y[0] == p[[1]], y'[0] == p[[2]]}, y, {x, 0, T}]},
ParametricPlot[Evaluate[{y[x], y'[x]} /. sol], {x, 0, T}, PlotRange -> {{-4, 4}, {-3, 3}}]],
{{p, {2, 1}}, Locator}, {{T, 10}, 0, 100}]