function dx=modeloSIR(t,x)
beta=5e-4;  % beta=1e-4 --> aislamiento social;
gamma=0.05;
N=1000;     % Poblacion total previa a la infección
dx=[-beta*x(1)*x(2);
    beta*x(1)*x(2)-gamma*x(2);
    gamma*x(2)];
end