Runge-Kutta 算法

Runge-Kutta 4阶算法是一种常用的数值计算方法，用于求解常微分方程组。在Matlab中，可以使用ode45函数来实现该算法。该函数可以自动选择合适的步长，并返回求解结果。具体使用方法如下：

1. 定义常微分方程组：

function dydt = myode(t,y)
dydt = [y(2); -y(1)];

2. 调用ode45函数：

[t,y] = ode45(@myode,[ 10],[1 ]);

其中，@myode表示将myode函数作为参数传递给ode45函数，[ 10]表示求解区间为到10，[1 ]表示初始条件为y()=1，y'()=。

3. 绘制结果：

plot(t,y(:,1),'-o',t,y(:,2),'-x');

其中，y(:,1)表示y的第一列，即y(1)，y(:,2)表示y的第二列，即y(2)。