¿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)-->?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)];

fin