Numerical DAE approach for solving a system dynamics problem