二阶线性非齐次微分方..最近在研究水星进动问题,在最开始的非进动部分,也就是牛顿力学的经典模型下,得到了如图所示的二阶线性非齐次微分方程,通过Odesolve函数求其数值解,其结果及图形如下图:而实际上,这个图形是可以解出解
求微分方程的解析解和其在上的数值解,计算它们之间的误差(用trapz函数计算数值积分),分别用不同的线型和颜色将解析解和数值解绘制在一张图上。计算数值解时,要求步长为0.01。提示:求解析解的数值结果时利用subs命令。trapz函数示例:计算x = 0:0.01:pi;y = sin( ) ...