I’m trying to use DiffEqOperators to define the Kermack-McKendrick model, which is basically a PDE generalization of the simple SIR epidemiological model:
The equation for I is a PDE, while that for S and R is an ODE with an integral of I over a at t.
The declarations seem straightforward:
@parameters t a β γ
@variables S(…) I(…) R(…)
@derivatives Dt’~t
@derivatives Da’~a
However, how does one describe integrals over one dimension (in this case, a)?
4 posts - 2 participants
