Running the EpidemiologicalModelsForInfluenzaAndCOVID-19--part_1.nb gives many errors.
I have downloaded the CompartmentalModeling.wl file and put it in the notebook folder and load it inside the notebook also.
Example of the first error:
Running the code:
{susceptibleSEIR, infectedSEIR, exposedSEIR,
recoveredSEIR} = {\[ScriptCapitalS], \[ScriptCapitalI], \
\[ScriptCapitalE], \[ScriptCapitalR]} /.
ParametricNDSolve[
Join[odesSEIR /.
forceOfInfectionSEIR, {\[ScriptCapitalS][
0] == \[ScriptCapitalN] - I0, \[ScriptCapitalE][0] ==
0, \[ScriptCapitalI][0] == I0, \[ScriptCapitalR][0] == 0}],
Head /@ varsSEIR, {t, 0, 100}, {\[ScriptCapitalN],
I0, \[Beta], \[Zeta], \[Gamma]}];
I get the errors:
ParametricNDSolve::dspar: 1.5` cannot be used as a parameter.
ReplaceAll::reps: {ParametricNDSolve[{(\[ScriptCapitalE]^\[Prime])[t]==0. -\[Zeta] \[ScriptCapitalE][t]+1.5 \[ScriptCapitalI][t] \[ScriptCapitalS][t],(\[ScriptCapitalI]^\[Prime])[t]==0. +\[Zeta] \[ScriptCapitalE][t]-0.0478 \[ScriptCapitalI][t],(\[ScriptCapitalR]^\[Prime])[t]==0. +0.0478 \[ScriptCapitalI][t],(\[ScriptCapitalS]^\[Prime])[t]==0. -1.5 \[ScriptCapitalI][t] \[ScriptCapitalS][t],\[ScriptCapitalS][0]==-I0+\[ScriptCapitalN],\[ScriptCapitalE][0]==0,\[ScriptCapitalI][0]==I0,\[ScriptCapitalR][0]==0},{\[ScriptCapitalE],\[ScriptCapitalI],\[ScriptCapitalR],\[ScriptCapitalS]},{t,0,100},{\[ScriptCapitalN],I0,1.5,\[Zeta],0.0478}]} is neither a list of replacement rules nor a valid dispatch table, and so cannot be used for replacing.
Set::shape: Lists {susceptibleSEIR,infectedSEIR,exposedSEIR,recoveredSEIR} and {\[ScriptCapitalS],\[ScriptCapitalI],\[ScriptCapitalE],\[ScriptCapitalR]}/. ParametricNDSolve[{(\[ScriptCapitalE]^\[Prime])[t]==0. -\[Zeta] \[ScriptCapitalE][t]+1.5 \[ScriptCapitalI][t] \[ScriptCapitalS][t],(\[ScriptCapitalI]^\[Prime])[t]==0. +\[Zeta] \[ScriptCapitalE][t]-0.0478 \[ScriptCapitalI][t],(\[ScriptCapitalR]^\[Prime])[t]==0. +0.0478 \[ScriptCapitalI][t],(\[ScriptCapitalS]^\[Prime])[t]==0. -1.5 \[ScriptCapitalI][t] \[ScriptCapitalS][t],\[ScriptCapitalS][0]==-I0+\[ScriptCapitalN],\[ScriptCapitalE][0]==0,\[ScriptCapitalI][0]==I0,\[ScriptCapitalR][0]==0},{\[ScriptCapitalE],\[ScriptCapitalI],\[ScriptCapitalR],\[ScriptCapitalS]},{t,0,100},{\[ScriptCapitalN],I0,1.5,\[Zeta],0.0478}] are not the same shape.
