课程设计报告-四阶Runge-Kutta方法.docVIP

《计算机数值方法》 课程设计报告 题目 四阶Runge-Kutta方法 学生姓名 班级 成绩 指导教师 年月日 四、计算公式或算法 输入（编写或调用计算的函数文件）， 3．For End 4．输出 五、Matlab 程序 x=[a:h:b]; y(1)=y1; n=(b-a)/h+1; for i=2:n fk1=f(x(i-1),y(i-1)); fk2=f(x(i-1)+h/2,y(i-1)+fk1*h/2); fk3=f(x(i-1)+h/2,y(i-1)+fk2*h/2); fk4=f(x(i-1)+h,y(i-1)+fk3*h); y(i)=y(i-1)+h*(fk1+2*fk2+2*fk3+fk4)/6; end y 六、测试数据及结果 用调试好的程序解决如下问题： 应用经典的四阶Runge-Kutta方法解初值问题 取 步骤一：编写函数具体程序. 1.求解解析解程序： dsolve(Dy=(y^2+y)/t,y(1)=-2,t) 结果： 2.综合编写程序如下： a=1; b=3; h=0.5; y(1)=-2; x(1)=a; n=(b-a)/h+1; yy(1)=-2; for i=2:n