鏁板�煎垎鏋愪綔涓歁ATLAB绋嬪簭棰?- 鐧惧害鏂囧簱

鏂板缓

涓婁紶

棣栭〉

鍔╂墜

鏈�杩?/div>

璧勬枡搴?/div>

宸ュ叿

 

function

 [x,y1,y2]=my(dyfun,xspan,y0,h)

 

x=xspan(1):h:xspan(2);

 

y1(1)=y0;

 

y2(1)=y0;

 

for

 n=1:length(x)-1

 

y1(n+1)=y1(n)+h*feval(dyfun,x(n),y1(n));

 

k1=feval(dyfun,x(n),y2(n));

 

k2=feval(dyfun,x(n)+h/2,y2(n)+h/2*k1);

 

k3=feval(dyfun,x(n)+h/2,y2(n)+h/2*k2);

 

k4=feval(dyfun,x(n)+h,y2(n)+h*k3);

 

y2(n+1)=y2(n)+h/6*(k1+2*k3+2*k2+k4);

 

end