This is a simple SIR model, built to reproduce the one that we made in InsightMaker.

Andy Long
Spring, 2018

Starting parameters from InsightMaker model

Clear["Global`*"]

RecoveryRate=0.1;
InfectionRate=0.001;
ResusceptibleRate=0.01;
infected=20;
pop=100;

Manipulate[
lv=
NDSolve[
{(*  susceptible: *)
sus'[t]==ResusceptibleRate  rec[t]- InfectionRate sus[t] inf[t],
(* infected: *)
inf'[t]==InfectionRate sus[t] inf[t]-RecoveryRate inf[t],
(* recovered: *)
rec'[t]==RecoveryRate inf[t]-ResusceptibleRate  rec[t],
sus[0]==pop,
inf[0]==infected,
rec[0]==0
},
{sus[t],inf[t],rec[t]},
{t,0,tend}];
Plot[Evaluate[{sus[t],inf[t],rec[t]}/.lv],{t,0,tend},
PlotRange->All,AxesOrigin->{0,0}],
{{lv,{}},None},
{{tend,50},2,2000,10},
{{RecoveryRate,0.09},0.01,.7,.01},
{{InfectionRate,0.008},0.00001,.1,.001},
{{ResusceptibleRate,0.027},0.00001,.1,.001}
]