最近在研究水星进动问题，在最开始的非进动部分，也就是牛顿力学的经典模型下，得到了如图所示的二阶线性非齐次微分方程，通过Odesolve函数求其数值解，其结果及图形如下图：
而实际上，这个图形是可以解出解析解的，其过程如下：
很明显它的解析解是一个椭圆方程
为什么数值解和解析解相差那么多，问题处在了那里？请高手帮我看一下，谢谢！
源文件在这里，版本用的是Mathcad 15 M040版
链接：百度网盘 密码：st5z

对了，上边是解析解的图像。

没有人帮个忙吗？我只是怀疑他在这个方面是有不足还是说我的操作有误啊

右键点击Odesolve，选择“固定”就可以了。
嗯，但得到也不是规则的椭圆。