Math ParametericPlot上没有输出

我正在用Mathematica的

NDSolve

求解和绘制双摆的运动方程

我已经使用标准绘图成功地绘制了角度位置。但是当我开始用参数图来表示每个物体的位置时。我没有错误，只是没有情节

eqn1 = 2 th''[t] + Sin[th[t] - ph[t]] (ph'[t])^2 + Cos[th[t] - ph[t]] (ph''[t]) + (2 g/l) Sin[th[t]]
eqn2 = ph''[t] + Sin[th[t] - ph[t]] (th'[t])^2 + Cos[th[t] - ph[t]] (th''[t]) + (g/l) Sin[th[t]]
eqnA = eqn1 /. {g -> 10, l -> 1}
eqnB = eqn2 /. {g -> 10, l -> 1}

sol = NDSolve[{eqnA == 0, eqnB == 0, th[0] == 0.859, th'[0] == 0, ph[0] == 0.437, ph'[0] == 0}, {th, ph}, {t, 0, 10}]
Plot[{th[t], ph[t]} /. sol, {t, 0, 10}]
r1 = {lSin[th[t]] + lSin[ph[t]], -lCos[th[t]] - lCos[ph[t]]} /. {l -> 1, g -> 10}

ParametricPlot[r1 /. sol, {t, 0, 10}]

替换

r1 = {lSin[th[t]] + lSin[ph[t]], -lCos[th[t]] - lCos[ph[t]]} /. {l->1, g->10}

与

您的

参数图应该出现

您可能还记得一个有用的技巧，当任何绘图都没有出现时，您可以尝试用
表
替换绘图，并查看它显示了什么。通常情况下，数据表会提供所需的提示，说明为什么绘图没有出现。
我投票将这个问题作为离题题来结束，因为它属于主题
r1 = {l*Sin[th[t]] + l*Sin[ph[t]], -l*Cos[th[t]] - l*Cos[ph[t]]} /. {l->1, g->10}