(Subscript[R, e])=2.6;
(Subscript[R, max])=5.7;
(Subscript[R, min])=1.3;
(Subscript[R, 0])=(Subscript[R, min]);
ν=5;
(Subscript[γ, 1])=1/14;
(Subscript[γ, 2])=1/14;
(Subscript[γ, 3])=1/14;
(Subscript[δ, 1])=1/5;
(Subscript[δ, 2])=1/5;
(Subscript[δ, 3])=1/5;
(Subscript[β, 1])=(Subscript[R, 0])(Subscript[δ, 1])/ν;
(Subscript[β, 2])=(Subscript[R, 0])(Subscript[δ, 2])/(ν^2);
(Subscript[β, 3])=(Subscript[R, 0])(Subscript[δ, 3])/(ν^2);
(Subscript[α, 1])=(Subscript[β, 1])/(ν^2);
(Subscript[α, 2])=(Subscript[β, 2])/(ν^2);
(Subscript[α, 3])=(Subscript[β, 3])/(ν^2);
(Subscript[κ, 1])=(Subscript[β, 1]);
(Subscript[κ, 2])=(Subscript[β, 2]);
(Subscript[κ, 3])=(Subscript[β, 3]);
(Subscript[μ, 1])=0.002;
(Subscript[μ, 2])=0.2;
(Subscript[μ, 3])=0.001;
M=18000000;
(Subscript[p, 1])=0.4;
(Subscript[p, 2])=0.05;
(Subscript[p, 3])=1-(Subscript[p, 1])-(Subscript[p, 2]);
(Subscript[n, 1])=(Subscript[p, 1])M;
(Subscript[n, 2])=(Subscript[p, 2])M;
(Subscript[n, 3])=(Subscript[p, 3])M;

system={
(Subscript[s, 1])'[t]==-(Subscript[β, 1]) (((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[s, 1])[t]-(Subscript[κ, 1]) (((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[s, 1])[t]+(Subscript[κ, 1])(Subscript[γ, 1]) (Subscript[c, 1])[t],
(Subscript[c, 1])'[t]==(Subscript[κ, 1]) (((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[s, 1])[t]-(Subscript[κ, 1])(Subscript[γ, 1]) (Subscript[c, 1])[t]-(Subscript[α, 1]) (((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[c, 1])[t],
(Subscript[i, 1])'[t]== (Subscript[β, 1]) (((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[s, 1])[t]-(Subscript[γ, 1])(1-(Subscript[μ, 1]))(Subscript[i, 1])[t]-(Subscript[γ, 1])(Subscript[μ, 1]) (Subscript[i, 1])[t]+(Subscript[α, 1]) (((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[c, 1])[t],
(Subscript[r, 1])'[t]== (Subscript[γ, 1])(1-(Subscript[μ, 1]))(Subscript[i, 1])[t],
(Subscript[d, 1])'[t]==(Subscript[γ, 1])(Subscript[μ, 1]) (Subscript[i, 1])[t],
(Subscript[s, 2])'[t]==-(Subscript[β, 2]) ((((Subscript[i, 1])[t])/(Subscript[n, 1]))) (Subscript[s, 2])[t]-(Subscript[κ, 2]) (((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[s, 2])[t]+(Subscript[κ, 2])(Subscript[γ, 2]) (Subscript[c, 2])[t],
(Subscript[c, 2])'[t]==(Subscript[κ, 2]) (((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[s, 2])[t]-(Subscript[κ, 2])(Subscript[γ, 2]) (Subscript[c, 2])[t]-(Subscript[α, 2]) (((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[c, 2])[t],
(Subscript[i, 2])'[t]== (Subscript[β, 2]) ((((Subscript[i, 1])[t])/(Subscript[n, 1]))) (Subscript[s, 2])[t]+(Subscript[α, 2])(((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[c, 2])[t]-(Subscript[γ, 2])(1-(Subscript[μ, 2]))(Subscript[i, 2])[t]-(Subscript[γ, 2])(Subscript[μ, 2])(Subscript[i, 2])[t],
(Subscript[r, 2])'[t]== (Subscript[γ, 2])(1-(Subscript[μ, 2]))(Subscript[i, 2])[t],
(Subscript[d, 2])'[t]==(Subscript[γ, 2])(Subscript[μ, 2])(Subscript[i, 2])[t],
(Subscript[s, 3])'[t]==-(Subscript[β, 3]) ((((Subscript[i, 1])[t])/(Subscript[n, 1]))) (Subscript[s, 3])[t]-(Subscript[κ, 3]) (((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[s, 3])[t]+(Subscript[κ, 3])(Subscript[γ, 3]) (Subscript[c, 3])[t],
(Subscript[c, 3])'[t]==(Subscript[κ, 3]) (((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[s, 3])[t]-(Subscript[κ, 3])(Subscript[γ, 3]) (Subscript[c, 3])[t]-(Subscript[α, 3]) (((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[c, 3])[t],
(Subscript[i, 3])'[t]== (Subscript[β, 3]) ((((Subscript[i, 1])[t])/(Subscript[n, 1]))) (Subscript[s, 3])[t]+(Subscript[α, 3])(((Subscript[i, 1])[t])/(Subscript[n, 1])) (Subscript[c, 3])[t]-(Subscript[γ, 3])(1-(Subscript[μ, 3]))(Subscript[i, 3])[t]-(Subscript[γ, 3])(Subscript[μ, 3])(Subscript[i, 3])[t],
(Subscript[r, 3])'[t]== (Subscript[γ, 3])(1-(Subscript[μ, 3]))(Subscript[i, 3])[t],
(Subscript[d, 3])'[t]==(Subscript[γ, 3])(Subscript[μ, 3])(Subscript[i, 3])[t],
(Subscript[c, 1])[0]==10000,
(Subscript[i, 1])[0]==196,
(Subscript[r, 1])[0]==20,
(Subscript[d, 1])[0]==5,
(Subscript[s, 1])[0]==(Subscript[n, 1])-(Subscript[c, 1])[0]-(Subscript[i, 1])[0]-(Subscript[r, 1])[0]-(Subscript[d, 1])[0],
(Subscript[c, 2])[0]==500,
(Subscript[i, 2])[0]==0,
(Subscript[r, 2])[0]==0,
(Subscript[d, 2])[0]==0,
(Subscript[s, 2])[0]==(Subscript[n, 2])-(Subscript[c, 2])[0]-(Subscript[i, 2])[0]-(Subscript[r, 2])[0]-(Subscript[d, 2])[0],
(Subscript[c, 3])[0]==0,
(Subscript[i, 3])[0]==0,
(Subscript[r, 3])[0]==0,
(Subscript[d, 3])[0]==0,
(Subscript[s, 3])[0]==(Subscript[n, 3])-(Subscript[c, 3])[0]-(Subscript[i, 3])[0]-(Subscript[r, 3])[0]-(Subscript[d, 3])[0]
};

solmin=NDSolve[system,{(Subscript[s, 1]),(Subscript[c, 1]),(Subscript[i, 1]),(Subscript[r, 1]), (Subscript[d, 1]),(Subscript[s, 2]),(Subscript[c, 2]),(Subscript[i, 2]),(Subscript[r, 2]), (Subscript[d, 2]),(Subscript[s, 3]),(Subscript[c, 3]),(Subscript[i, 3]),(Subscript[r, 3]), (Subscript[d, 3])},{t,0,500}];

Plot[Evaluate[{(Subscript[i, 1])[t]+(Subscript[i, 2])[t]+(Subscript[i, 3])[t],(Subscript[d, 1])[t]+(Subscript[d, 2])[t]+(Subscript[d, 3])[t]}/.sol],{t,0,500}, PlotRange->All, PlotLabels->{"Infectados","Muertes"}]

Evaluate[{(Subscript[i, 1])[t]+(Subscript[i, 2])[t]+(Subscript[i, 3])[t],(Subscript[d, 1])[t]+(Subscript[d, 2])[t]+(Subscript[d, 3])[t]}/.sol]/.t->300

Plot[Evaluate[{(Subscript[i, 1])[t]+(Subscript[i, 2])[t]+(Subscript[i, 3])[t],(Subscript[d, 1])[t]+(Subscript[d, 2])[t]+(Subscript[d, 3])[t]}/.solmax],{t,0,500}, PlotRange->All, PlotLabels->{"Infectados","Muertes"}]

Plot[Evaluate[{(Subscript[i, 1])[t]+(Subscript[i, 2])[t]+(Subscript[i, 3])[t],(Subscript[d, 1])[t]+(Subscript[d, 2])[t]+(Subscript[d, 3])[t]}/.solmin],{t,0,500}, PlotRange->All, PlotLabels->{"Infectados","Muertes"}]

Plot[{Evaluate[(((Subscript[i, 1])[t]+(Subscript[i, 2])[t]+(Subscript[i, 3])[t])/M)100/.solmax],Evaluate[(((Subscript[i, 1])[t]+(Subscript[i, 2])[t]+(Subscript[i, 3])[t])/M)100/.sol],Evaluate[(((Subscript[i, 1])[t]+(Subscript[i, 2])[t]+(Subscript[i, 3])[t])/M)100/.solmin]},{t,0,500},PlotRange->All, PlotLabels->Placed[{"Máximo","Estimado","Mínimo"},Above], PlotLabel->"Diferentes Escenarios \n Infectados en Guatemala", AxesLabel->{"Días","(%)"}, AxesOrigin->{0,0}]

Plot[{Evaluate[(((Subscript[d, 1])[t]+(Subscript[d, 2])[t]+(Subscript[d, 3])[t])/M)100/.solmax],Evaluate[(((Subscript[d, 1])[t]+(Subscript[d, 2])[t]+(Subscript[d, 3])[t])/M)100/.sol],Evaluate[(((Subscript[d, 1])[t]+(Subscript[d, 2])[t]+(Subscript[d, 3])[t])/M)100/.solmin]},{t,0,500},PlotRange->All, PlotLabels->Placed[{"Máximo","Estimado","Mínimo"},Above], PlotLabel->"Diferentes Escenarios \n Muertes en Guatemala", AxesLabel->{"Días","(%)"}]

Evaluate[{(Subscript[d, 1])[t]+(Subscript[d, 2])[t]+(Subscript[d, 3])[t]}/.solmax]/.t->300

Evaluate[{(Subscript[d, 1])[t]+(Subscript[d, 2])[t]+(Subscript[d, 3])[t]}/.solmin]/.t->300

Evaluate[{(Subscript[d, 1])[t]+(Subscript[d, 2])[t]+(Subscript[d, 3])[t]}/.sol]/.t->300