SIR模型相軌線如何畫？matlab代碼

i0?=?0.02;

s0?=?0.98;

x0=[i0;s0];

[t,x]=ode45(@sirfun,ts,x0);

figure;?plot(x(:,2),x(:,1));grid?on;

function?dy?=?sirfun(t,x)

%?x(1)-->?i,?x(2)-->s

%?di/dt?=?λsi-μi

%?ds/dt?=?-λsi

lambda?=?1;

mu?=?0.3;

dy?=?[lambda*x(1)*x(2)-mu*x(1);

-lambda*x(1)*x(2)];

end