数学建模第4章.docVIP

PAGE PAGE 23 4.1 欧拉方法和龙格－库塔方法求微分方程数值解，画出解的图形，对结果进行比较分析。以下方程供选择： d) 幂级数解 解： 原方程化为 记，并用t代替x，则原方程化为： ；且； 于是可以用龙格－库塔法求解。 【模型求解】： 用Matlab作龙格－库塔法求解： % chapter 4—1.d %此函数是微分方程组 function Xdot=ch41dfun(t,x) Xdot=[x(2),-cos(t)*x(1)]; %ch41d.m function I=ch41d(a) x0=[1,0]; [t,x]=ode45(ch41dfun,[0,a],x0); y