${\textstyle S(0)=997,I(0)=3,R(0)=0}$

${\textstyle \beta =0,0005}$

${\textstyle \gamma =0,035}$

${\textstyle \mu =0,005}$

${\textstyle \beta =0,0005}$

${\textstyle \beta =0,00012}$

${\displaystyle {\frac {\mathrm {d} S}{\mathrm {d} t}}=-\beta \cdot S(t)\cdot I(t)}$

${\displaystyle {\frac {\mathrm {d} I}{\mathrm {d} t}}=\beta \cdot S(t)\cdot I(t)-\gamma \cdot I(t)-\mu \cdot I(t)}$

${\displaystyle {\frac {\mathrm {d} R}{\mathrm {d} t}}=\gamma \cdot I(t)}$

${\displaystyle {\frac {\mathrm {d} D}{\mathrm {d} t}}=\mu \cdot I(t)}$

Diese Animation wurde mit de:GeoGebra erstellt, wobei die numerischen Lösungen des genannten Systems der Differentialgleichungen ermittelt wurden. Der entscheidende Teil der Konstruktion lautet:

##########################################
# das System der Differentialgleichungen #
##########################################
S'(t, S, I, R, D) = -ß I S
I'(t, S, I, R, D) = ß S I - γ I - μ I
R'(t, S, I, R, D) = γ I
D'(t, S, I, R, D) = μ I

#######################################################################
# numerische Lösung des Systems der Differentialgleichungen bestimmen #
#######################################################################
NLöseDgl[{S', I', R', D'}, 0, {s_0, i_0, r_0, d_0}, T_{max}]

# liefert:
# NumerischesIntegral1 -> S
# NumerischesIntegral2 -> I
# NumerischesIntegral3 -> R
# NumerischesIntegral4 -> D

Anmerkung: In der Literatur werden teilweise modifizierte Formen dieser DGL benutzt, die aber gleichwertig sind. Leider werden oft dieselben Parameter benutzt. Beispiel für die erste DGL, wobei die Infektionsrate der Klarheit halber ${\textstyle {\tilde {\beta }}}$ benannt ist:

${\displaystyle {\frac {\mathrm {d} S}{\mathrm {d} t}}=-{\frac {\tilde {\beta }}{N}}\cdot S(t)\cdot I(t)}$

${\textstyle \beta ={\frac {\tilde {\beta }}{N}}}$

${\textstyle N}$

${\textstyle N}$