求助，用 MATLAB 的 ode45 求解微分方程组

Matlab软件中如何使用ode 45 解微分方程组？求解过程应按下列步骤来计算：

第一步：根据已经微分方程组和相关系数，自定义求解微分方程组的函数，其函数名 odefun，其参数为【t,z】

第二步：由于未知初始条件，用随机数初定，即z0=rand(1,6)/1000

第三步：确定时间t的范围，如tspan=[0 50]

第四步：使用ode45函数，求其数值解，即

[t,z]=ode45(@(t,z)odefun(t,z),tspan,z0)

第五步：使用plot函数，绘制其相关图形。

第六步：编程后执行，得到结果。

图1，耦合14－16方程组的结果。

ode45 求解非刚性微分方程，中阶方法。

ode45可以解决问题 M(t,y)*y' = f(t,y) ，质量矩阵M是非奇异的。

Example

解决系统 y' = vdp1(t,y)

使用默认的相对误差容差1e-3和每个组件的默认绝对容差1e-6，并绘制解决方案的第一个组件。