Linear ODE systems with SymPy's dsolve: the Bateman function of an oral dose

Afterwards you can solve linear ODE systems with initial values using SymPy's dsolve, handle equal rates with limit, and sum repeated doses in closed form.