¿Cómo dibujar la trayectoria de fase del modelo SIR? código matlab

ts?=?0:50;

i0?=?0.02;

s0?=?0.98;

x0= [i0;s0];

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

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

función?dy?=?sirfun(t,x)

x(1)--gt;?i,?x (2)--gt;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)];

fin