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，并绘制解决方案的第一个组件。
