“…After compilation, the program is executed. The evolution of the state variables of the model is obtained by numerical integration of differential equations described in the previous sections [15]. There are several methods for the integration of differential systems (Runge-Kutta, Adams, etc.).…”