用mathematica绘制弹簧摆的运动轨迹

对于如图的弹簧摆，可以写出常微分方程组如下。

这个方程组没有解析解，我们用mathematica求方程组的数值解。代码如下。

摆球的运动轨迹如图

如果弹簧在天花板上的连接点也在做水平运动，运动方程为Asin(w*t) ，则原微分方程变成以下形式。

绘制摆球轨迹如图。