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.

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

Diele Fasma;Marangi Carmela
2010

Abstract

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.
2010
Istituto Applicazioni del Calcolo ''Mauro Picone''
Optimal growth models
Steady-state invariance
Partitioned symplectic Runge-Kutta methods
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.14243/116591
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 9
  • ???jsp.display-item.citation.isi??? 10
social impact