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
