I have the next code:
sirds[[Alpha], [Beta], [Delta], [Mu]] := {S[t], SS[t], I[t], R[t],
D[t]} /.
First@NDSolve[{
S'[t] == -[Alpha]*S[t]*I[t] - [Delta]*SS[t],
I'[t] == [Alpha]*S[t]*I[t] - [Beta]*I[t] - [Mu]*I[t],
R'[t] == [Beta]*I[t],
SS'[t] = [Delta]*SS[t],
D'[t] == [Mu]*I[t], S[0] == 911000, I[0] == 22894, R[0] = 0,
D[0] == 196,
SS[0] = 0}, {[Alpha], [Beta], [Delta], [Mu]}, {t, 365}]
Manipulate[Plot[
Evaluate[sirds[[Alpha], [Beta], [Delta], [Mu]]],
{t, 0, 36},
PlotLegends -> {"Susceptible", "Vacunada", "Infectada",
"Recuperada", "Muerta"},
Frame -> True, PlotRange -> {0, 50000},
PlotStyle -> {
{Blue, Thickness[0.005]},
{Red, Thickness[0.005]},
{Green, Thickness[0.005]},
{Black, Thickness[0.005]},
{Yellow, Thickness[0.005]}}],
{{[Alpha], .017428, "Infectividad"}, 0, .2},
{{[Beta], .3, "Recuperación"}, 0, .8},
{{[Delta], .0213577, "Letalidad"}, 0, .2},
{{[Mu], .002, "Eficacia"}, 0, .2}]
When I run it, it makes the table that I want, but it doesn't show the graphics that should show. Someone can help me? I took the same structure of this original code that runs perfectly:
sir[A_, B_, P_ ] := {a[t], f[t], r[t], d[t]} /.
First@NDSolve[{
a'[t] == -Aa[t] f[t],
r'[t] == Bf[t],
d'[t] == Pf[t],
f'[t] == Aa[t] f[t] - Bf[t] - Pf[t],
a[0] == 478530 - 36, f[0] == 36, r[0] == 0, d[0] == 0}, {a, f, r,
d}, {t, 112}]
Manipulate[Plot[
Evaluate[sir[A, B, P]],
{t, 0, 112},
PlotLegends -> {"Susceptible", "Infectada", "Recuperada", "Muerta"},
Frame -> True, PlotRange -> {0, 478530}, PlotStyle -> {
{Blue, Thickness[0.005]},
{Red, Thickness[0.005]},
{Green, Thickness[0.005]},
{Black, Thickness[0.005]}}],
{{A, .00000246, "Infectividad"}, 0, 2},
{{B, .69, "Recuperación"}, 0, 2},
{{P, .0111, "Letalidad"}, 0, 2}]