不知道咋解了
m*g*l*Sin[\[Theta][t]]/2 + \[Delta]*
Sin[\[Theta][t] -
arctan[l*Sin[\[CapitalTheta][t]]/[l*Cos[\[Theta][t]] + d]]]*
Sin [ArcTan[l*Sin[\[Theta][t]]/[l*Cos[\[Theta][t]] + d]]]^2/[l*
sin[\[Theta][t]]^2] - \[Delta]*
Sin[\[Theta][t] -
ArcTan[l*Sin[\[Theta][t]]/[l*Cos[\[Theta][t]] + d + h]]]*
Sin[ArcTan[l*Sin[\[Theta][t]]/[l*Cos[\[Theta][t]] + d]]]^2/[l*
Sin[\[Theta][t]]^2] == m*l^2*D[\[Theta][t], {t, 2}]/3
NDSolve[eqn, \[Theta][0] == Pi/2]

微分方程 : dv切 = v切 cos dθ - v径 sin dθ - v切 这是 二体问题 的 一个 微分方程, 也可以说是 一体问题 的 一个 微分方程 。