The code was written in C++ because of performance considerations of the ODE solvers.