について一度微分した最初の方程式をt
使用して を代用できることに注意する必要がありv[t]
ます。しかし、2 番目の方程式は 2 次の ODE になり、別の追加の初期条件を指定する必要があります。あげます
v[0]==x'[0]==some number
次に、この ODE を解いてx
回復できるようにした後、幾何学的に状況が明確になるようv[t]==x'[t]
に、a の観点から解を示します。Manipulate
(* First equation *)
v[t] = x'[t];
(*
Differentiate this equation once and substitute
for v[t] in the second equation
*)
Manipulate[
With[{u = Constant, der = derval},
res = NDSolve[{x''[t] == -x[t] - u*x'[t]^3, x[0.] == 2,x'[0.] == der},
x, {t, 0., 30.}] // First;
Plot[Evaluate[{x[t], v[t]} /. res], {t, 0, 30}, PlotRange -> All,
Frame -> True,Axes -> None, ImageSize -> 600]
],
{{Constant, 0.,TraditionalForm@(u)}, 0.,3, .1},
{{derval, -3., TraditionalForm@(v[0] == x'[0])}, -3, 3, .1}
]

これがお役に立てば幸いですが、次回質問する前に、まず理論をブラッシュアップする必要があります。あなたが質問した質問は、Mathematica プログラミングではなく、非常に基本的で初歩的な数学に関するものであることがわかります。幸運を!!