<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">TEL</journal-id><journal-title-group><journal-title>Theoretical Economics Letters</journal-title></journal-title-group><issn pub-type="epub">2162-2078</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/tel.2012.21002</article-id><article-id pub-id-type="publisher-id">TEL-17348</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Business&amp;Economics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Time Representation in Economics
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>tefano</surname><given-names>Bosi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Lionel</surname><given-names>Ragot</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>University of Evry, Evry, France</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>stefano.bosi@univ-evry.fr(TB)</email>;<email>lionel.ragot@univ-evry.fr(LR)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>23</day><month>02</month><year>2012</year></pub-date><volume>02</volume><issue>01</issue><fpage>10</fpage><lpage>15</lpage><history><date date-type="received"><day>October</day>	<month>31,</month>	<year>2011</year></date><date date-type="rev-recd"><day>November</day>	<month>27,</month>	<year>2011</year>	</date><date date-type="accepted"><day>December</day>	<month>5,</month>	<year>2011</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  In this paper, we study general polynomial discretizations in backward and forward looking, and the preservation of stability properties. We apply these results to the Ramsey model [4]. Its discrete-time version is a hybrid discretizations of a backward-looking budget constraint and a forward-looking Euler equation. Saddle-path stability is a robust property under discretization.
 
</p></abstract><kwd-group><kwd>Discretization; Growth</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Discretizations</title><p>Continuous-time systems can be approximated by discretetime systems. In the spirit of Krivine, Lesne and Treiner [<xref ref-type="bibr" rid="scirp.17348-ref2">2</xref>], we bridge continuous and discrete-time dynamics through general polynomial discretizations.</p><p>Discretizations can differ according to the step, the order and the direction of discretization. The step gives the length of the period in discrete time. The order is that of the Taylor expansion of a continuous-time model. The direction depends on the backward or forward-looking nature of this Taylor expansion. A hybrid discretization mixes backward and forward-looking approximations.</p><p>We want to show that the steady state is invariant to the step, the order and the direction of discretization and its continuous-time stability properties (sink, saddle, source) are preserved under a sufficiently small discretization step in any case (backward, forward or hybrid).</p><p>Instead of considering a continuous variable t and the corresponding position <img src="2-1500065\f6297493-702b-4e46-a53f-3f19095732d9.jpg" /> determined by an m-dimensional system of ordinary differential equations: <img src="2-1500065\e057f299-1ccf-4e04-a0e6-370f9c27086d.jpg" /> where <img src="2-1500065\5196783f-7e06-4d2b-b8e9-2f4542f29f74.jpg" /> jointly with the initial condition <img src="2-1500065\50d843fb-e79b-4d78-a73e-c81e39d506b3.jpg" /> let us pick up a regular sequence of time values:</p><p><img src="2-1500065\9f81c5a4-5930-4be8-b8f7-3153f503c3e9.jpg" />where h is a (possibly small) positive constant (discretization step), and the associated values:</p><p><img src="2-1500065\41e9f171-cfc0-4ae8-ad7e-db1159786dfb.jpg" />.</p><p>The path from <img src="2-1500065\79e2acb0-e530-45da-93b3-4d3b5d3547ab.jpg" /> to <img src="2-1500065\efd7457d-9814-4be1-82f8-15ca668673d2.jpg" /> can be reconstructed component by component through an appropriate integration of<img src="2-1500065\749a20ea-d5c9-46dd-9ee6-25260445c92f.jpg" />. Focusing on the ith component of the vector<img src="2-1500065\4acdb276-6c44-49be-bfd2-caa026133342.jpg" />, we can integrate the time derivative on the right or on the left to obtain, respectively,</p><p><img src="2-1500065\b3aef33e-c60e-4ef8-9587-611345a068cc.jpg" /></p><p>with<img src="2-1500065\7e124402-9356-490c-acce-8c054fb627f3.jpg" />. Defining</p><p><img src="2-1500065\15e7b7d4-1a00-40c4-adda-4d6de9e9b5d6.jpg" /></p><p>we get<img src="2-1500065\04716d6e-8389-4795-9391-d607709f74a0.jpg" />.</p><p>A discretization is an approximation of <img src="2-1500065\ac20ecfb-1b32-4298-914b-2fa0696f10f3.jpg" /> (<img src="2-1500065\e39f861c-ac87-43c1-9750-51b17da200f0.jpg" />) through a simpler function evaluated at <img src="2-1500065\7b9201a5-60dc-42ab-8789-4e16f9252c36.jpg" /> (<img src="2-1500065\92419e4c-ff09-4c42-8d5a-ff8eaeb1b503.jpg" />). The Euler-Taylor discretization is a polynomial approximation. Assuming that <img src="2-1500065\ae2da189-e9f1-4274-b50b-7bf4fed2c779.jpg" /> and considering the qth order polynomial, we obtain a backward or a forward-looking discretization:</p><disp-formula id="scirp.17348-formula51934"><label>(1)</label><graphic position="anchor" xlink:href="2-1500065\6efde27c-df1f-42a5-a223-2fe2710cee40.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.17348-formula51935"><label>(2)</label><graphic position="anchor" xlink:href="2-1500065\89365aa4-bb61-4763-ba26-859396b35c6f.jpg"  xlink:type="simple"/></disp-formula><p>because<img src="2-1500065\d740db86-bdb8-4231-8046-cdcd64198d75.jpg" />. A discretization is said to be hybrid if (1) holds for some components of the vector x and (2) holds for the others.</p><p>Setting<img src="2-1500065\d12cb70d-eff6-48cb-a21a-ec6961bdaa42.jpg" />, we obtain from (1) and (2) a first-order discretization:</p><p><img src="2-1500065\ff15c751-4e2f-435f-aa29-739363bbdd37.jpg" /></p><p>that is</p><disp-formula id="scirp.17348-formula51936"><label>(3)</label><graphic position="anchor" xlink:href="2-1500065\56e91a89-06d2-4e1f-934d-7f10e4705827.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.17348-formula51937"><label>(4)</label><graphic position="anchor" xlink:href="2-1500065\6c59a8a5-ce9e-4ff0-8087-69c5ce1da35a.jpg"  xlink:type="simple"/></disp-formula><p>where the subscript i denotes the ith component of the vector.</p><p>Equation (3) (respectively, (4)) constitutes a backwardlooking (forward-looking) discretization, because the variation <img src="2-1500065\2ce1b384-906a-4715-b908-e6a3d3ea3176.jpg" /> depends on the past value <img src="2-1500065\6bcb2aab-0725-4e58-b72f-e1943ca9526c.jpg" /> (future value<img src="2-1500065\8f5f0512-bbc1-47d3-b723-0e81df24074f.jpg" />) on the right-hand side. Equation (3) is the classical Euler discretization. In economics, forward-looking discretizations are of interest because agents behave according to their expectations.</p><p>The sequences <img src="2-1500065\9423a48a-7e42-4da4-b8ab-13585b819dcf.jpg" /> are approximations of the true sequence<img src="2-1500065\779c5848-bd9f-40c4-b82b-80adaaa35764.jpg" />, exact solution to system<img src="2-1500065\dc44cffd-1454-4704-8c43-b8fdaadd0cf7.jpg" />: the smaller h, the more accurate the representation.</p><p>Higher-order discretizations are also possible. Let us discretize the continuous-time dynamical system <img src="2-1500065\a2ad7260-22ad-4947-b5fa-7746400f7ffd.jpg" /> with <img src="2-1500065\dbb2c4d3-4188-4ae9-8048-36f1b23fb6b8.jpg" /> by second-order Taylor polynomials, that is approximate the ith component of <img src="2-1500065\f7c3906f-67fd-4c45-af20-6e3117fe8a54.jpg" /> with a quadratic form. Using (1) and (2), we obtain in backward and forward-looking, respectively:</p><p><img src="2-1500065\9434b000-ff7e-4b99-a686-b4bfadc4deec.jpg" /></p><p>where the subscript i denotes the ith component of the vector.</p><p>If f is an analytic function, infinite-order backward or forward discretizations converge exactly to <img src="2-1500065\e68f4880-976e-4fcd-af1b-b97a62c9287e.jpg" /> and (1) and (2) now hold with equality:</p><p><img src="2-1500065\04a44e63-ff27-449d-91af-5412c9114bc4.jpg" /></p><p>In this case, the Taylor polynomials become a convergent series and the discretized dynamics represent exactly the continuous-time system whatever the step h.</p><p>In general, a discretization is a closer approximation of a continuous-time system when the step h is smaller or the order of discretization q higher. The dynamic properties of a continuous-time system can be preserved lowering h or increasing q.</p></sec><sec id="s2"><title>2. Dynamic Equivalence</title><p>To compare continuous-time and discrete-time system, we study approximations in a neighborhood of the steady state and focus on the persistence of local stability properties.</p><p>Focus first on the steady state. The system <img src="2-1500065\29cd3542-d92e-4fef-9925-a73eb0938e43.jpg" /> and its discrete-time approximation <img src="2-1500065\5c5c34df-d1f3-4586-a7f6-241ace1de0dd.jpg" /> have the same steady state. Indeed, in both the cases, we require <img src="2-1500065\0c473f16-70b5-4c86-8eac-ea3c56e12c8c.jpg" /> (respectively, <img src="2-1500065\ac5b2b6f-7bb5-4cf1-adac-c7150734e085.jpg" />and<img src="2-1500065\980bb2a1-793c-4082-8d3f-55980f206c18.jpg" />). We further notice that the system of m equations <img src="2-1500065\ba9409bd-97ad-4159-ac06-99800cec6027.jpg" /> neither depends on the discretization degree h nor on the discretization method (forward or backward-looking). Therefore, the steady state is invariant to discretization.</p><p>Focus now on the stability properties. Are they preserved under discretization in a neighborhood of the steady state?</p><p>Without loss of generality, we consider two-dimensional dynamics. In the spirit of Samuelson [<xref ref-type="bibr" rid="scirp.17348-ref3">3</xref>], we can represent the stability properties in the plane of trace T and determinant D of the Jacobian matrix J of the system evaluated at the steady state.</p><p>In the following, the subscripts <img src="2-1500065\84dae581-53a2-4539-b684-4a0c24da77e3.jpg" /> and 1 will denote variables in continuous or discrete time respectively.</p><p>1) In continuous time, stability depends on the real part of these eigenvalues. If both the real parts are negative (positive), the steady state is a sink (source) (in this case, the trace of <img src="2-1500065\788346f3-4cd1-4798-9996-b3783b313314.jpg" /> is negative (positive) and the determinant of <img src="2-1500065\eebbe52f-a48f-40df-81b0-3a9427e75ac4.jpg" /> is positive (positive)). If the signs of the real parts are different, the eigenvalues are real and the steady state is a saddle point (in this case, the determinant is negative).</p><p>2) In discrete time, the modulus of an eigenvalue <img src="2-1500065\6d6c24e2-7ad7-4226-8c3a-7b2e5828de86.jpg" /> matters. When <img src="2-1500065\dfd59cff-87e8-479e-95ca-19d168a737a4.jpg" /> (<img src="2-1500065\2041e11c-557b-4101-b4b3-f496aed2783a.jpg" />) the eigenvalue is inside (outside) the unit circle. If both the eigenvalues are inside (outside) the unit circle, the steady state is a sink (source). If one is inside and the other outside the unit circle, the steady state is a saddle point.</p><p>We can evaluate the characteristic polynomial</p><p><img src="2-1500065\e874e1f7-41d2-4c44-bd36-8dc50ee17e68.jpg" /></p><p>at –1 and 1. Focus on the <img src="2-1500065\cabbb0c6-9751-4df6-ba5c-af1f4e83e894.jpg" />-plane. Along the line<img src="2-1500065\ba20b1c7-38f4-4ec3-a705-2c5e0f1f8e67.jpg" />, one eigenvalue is equal to 1 because <img src="2-1500065\b463563e-3b92-4393-84f3-dc6376698d0a.jpg" /> Along the line<img src="2-1500065\4b6899d9-da5f-46e6-95a0-a73e36638afe.jpg" />, one eigenvalue is equal to –1 because <img src="2-1500065\b14eb457-4227-4601-b366-07587eed7563.jpg" /> &#160;On the segment defined by <img src="2-1500065\c4847c8a-bd05-4c88-bb9d-bc0da644c50e.jpg" /> and<img src="2-1500065\aa362e15-bc64-4c14-ae58-b14b130d9515.jpg" />, the two eigenvalues are nonreal and conjugate with unit modulus. Consider first the points that neither belong to these lines nor to the segment. Inside the triangle defined by <img src="2-1500065\1abc26d0-8d03-463d-a6f3-2c1c8216a427.jpg" /> and <img src="2-1500065\bc434744-7c61-497f-95a8-579dac920fa1.jpg" /> the steady state is a sink. It is a saddle point if <img src="2-1500065\5b1a9795-dd7f-4ea7-92a9-602e40202945.jpg" /> lies on the left sides of both the lines <img src="2-1500065\fe764164-ba92-45ee-8fb3-5dbe03c5f7ae.jpg" /> and<img src="2-1500065\fd937c32-3c5f-4fdc-8430-8ceb67d024d5.jpg" />, or on the right sides of both of these lines (<img src="2-1500065\fe291750-644f-4eb3-9c5e-1fffc53750cb.jpg" />). It is a source otherwise.</p><p>At least, a two-dimensional system is required to study the three cases (sink, saddle and source) together and to consider hybrid discretizations. Without loss of generality, we linearize the following system of ordinary differential equations</p><disp-formula id="scirp.17348-formula51938"><label>(5)</label><graphic position="anchor" xlink:href="2-1500065\09338a07-97f6-42a3-877e-c5483b37ba45.jpg"  xlink:type="simple"/></disp-formula><p>Local dynamics around the steady state are represented by the Jacobian matrix <img src="2-1500065\27cb338e-a1db-4d96-9d64-0714635aba0c.jpg" /> evaluated at the steady state (<img src="2-1500065\c64c6648-1e28-4aae-a2db-0daa47f982ab.jpg" />).</p><p>We focus on first-order discretizations, but our equivalence results hold also for higher-order discretizations (see Bosi and Ragot [<xref ref-type="bibr" rid="scirp.17348-ref4">4</xref>]).</p><sec id="s2_1"><title>2.1. Backward-Looking Discretizations</title><p>We linearize the backward-looking discretization</p><disp-formula id="scirp.17348-formula51939"><label>(6)</label><graphic position="anchor" xlink:href="2-1500065\bc058daf-65ac-4fda-ba52-cd05a418c7c7.jpg"  xlink:type="simple"/></disp-formula><p>of the system (5) around the common steady state <img src="2-1500065\a7274670-cfa1-4c98-a507-4208be9e621b.jpg" /> and we obtain<img src="2-1500065\8408db4d-0bc9-45aa-b825-8491096f22a7.jpg" />, where I and <img src="2-1500065\a419074f-4f7a-43e5-ae54-443962f1c391.jpg" /> are the two-dimensional identity matrix and Jacobian matrix of system (6). We observe that <img src="2-1500065\7e85fcbd-fd93-489d-9f23-fb0fd4d52281.jpg" /> depends on the steady state x which, in turn, does not depend on h. Then, <img src="2-1500065\7ab110d5-54fe-4013-b9da-eedc45704d9e.jpg" />depends only linearly on h.</p><p>As above, let us denote the trace and determinant of <img src="2-1500065\501355f3-8271-4cb1-ad1c-c2bc53c8d0d7.jpg" /> and <img src="2-1500065\de253286-34fd-4dde-8bd0-6af24a641f11.jpg" /> by <img src="2-1500065\ec529020-8c1d-4346-a7cd-96746435224d.jpg" /> and <img src="2-1500065\f741c0e9-2ad6-4995-afbd-0c564650b48b.jpg" /> respectively. The characteristic polynomial in discrete time is given by<img src="2-1500065\e45b0d5c-ca1a-44ae-8a86-9f673d41e722.jpg" />, where</p><disp-formula id="scirp.17348-formula51940"><label>(7)</label><graphic position="anchor" xlink:href="2-1500065\c9b2ae03-79c2-45a6-80fe-f1d92acd4ce5.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.17348-formula51941"><label>(8)</label><graphic position="anchor" xlink:href="2-1500065\d8b7812f-d2a7-46ca-8d14-e7118ef9cc37.jpg"  xlink:type="simple"/></disp-formula><p>There are three critical values of the discretization step that determine the intervals of equivalence between the continuous and the discrete-time dynamics:</p><p><img src="2-1500065\2b0451f3-b937-4221-9828-041b9d2f36f8.jpg" /></p><p>Proposition 1 Consider<img src="2-1500065\dc246d98-ef04-45d0-b417-9ff3f711dc80.jpg" />.</p><p>1) Let the steady state be a sink in continuous time (<xref ref-type="fig" rid="fig1">Figure 1</xref>).</p><p>1.1) If<img src="2-1500065\2879ccfe-ee51-466b-83ba-6931693bb404.jpg" />, then the steady state is a sink in discrete time if <img src="2-1500065\cfc5529d-c9f4-4897-82c3-1775841c7ada.jpg" /> and a source if<img src="2-1500065\c725a241-07fd-4ddf-8d57-70f5eb4f4088.jpg" />.</p><p>1.2) If <img src="2-1500065\794e95f0-848d-4e98-9c52-5e8ee2f11d26.jpg" />then the steady state is a sink if <img src="2-1500065\42885dfd-6f0d-4e2a-877d-8bd2c4d6ca37.jpg" /> <img src="2-1500065\b07e638a-4df5-427b-b809-a3667f3d0ed2.jpg" />, a saddle if<img src="2-1500065\d9253e9e-2f58-460e-b868-dcf9e9eaf13d.jpg" />and source if<img src="2-1500065\26ae8054-d2bb-4edb-9fb1-eea941801069.jpg" />.</p><p>2) If the steady state is a saddle in continuous time, then the steady state is a saddle in discrete time if <img src="2-1500065\d0395d3a-2aa1-4c03-98cf-80685f9833d8.jpg" /> and source if <img src="2-1500065\4c9a36b4-93c3-4b40-9b4b-17fdca21c550.jpg" /> (<xref ref-type="fig" rid="fig2">Figure 2</xref>).</p><p>3) If the steady state is a source in continuous time, then the source property is preserved whatever <img src="2-1500065\2e34cda9-31f2-4f5a-905a-691c9b16ce2b.jpg" /> (<xref ref-type="fig" rid="fig3">Figure 3</xref>).</p><p>The system generically undergoes a Hopf bifurcation at <img src="2-1500065\01f83424-af55-4d86-a908-075a4b301d34.jpg" /> and flip bifurcations at<img src="2-1500065\4b6ef0ed-b90a-4f23-ac30-e7dad893a456.jpg" />,<img src="2-1500065\82785c84-be86-4387-bc68-dcd131d395cf.jpg" />.</p><p>Proof From (7) and (8), it is possible to plot a curve <img src="2-1500065\a7ed8171-747e-4e93-a400-14d0b5e98b5e.jpg" /> for each one of these different cases: <img src="2-1500065\b7f71a49-0d7e-4b8a-9678-fa98d58fdc1e.jpg" />given<img src="2-1500065\a4cba7e2-692b-4002-9c03-77d9ad707487.jpg" />.</p><p>1) Assume that the steady state is a sink in continuous time:<img src="2-1500065\f39b0a58-732c-4d08-a053-5ab3ecd986b6.jpg" />. According to (8),<img src="2-1500065\8fc2ae5f-bda0-4815-85f3-6db3dce6ce23.jpg" />. Focus on two cases: (1.1) <img src="2-1500065\922fe5f4-f826-4c36-b9cd-74e574b90847.jpg" />and (1.2)<img src="2-1500065\2e29dd2d-0f35-44cc-9a8c-c21edda011bc.jpg" />.</p><p>1.1) If <img src="2-1500065\e0bce767-67b9-411a-89d6-ffcc1b42e0dd.jpg" /> then always <img src="2-1500065\33b822df-0ff9-458c-a988-4c9d275512fe.jpg" /> that is <img src="2-1500065\7de81e62-af7d-497b-8324-9dd5fc7c1cb6.jpg" /> So, the steady state is a sink if<img src="2-1500065\18e41277-d69f-41d9-9406-b67591a4c046.jpg" />, that is if<img src="2-1500065\9df4689c-f992-46d3-8bf1-4cbcd4861368.jpg" />, and a source if<img src="2-1500065\62da1020-1d9d-4a6c-83a4-0db165c2a024.jpg" />. This case corresponds to the upper parabola in <xref ref-type="fig" rid="fig1">Figure 1</xref>. Increasing h away from zero means moving away from the point where<img src="2-1500065\1ff01ac4-d1c2-49e7-ba8f-15abe97f5055.jpg" />, along the parabola.</p><p>1.2) If <img src="2-1500065\273a8240-cbb0-4eef-a4b6-3571ae3ceb30.jpg" /> then <img src="2-1500065\c6582ce8-6012-4821-b0cd-f8e159e65176.jpg" /> if and only if <img src="2-1500065\72bceacd-8121-4baa-9cf9-d4b5dbc4d688.jpg" /> <img src="2-1500065\43ada512-5431-4015-95e9-6c4df3ed89a2.jpg" />. In addition, <img src="2-1500065\073c27e1-0273-46dd-8695-f70477637464.jpg" />if and only if<img src="2-1500065\ac18d2e3-63b9-43a8-9e53-22a4abe57ad1.jpg" />. We notice also that<img src="2-1500065\2651c8fc-540c-4a5d-a4cd-a02ace7c6b79.jpg" />. Then, the steady state is a sink if<img src="2-1500065\8ecffb28-f421-4973-9c17-8e7a2ab805e0.jpg" />, a saddle if <img src="2-1500065\a226df5f-31fa-4ed1-98ea-3b4aba901162.jpg" /> and a source if<img src="2-1500065\889909b4-3e86-44ad-9077-61a6e0de4c3f.jpg" />. This case corresponds to the lower parabola in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p>2) Assume now that the steady state is a saddle in continuous time: <img src="2-1500065\f44a17ae-75c3-48c4-8c10-2b8efb9afe1d.jpg" />According to (8), <img src="2-1500065\d69deca6-eefe-4b34-8bfe-4aea6a28df24.jpg" />We observe that <img src="2-1500065\5585e2cf-ba42-48f5-93e2-e9954e85f638.jpg" /> and that <img src="2-1500065\60751c52-e5b1-483c-ba1a-f5ee2302fcc8.jpg" /> if and only if <img src="2-1500065\1093e487-3782-43b3-a9d6-fe5d5c672041.jpg" /> Thus, the steady state is a saddle if <img src="2-1500065\3e81825f-0e57-46a3-aa3f-91d0f30ea0d0.jpg" /> and a source if <img src="2-1500065\b7351b08-22c4-4175-a063-86931c7d57fd.jpg" /> If <img src="2-1500065\112ce58d-f4c7-4f61-b325-d579fc9aa7c6.jpg" /> (<img src="2-1500065\3eb6393e-8ef8-4d41-b26d-a97244ac7d0d.jpg" />)the curve <img src="2-1500065\bb2e339d-8d14-4a5a-9591-5bef676fc9af.jpg" /> is represented by the leftward (rightward) branch of parabola in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><p>3) Assume now that the steady state is a source in continuous time: <img src="2-1500065\d5494af3-2faa-40ce-b3b7-c3eb3cce461c.jpg" />and<img src="2-1500065\f7da4694-dae6-45f1-953d-cd48fff79e19.jpg" />. (7) and (8) imply <img src="2-1500065\fceab324-e821-4388-8019-f9917dba17fc.jpg" /> and <img src="2-1500065\45bad8ac-ec89-407a-a509-72ea60c235f9.jpg" /> for every<img src="2-1500065\5fc93748-3a48-4f2f-8d14-8d42109b5620.jpg" />. Therefore the source property is preserved whatever<img src="2-1500065\2bc5ab0b-c04d-4c20-8cd0-5d963bea070b.jpg" />. The branch of parabola in <xref ref-type="fig" rid="fig3">Figure 3</xref> represents this case.</p><p>Corollary 2 (topological equivalence in backward looking) In any case of Proposition 1, there exists a nonempty interval <img src="2-1500065\d1a7ca0b-2e13-4e15-b258-f78e8b823eb4.jpg" /> for the discretization step h where the stability properties of the continuous-time system are preserved.</p><p>Proof Straightforward. Simply observe that, in the case (3),<img src="2-1500065\644af0a3-a8bd-41a4-97ea-2333d28b033b.jpg" />.</p></sec><sec id="s2_2"><title>2.2. Forward-Looking Discretizations</title><p>We linearize now the forward-looking discretization</p><disp-formula id="scirp.17348-formula51942"><label>(9)</label><graphic position="anchor" xlink:href="2-1500065\5b7acaac-da8f-4dbf-8db5-59da8964fb23.jpg"  xlink:type="simple"/></disp-formula><p>of system (5) around the common steady state <img src="2-1500065\2a6da51a-f968-4a2c-baeb-baba055d60ec.jpg" /> to obtain</p><p><img src="2-1500065\1cfe509f-fd0c-4434-857c-c30d575e0f04.jpg" /></p><p>Differently from the previous case, the Jacobian matrix of system (9) <img src="2-1500065\4a98f8b7-ac15-4c9d-b915-861ac6a5b271.jpg" />is no longer linear in h. The trace and the determinant of <img src="2-1500065\a25fe640-23e1-445c-978a-88638a6ee1df.jpg" /> are now given by</p><p><img src="2-1500065\30c0a3c8-d6bc-4d63-a48d-046107333a24.jpg" /></p><p><img src="2-1500065\89358292-73d7-401a-a72e-6a9f932e2cfd.jpg" /></p><p>As above, we set three critical values:<img src="2-1500065\47a0d85e-3074-41a2-87ac-389fac197a31.jpg" />,</p><p><img src="2-1500065\2a548cc7-268c-4553-a9a4-8e7d8e6ac1ea.jpg" /></p><p>Proposition 3 Consider<img src="2-1500065\80bc9a5b-9b51-43f2-b74c-a1ba30d408c9.jpg" />.</p><p>1) If the steady state is a sink in continuous time, then the sink property is preserved in discrete time whatever<img src="2-1500065\b5814f1e-a6cc-4ce5-b3eb-fe4d9417836b.jpg" />.</p><p>2) Let the steady state be a saddle in continuous time.</p><p>2.1) If<img src="2-1500065\7a4b976e-f08b-4b03-a6bc-072395804b63.jpg" />, then the steady state is a saddle.</p><p>2.2) If<img src="2-1500065\63cc8e02-6919-4ac1-b031-ff62063abb0b.jpg" />, then the steady state is a saddle if <img src="2-1500065\d784b1c1-9a25-44b4-8762-d36708366bcc.jpg" /> and a sink if<img src="2-1500065\d59a3174-c0ba-4fa7-b75d-a722323e10f8.jpg" />.</p><p>3) Let the steady state be a source in continuous time.</p><p>3.1) Let<img src="2-1500065\f2ddcae5-cbc3-4d9e-9812-a5984fd39b0f.jpg" />. If<img src="2-1500065\425be5e5-d77c-4814-a604-7f47c011f700.jpg" />, then the source property is preserved whatever <img src="2-1500065\57bfa040-6f31-4b20-a676-3851b9d82141.jpg" /> If<img src="2-1500065\4d18f480-c1e9-4454-a771-cce5fb8194dc.jpg" />, then the steady state is a source if <img src="2-1500065\20d76c3c-0a2c-4d62-b717-17499744e933.jpg" /> or<img src="2-1500065\fe03bf14-a0af-447e-92ae-85f38920263b.jpg" />, and a saddle if<img src="2-1500065\8087a03e-1d80-494d-a237-f23ed9dc91a6.jpg" />.</p><p>3.2) Let<img src="2-1500065\4e96b6fb-ae40-4d6f-b336-a8e86b9fea05.jpg" />. If<img src="2-1500065\6ad0c200-8f99-41b7-9ba8-50079c65d1a1.jpg" />, then the steady state is a source if <img src="2-1500065\45950bcb-324c-4357-abbf-ac998f5caa6d.jpg" /> and a sink if<img src="2-1500065\f4625b97-06a7-4d38-ac55-1b38e1d497d2.jpg" />. If<img src="2-1500065\ba5ff5c4-30cf-4044-8214-9a64384c130f.jpg" />, then the steady state is a source if <img src="2-1500065\43e4b5bd-b92a-4e1f-93ae-dedfe5ebbdcf.jpg" /> a saddle if <img src="2-1500065\282e3dea-f754-49ff-89d1-004d07f47883.jpg" /> and a sink if<img src="2-1500065\e8750e52-6eb4-42c3-9a17-a5a80ba0a4b4.jpg" />.</p><p>The system generically undergoes a Hopf bifurcation at <img src="2-1500065\02537be1-1465-44c0-a2ec-0fb201721cec.jpg" /> and flip bifurcations at<img src="2-1500065\861d2b9f-3a4c-47fc-9a73-37ea68146bc4.jpg" />,<img src="2-1500065\49be2bbb-a83b-449c-9f65-32f5afcb4e62.jpg" />.</p><p>Proof The proof is similar to that of Proposition 1. See Bosi and Ragot [<xref ref-type="bibr" rid="scirp.17348-ref4">4</xref>] for more details.</p><p>Corollary 4 (topological equivalence in forward looking) In every case of Proposition 3, there exists a nonempty interval <img src="2-1500065\0f697403-414f-4ab3-bceb-7b9b649f6180.jpg" /> for the discretization step h where the stability properties of the continuous-time system are preserved.</p><p>Proof Straightforward. Simply observe that, in cases (1) and (2.1),<img src="2-1500065\1496e2c0-2381-4414-984e-621e25dd74b0.jpg" />. The same happens in the case (3.1) if<img src="2-1500065\e8b2a35d-2969-410e-9dfd-f05da0a0f438.jpg" />.</p></sec><sec id="s2_3"><title>2.3. Hybrid Discretizations</title><p>In economics, many higher-dimensional models require a hybrid discretization to recover the equivalence between discrete and continuous time, that is a mix of discretization in backward and forward looking. Without loss of generality, we consider a system where the first equation is discretized backward and the second one forward. Thus, the system of differential equations (5) becomes:</p><disp-formula id="scirp.17348-formula51943"><label>(10)</label><graphic position="anchor" xlink:href="2-1500065\5e907a11-26ce-45d8-a11c-d9c2dbbd2469.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.17348-formula51944"><label>(11)</label><graphic position="anchor" xlink:href="2-1500065\379f1fbc-bf2e-47ee-9ee3-dddddf4fefc4.jpg"  xlink:type="simple"/></disp-formula><p>The steady state is invariant to the choice of time and to the type of discretization (backward/forward). The trace and the determinant of the Jacobian matrix <img src="2-1500065\ec3d4964-350a-48c1-9d65-0ede87db4c4e.jpg" /> of the hybrid system (10)-(11) become</p><disp-formula id="scirp.17348-formula51945"><label>(12)</label><graphic position="anchor" xlink:href="2-1500065\f3ac600b-1f10-424d-92a7-82f295f62f7d.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.17348-formula51946"><label>(13)</label><graphic position="anchor" xlink:href="2-1500065\f7bed24a-ca43-4086-a946-53a9a973a3aa.jpg"  xlink:type="simple"/></disp-formula><p>Notice that, in the particular case<img src="2-1500065\ba5813f5-46c0-4a8a-9ae4-559f72e7a158.jpg" />, (12) and (13) write</p><p><img src="2-1500065\5649681b-4c68-479e-a47d-4657dca30e32.jpg" /></p><p>Let&#160;&#160; <img src="2-1500065\f1fb351e-480b-492d-b87b-984cd82e2bc8.jpg" /></p><p><img src="2-1500065\2a8c9178-3c5a-4f9f-b1e1-12f4a2b6ebf1.jpg" /></p><p>where<img src="2-1500065\cf8b83fe-8427-4c34-afa0-ff0828b9ac08.jpg" />.</p><p>Proposition 5 Consider<img src="2-1500065\0976b233-70b8-4658-9e9d-502bcddc839f.jpg" />.</p><p>1) Let<img src="2-1500065\9526aeec-357e-459a-985c-de9a53f653f7.jpg" />.</p><p>1.1) If the steady state is a sink in continuous time, then the steady state in discrete time is a sink if <img src="2-1500065\24b20d59-291c-457f-8680-24c04622ca39.jpg" />, and a saddle if<img src="2-1500065\aca600fc-342b-48af-ac17-cf361e42ec99.jpg" />.</p><p>1.2) Let the steady state be a saddle in continuous time.</p><p>1.2.1) If <img src="2-1500065\6d95472d-bfef-417d-9c48-f118b2f8cd28.jpg" /> or<img src="2-1500065\a7924ce2-6286-4115-9899-1fe5b6bb9117.jpg" />, then the steady state is a saddle point.</p><p>1.2.2) If <img src="2-1500065\a4d3615c-7653-4186-88a7-d2cc3bc59744.jpg" /> and<img src="2-1500065\a5cbaaf8-eb48-4fe8-bf0a-18e37a0ec468.jpg" />, then the steady state is a saddle if <img src="2-1500065\2cd6d94a-d08f-45d1-b980-c15b15dad8a9.jpg" /> or<img src="2-1500065\f805688d-c2d8-4810-891a-8cb999a4dc93.jpg" />, and a source if<img src="2-1500065\83d086a1-76b3-41fd-b03b-b79cac7a009d.jpg" />.</p><p>1.3) If the steady state is a source in continuous time, then the steady state is a source if <img src="2-1500065\195a9d19-5f76-49f0-9607-913cc22c8d7c.jpg" /> and a saddle if<img src="2-1500065\7f504535-3e41-4e2c-80dd-9edd2a6bc048.jpg" />.</p><p>2) Let <img src="2-1500065\00c1c0ef-62c4-4ca4-8a73-b6c999e29662.jpg" /> with<img src="2-1500065\211c73bb-1750-42e8-985f-3e9f522bf7f3.jpg" />. All the previous cases hold, provided we restrict the analysis to the interval<img src="2-1500065\cb20f94d-0a04-4e23-b45d-df7bd2a602c6.jpg" />.</p><p>The system generically undergoes a Hopf bifurcation at <img src="2-1500065\547ef680-d16f-4361-8e69-3c5317cf3fe5.jpg" /> and a flip bifurcation at<img src="2-1500065\2a21b0f2-81c1-414b-98bf-5efaa9ae0935.jpg" />,<img src="2-1500065\423a09c1-f305-4f6c-93af-48109c1ef72c.jpg" />.</p><p>Proof The proof is similar to that of Proposition 1. See Bosi and Ragot [<xref ref-type="bibr" rid="scirp.17348-ref4">4</xref>] for more details.</p><p>Corollary 6 (topological equivalence in hybrid looking). In every case of Proposition 5, there exists a nonempty interval <img src="2-1500065\dbacfffe-9619-4ada-aa1f-fba41016bca6.jpg" /> for the discretization step h where the stability properties of the continuous-time system are preserved.</p><p>Proof Straightforward. Simply observe that, in the case (1.2.1),<img src="2-1500065\dc14031e-ba94-4046-bb81-2398eafb22ec.jpg" />.</p></sec></sec><sec id="s3"><title>3. Ramsey Model</title><p>In the seminal Ramsey [<xref ref-type="bibr" rid="scirp.17348-ref1">1</xref>], the planner maximizes the undiscounted dynastic utility:<img src="2-1500065\bc3416dc-e5c8-4397-bbbd-884409c6cb02.jpg" />, under a resource constraint <img src="2-1500065\085903a3-38ff-4a96-9070-5054bd4a2667.jpg" /> where <img src="2-1500065\0631224b-8d8b-479d-a360-15374db414fa.jpg" /> and <img src="2-1500065\191a7bca-a559-415f-b206-6b5ede7dfa6d.jpg" /> denote the individual capital and consumption. The initial endowment <img src="2-1500065\b01ecd4d-f727-4ed1-97bc-cbd491a275db.jpg" /> is given.</p><p>The intensive production function <img src="2-1500065\fd419c69-36b3-45c7-b514-905ef9db580d.jpg" /> is strictly increasing and strictly concave in the capital intensity and satisfies the Inada conditions. The felicity <img src="2-1500065\01b35ad1-004e-4dfe-b3bd-e0595ec91b99.jpg" /> is also strictly increasing and strictly concave in the consumption level. c denotes the bliss point, that is the steady state value of consumption: <img src="2-1500065\33621038-72c2-42a6-b997-22d0f96dc6f9.jpg" />with<img src="2-1500065\c270eab1-e5e8-4080-8aff-d79097763b66.jpg" />.</p><p>The planner maximizes the Hamiltonian:</p><p><img src="2-1500065\3bfe058d-f668-4299-99c2-595812b661ea.jpg" /></p><p>to find the first-order conditions:</p><disp-formula id="scirp.17348-formula51947"><label>(14)</label><graphic position="anchor" xlink:href="2-1500065\5b19e404-61b7-4d72-aa63-074dc17ce46d.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.17348-formula51948"><label>(15)</label><graphic position="anchor" xlink:href="2-1500065\99b51de1-5328-4e82-8450-8b9f0f746cff.jpg"  xlink:type="simple"/></disp-formula><p>where<img src="2-1500065\eaf997a4-af37-42c3-8609-aada8ebdd59b.jpg" />. The strict concavity of u ensures that <img src="2-1500065\697383c1-e78b-44e5-befd-9c9714cbdfcb.jpg" /> is a well-defined function of the multiplier<img src="2-1500065\adf273ec-a870-4893-ad24-edb6e40db0a3.jpg" />.</p><p>In discrete time, the planner maximizes <img src="2-1500065\889dc6f4-8b5e-4c05-9d5a-7a02e61445bb.jpg" /> under a sequence of resource constraints:<img src="2-1500065\ff27386e-efba-425c-9f5b-412c32093ec0.jpg" />, to obtain the firstorder conditions:</p><disp-formula id="scirp.17348-formula51949"><label>(16)</label><graphic position="anchor" xlink:href="2-1500065\f74b3ead-2c92-4453-a43e-a20c8ab4d4d0.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.17348-formula51950"><label>(17)</label><graphic position="anchor" xlink:href="2-1500065\747af3d1-a908-436d-9432-bdbc818579d0.jpg"  xlink:type="simple"/></disp-formula><p>We want to prove that the discrete-time system (16)- (17) is a discretization of the continuous-time system (14)- (15).</p><p>Proposition 7 The discrete-time Ramsey model comes from a first-order hybrid Euler discretization of the continuous-time model, that is a backward-looking discretization of the resource constraint (14) and a forwardlooking discretization of the Euler Equation (15), with a unit step.</p><p>Proof Under the backward-looking linear discretization of the continuous-time resource constraint (14):</p><disp-formula id="scirp.17348-formula51951"><label>(18)</label><graphic position="anchor" xlink:href="2-1500065\edc30108-1cd8-4644-83fd-03e855b08926.jpg"  xlink:type="simple"/></disp-formula><p>we recover exactly the discrete-time resource constraint (16) with a unit discretization step (<img src="2-1500065\708872cf-98f3-416a-b59d-061e063ff1fa.jpg" />). However, the intertemporal arbitrage requires a forward-looking discretization. Focus on (15) and apply (4):</p><p><img src="2-1500065\fe5a5d30-9c8a-49d9-a636-8e149984ad33.jpg" /></p><p>to obtain</p><disp-formula id="scirp.17348-formula51952"><label>(19)</label><graphic position="anchor" xlink:href="2-1500065\bdf7dc60-e9a2-4b82-8f51-48e4c31cfcb0.jpg"  xlink:type="simple"/></disp-formula><p>which gives exactly the discrete-time Euler Equation (17) under a unit discretization step<img src="2-1500065\d484eb15-270b-45f8-9b2f-721ac1965ece.jpg" />.</p><p>The forward-looking discretization of (15) is more suitable to capture saving decisions. Indeed, the expected productivity affects the arbitrage between consumption today and consumption tomorrow.</p><p>Let us consider the steady state. For all the three dynamical systems (14)-(15), (16)-(17) and (18)-(19) the steady state is defined by: <img src="2-1500065\872f0661-502a-4e73-9b59-5ddf4da0a271.jpg" />and <img src="2-1500065\b67fd381-e4a6-48a8-94f7-6af97a9a398a.jpg" /> (assumptions on technology and preferences ensure its existence and uniqueness).</p><p>Focus now on the stability properties. The Jacobian matrix <img src="2-1500065\c0182148-13f2-43ce-8a1d-da10bbc1b72f.jpg" /> of the continuous time system (14)-(15) is given by:</p><p><img src="2-1500065\7dd5d44a-e1b1-4f02-bd94-ba7d98cd5946.jpg" /></p><p>where<img src="2-1500065\0670b0c4-5f07-4841-be0d-c95e4e422863.jpg" />, <img src="2-1500065\f4c1170a-3f0c-4a87-9678-e28e922110b7.jpg" />with</p><p><img src="2-1500065\7fc1cd81-a89a-4b87-bb53-33002dfb81a8.jpg" />, <img src="2-1500065\5fa11229-75fb-4e27-b4e4-5ccf74ebf82a.jpg" /></p><p>and<img src="2-1500065\68c012e5-0a7f-4641-a582-7d47e7e1490b.jpg" />. The trace and the determinant in continuous time are <img src="2-1500065\b830680c-9eb2-4571-9eb5-254b47700896.jpg" /> and</p><p><img src="2-1500065\2e5983a9-a5b6-406f-86c4-edadaf0f1eca.jpg" />.</p><p>Notice that <img src="2-1500065\26ff442e-b47c-41a2-bb2b-c75f54a9e65f.jpg" /> implies the saddle-path stability property.</p><p>The hybrid Euler discretization (18)-(19) is consistent with the continuous-time case.</p><p>Proposition 8 The steady state of the discretized model is a saddle point (as in the continuous-time case) whatever the discretization step h.</p><p>Proof The Jacobian matrix <img src="2-1500065\e3da78e3-f282-435c-9402-ab61ce5e72b3.jpg" /> of the hybrid Euler discretization (18)-(19) is:</p><p><img src="2-1500065\f4decc2a-40ae-4c78-a137-b2d22e2abe7f.jpg" /></p><p>where A and B are defined above. The trace and determinant become <img src="2-1500065\ada9a75f-a5b4-461c-a86c-c7eb3bd0c718.jpg" /> and<img src="2-1500065\3646469e-5fbb-47ce-998c-cb166e730f09.jpg" />. We obtain <img src="2-1500065\5d7153ba-5dba-4bfb-86ec-972e70f2f4b0.jpg" /> and we recover the saddle-path stability property, whatever the discretization step h. There is no room for bifurcations, as in the continuous-time case.</p><p>Therefore, the saddle-path stability is a robust property of the Ramsey model because it holds whatever the discretization step.</p></sec><sec id="s4"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.17348-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">H. Krivine, A. Lesne and J. Treiner, “Discrete-Time and Con-tinuous-Time Modelling: Some Bridges and Gaps,” Mathe-matical Structures in Computer Science, Vol. 17, No. 2, 2007, pp. 261-276.  
doi:10.1017/S0960129507005981</mixed-citation></ref><ref id="scirp.17348-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">P. A. Samuelson, “The Stability of Equilibrium: Comparative Statics and Dynamics,” Econometrica, Vol. 9, No. 2, 1941, pp. 97-120. doi:10.2307/1906872</mixed-citation></ref><ref id="scirp.17348-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">S. Bosi and L. Ragot, “Time, Bifurcations and Economic Applications,” CES Working Papers, 2009.28, University of Paris 1, 2009. </mixed-citation></ref><ref id="scirp.17348-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">F. Ramsey, “A Mathematical Theory of Saving,” Economic Journal, Vol. 38, No. 152, 1928, pp. 543-559.  
doi:10.2307/2224098</mixed-citation></ref></ref-list></back></article>