Steady-state invariance in high-order Runge-Kutta discretization of optimal growth models

This work deals with infinite horizon optimal growth models and uses the results in the Mercenier and Michel (1994a) paper as a starting point. Mercenier and Michel (1994a) provide a one-stage Runge-Kutta discretization of the above-mentioned models which preserves the steady state of the theoretical solution. They call this feature the "steady-state invariance property". We generalize the result of their study by considering discrete models arising from the adoption of s-stage Runge-Kutta schemes. We show that the steady-state invariance property requires two different Runge-Kutta schemes for approximating the state variables and the exponential term in the objective function. This kind of discretization is well-known in literature as a partitioned symplectic Runge-Kutta scheme. Its main consequence is that it is possible to rely on the well-stated theory of order for considering more accurate methods which generalize the first order Mercenier and Michel algorithm. Numerical examples show the efficiency and accuracy of the proposed methods up to the fourth order, when applied to test models.