I’d like to know how/if the following integral equation can be solved in ModelingToolkit:
using ModelingToolkit using DomainSets @parameters t β γ @variables S(..) I(..) R(..) It = Integral(t in DomainSets.ClosedInterval(0,t)) eqs = [S(t) ~ S(0.) - It(β*S(t)*I(t)), I(t) ~ I(0.) + It(β*S(t)*I(t)) - It(γ*I(t)), R(t) ~ R(0.) + It(γ*I(t))];
I know it’s a trivial example that can be rewritten as ODEs, but knowing how to do this would open up a lot.