求助：关于用matlab求解龙格库塔的方法

gwcrow

请问怎么用四阶龙格库塔求解微分方程，用matlab的方法求dV/dt=KV。
谢谢各位高手啦

蜜雪儿

function [x,y]=rk4(dfun,xspan,y0,h)
x=xspan(1):h:xspan(2);
y(1)=y0;
for k=1:length(x)-1
    k1=feval(dfun,x(k),y(k));
    k2=feval(dfun,x(k)+h/2,y(k)+h/2*k1);
    k3=feval(dfun,x(k)+h/2,y(k)+h/2*k2);
    k4=feval(dfun,x(k+1),y(k)+h*k3);
    y(k+1)=y(k)+h*(k1+2*k2+2*k3+k4)/6;
end
x=x';y=y';

gwcrow

万分感谢！！！！:victory: :victory:

happy

http://forum.vibunion.com/forum/viewthread.php?tid=17615

相较二楼的要完善很多

gwcrow

确实很完善，万分感谢！！！！！！！！