Modelo de epidemias SIR con ecuaciones diferenciales Nacimientos y muertes
Publicado por Martha B. (3 intervenciones) el 30/11/2015 08:53:24
Buenas noches, el mes pasado pedí ayuda con un código para modelar el modelo SIR de epidemias y Daniel me ayudo con el código de aquí abajo sin embargo, no se como a nadir las muertes y nacimientos,me podrían ayudar de nuevo?
Gracias por su atención
Gracias por su atención
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
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
Valora esta pregunta
1