常微分方程数值解实验报告2.docVIP

实验目的： 1．熟练编程技巧，体会Runge-kutta方法。 2．与前面的实验方法进行对比，了解数值计算方法区别和联系。 实验内容： 运用Runge-kutta方法求解常微分方程初值问题 实验题目：，。 实验原理： 取四阶Runge-kutta方法进行计算 实验程序： %fun为目标函数字符串 %x0为自变量初始值。 %y0为fun(x0); %bou=[a,b]自变量区间 %h为步长 fun=x*y+x; bou=[3,6]; a=bou(1); b=bou(2); x0=3; y0=8; h=0.1; n=ceil((b-a)/h); xx=linspace(a,b,n+1); yy=zeros(1,n+1); lengthx=length(xx); xx(1)=x0;yy(1)=y0; for i=2:n+1 x=xx(i-1);y=yy(i-1); k1=eval(fun); x=xx(i-1)+h/2;y=yy(i-1)+h*k1/2; k2=eval(fun); x=xx(i-1)+h/2;y=yy(i-1)+h*k2/2; k3=eval(fun); x=xx(i-1)+h;y=yy(i-1)+h*k3; k4=eval(fun); yy(i)=yy(i-1)+h/6*(k1+