ABSTRACT
We want to determine the solution X(t) of an ordinary differential equation (ODE) model given by a general function f t p X t( , , ( )), where p denotes the parameters,
d
d
X t
t f t p X t
( ) ( , , ( ))= (4.1)
and with a known initial condition, X0. The generality of this function emphasizes that the rate may depend on the time itself, a set of parameters, and the dependent variable. There are several numerical methods to solve ODEs such as Euler and Runge-Kutta (Swartzman and Kaluzny, 1987).