Matlab - Modelo de epidemias SIR con ecuaciones diferenciales

function epidemia

beta=.03;

gamma=2;

x0=100;

y0=10;

z0=0;

tmax=6;

[t,x]=ode45(@epi,[0 tmax],[x0,y0,z0]);

subplot(3,1,1),plot(t,x(:,1)),title('Susceptibles')

subplot(3,1,2),plot(t,x(:,2)),title('Infectados')

subplot(3,1,3),plot(t,x(:,3)),title('Recobrados')

    function dx = epi(~,x)

        dx = zeros(3,1);

        dx(1) = -beta*x(1)*x(2);

        dx(2) = +beta*x(1)*x(2)-gamma*x(2);

        dx(3) = +gamma*x(2);

    end

end