kill(a,b,c,d,e,x,A,B,C,S,I,R)$
a: 2.6$
b : 1$
epsilon: 0.1$
S0 : 0.34+epsilon$
I0 : 0.01$
R0 : 0.65-epsilon$
SIR: rk([ -a*S*I,a*S*I-b*I, b*I],[S,I,R],[S0,I0,R0],[t,0,10,0.1])$
SIRs: makelist([f[1],f[2]],f,SIR)$
SIRi: makelist([f[1],f[3]],f,SIR)$
SIRr: makelist([f[1],f[4]],f,SIR)$
wxplot2d([[discrete,SIRs],[discrete,SIRi],[discrete,SIRr]],[legend,"S","I","R"],[color,green,red,black]);
wxplot2d([[discrete,SIRs],[discrete,SIRi],[discrete,SIRr]],[legend,"S","I","R"],[color,green,red,black],[y,0,0.05]);

# Influenza

kill(a,b,c,d,e,x,A,B,C,S,I,R)$
a: 1.3$
b : 1$
epsilon: 0.1$
S0 : 0.76+epsilon$
I0 : 0.01$
R0 : 0.23-epsilon$
SIR: rk([ -a*S*I,a*S*I-b*I, b*I],[S,I,R],[S0,I0,R0],[t,0,10,0.1])$
SIRs: makelist([f[1],f[2]],f,SIR)$
SIRi: makelist([f[1],f[3]],f,SIR)$
SIRr: makelist([f[1],f[4]],f,SIR)$
wxplot2d([[discrete,SIRs],[discrete,SIRi],[discrete,SIRr]],[legend,"S","I","R"],[color,green,red,black]);
wxplot2d([[discrete,SIRs],[discrete,SIRi],[discrete,SIRr]],[legend,"S","I","R"],[color,green,red,black],[y,0,0.05]);

# Boginje (variola vera)

kill(a,b,c,d,e,x,A,B,C,S,I,R)$
a: 3.5$
b : 1$
epsilon: 0.1$
S0 : 0.28+epsilon$
I0 : 0.01$
R0 : 0.71-epsilon$
SIR: rk([ -a*S*I,a*S*I-b*I, b*I],[S,I,R],[S0,I0,R0],[t,0,10,0.1])$
SIRs: makelist([f[1],f[2]],f,SIR)$
SIRi: makelist([f[1],f[3]],f,SIR)$
SIRr: makelist([f[1],f[4]],f,SIR)$
wxplot2d([[discrete,SIRs],[discrete,SIRi],[discrete,SIRr]],[legend,"S","I","R"],[color,green,red,black]);
wxplot2d([[discrete,SIRs],[discrete,SIRi],[discrete,SIRr]],[legend,"S","I","R"],[color,green,red,black],[y,0,0.05]);

# Ospice

kill(a,b,c,d,e,x,A,B,C,S,I,R)$
a: 15$
b : 1$
epsilon: 0.02$
S0 : 0.06+epsilon$
I0 : 0.01$
R0 : 0.93-epsilon$
SIR: rk([ -a*S*I,a*S*I-b*I, b*I],[S,I,R],[S0,I0,R0],[t,0,10,0.1])$
SIRs: makelist([f[1],f[2]],f,SIR)$
SIRi: makelist([f[1],f[3]],f,SIR)$
SIRr: makelist([f[1],f[4]],f,SIR)$
wxplot2d([[discrete,SIRs],[discrete,SIRi],[discrete,SIRr]],[legend,"S","I","R"],[color,green,red,black]);
wxplot2d([[discrete,SIRs],[discrete,SIRi],[discrete,SIRr]],[legend,"S","I","R"],[color,green,red,black],[y,0,0.05]);