计算物理作业(龙格库塔四阶求解初值)

《计算物理作业(龙格库塔四阶求解初值)》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、EASTCHINAINSTITUTEOFTECHNOLOGY计算物理作业学生学号：2014203701学生姓名：专业：物理学班级:14237012017年1月6日-一编写程序求解微分方程用四阶龙格一库塔法求解微分方程的初值问题dydxy

2、，我们求解微分方程的方法是四阶Runge-Kutta法y(k+1)=y(k)+h*(k1+2*k2+2*k3+k4)/6;kl=f(x(k),y(k));k2=f(x(k)+0.5*h,y(k)+0.5*h*k1);k3=f(x(k)+0.5*h,y(k)+0.5*h*k2);k4=f(x(k)+h,y(k)+h*k3);采用这种方法计算的结果一般比较精确，并且编写的程序具有占用内存小的优点。源代码:functionSDeqRunKu4(dfun,a,b,yO,Nm)h=(b-a)/Nm;x=a:h:b;y二z

3、eros(1,Nm+1);y(l)二y0；fori=1:Nmkl=feval(dfun,x(i),y(i));k2=feval(dfun,x(i)+h/2,y(i)+h*kl/2);k3=feval(dfun,x(i)+h/2,y(i)+h*k2/2);k4二feval(dfun,x(i)+h,y(i)+h*k3);y(i+l)=y(i)+(kl+2*k2+2*k3+k4)*h/6;endS二[x‘y'];运行的结果:当h=0・06时，运行的结果如下图所示:11ir>lriiAkvm命令葩口不熟慈MATLAB?

4、活参闫有关快速入门的资源。»x=2:0.062.6;>>dfun=inlineCx/y'J>JJ,x,y)：>>S=DeqRunKu4(dfun,2,2.6,1,10);»[S]s=2.00001.00002.06001.11522.12001.22252.18001.32382.24001.42042.30001.51332.36001.60302.42001.69012.48001.77492.54001.8578打2.60001.9391误差分析与绘图当h=0・06时,运行的结果:目Graphl