matlab中ode23的基本使用方法（秒懂）

题目如下图

解题过程以及注释如下

function [y1,y2]= myfun3
tspan = 0:0.1:20;%计算结果的数值区间
y0 =[20;12];%y1和y2的初值，即y1(0)和y2(0)
[t,r]= ode23(@dy, tspan,y0);%调用ode23函数
%注意，这里使用的调用函数的方法是@函数名，而非'函数名'，因为该函数并非.m文件，不能用引号调用。
y1 = r(:,1);%结果赋值
y2 = r(:,2);%结果赋值

function z =dy(t,y)%因为是方程组，必须用矩阵哦
z = [(1 - 0.013*y(2))*y(1);%存放y1't
    (-1 + 0.011*y(1))*y(2)];%存放y2't