微分、差分方程习题.doc

微分方程和差分方程作业题参考答案 一、微分方程初值问题 （1）用四阶 Runge-Kutta 法求解微分方程初值问题的数值解(步长h取0.1)，求解范围为区间[0,3]． （2）用 ode45 方法常微分方程初值问题的数值解(近似解)，然后利用画图来比较两者间的差异． 解 （1）代码 clear f=sym(y-exp(x)*cos(x)); a=0; b=3; h=0.1; n=(b-a)/h+1; % n=(b-a)/h; x=0; y=1; szj=[x,y]; for i=1:n-1 % i=1:n l1=subs(f,{x,y l2=subs(f,{x,y