四阶龙格库塔法的matlab程序

四阶龙格库塔法是一种数值积分方法，用于解决常微分方程。下面是一个简单的Matlab代码，实现了四阶龙格库塔法：

function [t,y]=rk4(f,tspan,y0,h)
t(1)=tspan(1);
y(1,:)=y0;
i=1;
while t(i)<tspan(2)
    k1=f(t(i),y(i,:));
    k2=f(t(i)+h/2,y(i,:)+h/2*k1);
    k3=f(t(i)+h/2,y(i,:)+h/2*k2);
    k4=f(t(i)+h,y(i,:)+h*k3);
    t(i+1)=t(i)+h;
    y(i+1,:)=y(