The model is implemented via a set of ordinary differential equations, defined by Equations (1)–(6).