四阶龙格库塔法的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(i,:)+h/6*(k1+2*k2+2*k3+k4);
    i=i+1;
end

在使用该代码前，需要自己定义函数f，其表示微分方程右侧的函数。其他的参数包括：
- `tspan`：数组，表示计算的时间范围；
- `y0`：数组，表示初始值；
- `h`：步长，表示每次迭代的时间间隔。

调用该函数时，可以得到计算结果，其中`t`为时间序列，`y`为该微分方程的解。