用Mathematica求常微分方程数值解？

Question

用Mathematica求解数值解，碰到了如下问题：

自觉没问题，但是运行不出来，请各位帮忙看看。

Answer 1

用Mathematica求常微分方程数值解的一般步骤是：

调用NDSolve函数，输入微分方程、初始条件、未知函数和自变量的范围。
从NDSolve的输出中提取插值函数作为数值解。
使用Plot或ParametricPlot等函数绘制数值解的图像。

sol = NDSolve[{theta''[t] == -Q*Sin[theta[t]], theta[0] == 0, theta'[0] == u0}, theta, {t, 0, 20}];

现在，我们可以定义一个新的函数 thetaSol 来代表解：
thetaSol[t_] = theta[t] /. First[sol];
你可以用 thetaSol 函数计算任何给定时刻 t 的数值解：
thetaSol[3]
最后，使用 Plot 函数绘制解的图形：
Plot[thetaSol[t], {t, 0, 20}]
这将绘制出 t 从 0 到 20 的数值解曲线。