大佬们，想用mathematica解一个二阶微分方程（微振动），画出x（t）和x'（t）的图像；但是t如果太大，图像会“飞”出去，所以想用for循环，利用能量守恒为判据，来确定最大的t，但不知道为什么会报错...第一次用mathematica...
x[0] = 0
x'[0] = 0.4
C = 0.5 x'[0]^2 + 1/3 x[0]^3 + 0.5 x[0]^2
For[i = 0, 0.5 x'[i]^2 + 1/3 x[i]^3 + 0.5 x[i]^2 <= C, i + 0.01,
NDSolve[{x''[t] + x[t] + x[t]^2 == 0, x[0] == 0, x'[0] == 0.6},
x, {t, 0, i}]]
ParametricPlot[Evaluate[{x[t], x'[t]} /. %], {t, 0, i}]

一般WhenEvent StopIntegration来终止计算

你这思路就离谱。全算完再找符合条件的t不就完了？
xsol = NDSolveValue[{x''[t] + x[t] + x[t]^2 == 0, x[0] == 0, x'[0] == 6/10},
x, {t, 0, 10}]
1/2 x'[t]^2 + 1/3 x[t]^3 + 1/2 x[t]^2 /. x -> xsol /. t -> 0
Plot[1/2 x'[t]^2 + 1/3 x[t]^3 + 1/2 x[t]^2 /. x -> xsol, {t, 0, 9.3}, PlotRange -> All]
(* 尽管图像是这个样子，但我并不觉得这是因为能量不守恒了，因为奇点附近的计算精度肯定会受影响。姑且试着算下能量超出1%的位置： *)
FindRoot[1.01 0.18 == 1/2 x'[t]^2 + 1/3 x[t]^3 + 1/2 x[t]^2 /. x -> xsol, {t, 8}]
(* {t -> 9.4571} *)