Next Article in Journal
Building Damage and Business Continuity Management in the Event of Natural Hazards: Case Study of the 2004 Tsunami in Sri Lanka
Next Article in Special Issue
Mind Sized World Models
Previous Article in Journal / Special Issue
Growth Is the Problem; Equality Is the Solution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotic Behavior of a Delay Differential Neoclassical Growth Model

by
Akio Matsumoto
1,* and
Ferenc Szidarovszky
2
1
Department of Economics, Chuo University, 742-1, Higashi-Nakano, Hachioji, Tokyo, 192-0393, Japan
2
Department of Applied Mathematics, University of Pecs, Ifjusag u. 6, H-7624, Pecs, Hungary
*
Author to whom correspondence should be addressed.
Sustainability 2013, 5(2), 440-455; https://doi.org/10.3390/su5020440
Submission received: 27 November 2012 / Revised: 17 January 2013 / Accepted: 22 January 2013 / Published: 31 January 2013
(This article belongs to the Special Issue 40th Anniversary of 'The Limits to Growth')

Abstract

:
A neoclassical growth model is examined with a special mound-shaped production function. Continuous time scales are assumed and a complete steady state and stability analysis is presented. Fixed delay is then assumed and it is shown how the asymptotic stability of the steady state is lost if the delay reaches a certain threshold, where Hopf bifurcation occurs. In the case of continuously distriubuted delays, we show that with small average delays stability is preserved, then lost at a threshold, then it is regained if the average delay becomes sufficiently large. The occurence of Hopf bifurcation is shown at both critical values.

1. Introduction

The examination of economic growth models is one of the most frequently discussed issues in mathematical economics. Day [1,2] has investigated a neoclassical growth model, and a productivity and population growth model and showed the emergence of complex behavior even under simple economic structure. His models were based on discrete time scales and a mound-shaped production function that represented the negative effect of pollution resulting from increasing capital. It was demonstrated by numerical computations that these models could generate cyclic and even chaotic behavior. Following Day’s pioneering works, a lot of effort has been given to the understanding of complex economic dynamics. Day [3], Puu [4] and Bischi et al. [5] present the earlier contributions of this field. A large number of studies assumed discrete time scales. Li and Yorke [6] have introduced the “period-three condition” to detect chaos, which has many applications in first-order nonlinear difference equations. The papers collected by Rosser [7] offer many applications. Only a few studies are devoted to the case of continuous time scales, since there is no general criterion to detect chaos and the system must have at least three dimensions.
In this paper, we will examine an extension of the neoclassical growth model, which can be traced back to the early works of Solow [8] and Swan [9]. The neoclassical growth model is constructed with the two (usually implicitly mentioned) assumptions; one is the full-employment of labor and capital and the other is instantaneous adjustment in the output market. Thus, it is suitable for describing the long-run behavior of the economy. Due to the well-behaved production function, the steady state of the model is usually asymptotically stable. However, it is often observed in reality that growth path exhibits persistent fluctuations. The neoclassical model could be the good point of departure to show how such persistent behavior can emerge when nonlinearities and a production lag are present. Matsumoto and Szidarovszky [10] attempt to fill the gap and have introduced a neoclassical model with a mound-shaped production function that was assumed to be a Cobb-Douglas type function of the form F(x)=Axa(1 - x)b with x being the capital per unit labor. Although they show emergence of erratic fluctuations in the capital accumulation process, the production function is restrictive in the sense that it is defined only in the unit interval. This paper modifies this drawback, considers another type of mound-shaped production function and will examine the stability of the steady state with and without time delays in the continuous-time framework. Two kinds of delays will be discussed, fixed and continuously distributed (continuously hereafter) delays. We keep the relatively simple model of Matsumoto and Szidarovszky [10] in order to be able to compare the results and to illustrate that complex dynamics can be generated under simple economic assumptions with both function types.
This paper develops as follows. First the mathematical model is formulated without time delays, and complete steady state and stability analysis is presented. Then, models with fixed delays and then with continuous delays are introduced and complete stability analysis is given. The last section concludes the paper.

2. The Mathematical Model

Matsumoto and Szidarovszky [10] have introduced a special growth model of the form
Sustainability 05 00440 i001
where x is the capital per labor, s and α are positive parameters where s ∈ (0,1) is the average propensity to save and α = n+sμ with μ being the depreciation ratio of capital and n the growth rate of labor. In applying the Cobb-Douglas type function F(x) = Axα (1 - x)b the value of x has to be normalized into the (0,1) interval. In this paper we will assume that
F(x) = εxγe-δx,
which is a mound-shaped function in which x can take any positive value. This function has zero value at x = 0, converges to zero as x → ∞, increases for x < γ / δ and decreases for x > γ / δ, so it has its maximum at x = γ / δ. The drawbacks of the neoclassical production function that ignores natural resource or energy are partly remedied in our function. δ of the exponential term reflecting a strength of a “negative effect” caused by increasing concentration of capital. The value of δ is exogenously given, however, can be thought to be determined by a damaging degree of natural environment or energy resources. With the modified function, the mathematical model becomes
Sustainability 05 00440 i069
where α, γ, δ and β = s ε are positive parameters. The number of steady states and their locations depend on the specific values of the model parameters. We will consider three different cases: γ < 1, γ = 1 and γ > 1. We can give the following interpretation for the value of the parameter γ of the production function; γ can be thought as a proxy for measuring returns to scale of the production function. Indeed, when x is small (i.e., the exponential term is close to unity), output increases more than unity, exactly unity and less than unity if γ > 1, γ = 1 and γ < 1, respectively. Let now f(x) denote the right hand side of equation (1). If x(0) = 0, then the identically zero function is a solution which case is not interesting from the economic point of view and is eliminated from further considerations
Case I.
Assume first that γ < 1. The steady states are the solutions of f(x) = 0. Notice that f(0) = 0, so zero is a steady state. f(x) converges to - ∞ as x → ∞ Since
(x) = -α + βγxγ-1e-δx- βδxγe-δx
(x) converges to ∞ as x tends to zero with positive values. Hence f(x) increases for small values of x > 0. The steady state equation f(x) = 0 can be written as
xγ (-αx1-γ + βe-δx) = 0
so the positive steady state is the unique solution of equation
αx1-γ = βe-δx
The left hand side is zero at x = 0 and strictly increasing, furthermore, converges to ∞ as x tends to infinity. The right hand side is β > 0 at x = 0, strictly decreases and converges to zero as x → ∞. Hence there is a unique solution Sustainability 05 00440 i065 > 0 of (3), and f(x) > 0 if x < Sustainability 05 00440 i065 and f(x) < 0 as x > Sustainability 05 00440 i065 These relations imply that if x(0) < Sustainability 05 00440 i065, then x (t) increases and converges to Sustainability 05 00440 i065, and if x (0) > Sustainability 05 00440 i065, then x (t) decreases and converges to Sustainability 05 00440 i065. If x(0) = Sustainability 05 00440 i065, then x (t) remains Sustainability 05 00440 i065 for all x > 0. Thus Sustainability 05 00440 i065 is globally asymptotically stable.
Case II.
Assume next that γ = 1. Then the steady state equation has the form
x(-α + βe-δx) = 0
so zero is a steady state and there is a unique root of the second factor,
Sustainability 05 00440 i003
If βα, then the value of f(x) is negative for all x > 0. Therefore x(t) is decreasing and converges to zero with arbitrary x (0) > 0. If β > α, then Sustainability 05 00440 i065 > 0, furthermore f(x) > 0 as x < Sustainability 05 00440 i065, and f(x) < 0 as x > Sustainability 05 00440 i065. If x (0) < Sustainability 05 00440 i065, then x (t) increases and if x (0) > Sustainability 05 00440 i065, then x decreases and converges to Sustainability 05 00440 i065, and if x (0) = Sustainability 05 00440 i065, then x (t) remains Sustainability 05 00440 i065 for all T > 0 Hence Sustainability 05 00440 i065 is globally asymptotically stable.
Case III.
Consider finally the case of γ > 1. The steady state equation has now the form
x(-α + βxγ-1e-δx) = 0
so zero is a steady state again, and any other steady state is the solution of equation
g(x) = -α + βxγ-1e-δx = 0
Notice that
Sustainability 05 00440 i004
and
g´(x) = β(γ-1)xγ-2e-δx - βδxγ-1e-δx
= βxγ-2e-δx(γ-1 - δx)
Therefore g(x) as its global maximum at
Sustainability 05 00440 i005
increases for x < Sustainability 05 00440 i066 and decreases for x > Sustainability 05 00440 i066. Now we have three sub-cases.
(i) if g( Sustainability 05 00440 i066) < 0, then there is no positive steady state and with arbitrary x(0)>0, x(t) decreases and converges to zero.
(ii) if g( Sustainability 05 00440 i066) < 0, then Sustainability 05 00440 i065 = Sustainability 05 00440 i066 is the unique positive steady state and f(x) < 0 for all 0 < x Sustainability 05 00440 i065. If x(0) < Sustainability 05 00440 i065, then x (t) decreases and converges to 0, and if x(0) > Sustainability 05 00440 i066 then x (t) decreases again and now converges to Sustainability 05 00440 i065. If x(0) = Sustainability 05 00440 i065, then x (t) = Sustainability 05 00440 i065 for all t > 0.
(iii) If g( Sustainability 05 00440 i066) > 0, then equation (7) has two positive solutions, Sustainability 05 00440 i0651 < Sustainability 05 00440 i066 and Sustainability 05 00440 i0652 > Sustainability 05 00440 i066. Relation (6) implies that f(x) < 0 as x < Sustainability 05 00440 i0651 or x > Sustainability 05 00440 i0652, and f(x) > 0 if Sustainability 05 00440 i0651 < x < Sustainability 05 00440 i0652. Therefore if x(0) < Sustainability 05 00440 i0651, then x (t) decreases and converges to zero, if Sustainability 05 00440 i0651 < x(0) < Sustainability 05 00440 i0652, then x (t) increases and converges to Sustainability 05 00440 i0652, if x (0) > Sustainability 05 00440 i0652, then x (t) decreases and converges to Sustainability 05 00440 i0652. That is, Sustainability 05 00440 i0651 is locally unstable and Sustainability 05 00440 i0652 is locally asymptotically stable. If x(0) = Sustainability 05 00440 i0651 or x (0) = Sustainability 05 00440 i0652, then x (t) remains at that steady state level for all t > 0.

3. Model with Fixed Delay

The fixed delay T > 0 is assumed in the second term of the right hand side of equation (1), so we have the following equation:
Sustainability 05 00440 i006
Where
h(x) = βxγe-δx
The local asymptotic behavior of the trajectory can be examined by linearization. Let Sustainability 05 00440 i065 be a positive steady state. Then the linearized equation has the form
Sustainability 05 00440 i007
where xδ(t) is the deviation of x (t) from the steady state level. Looking for the solution in the usual form xδ(t) = eλtu, we have
Sustainability 05 00440 i008
which gives the characteristic equation
Sustainability 05 00440 i009
or
Sustainability 05 00440 i010
Lemma 1 Assume that |h´( Sustainability 05 00440 i065)| < α. Then Sustainability 05 00440 i065 is locally asymptotically stable.
Proof. Assume that Re λ ≥ 0. Then
|λ + α| ≥ α
and since
eλT = eT(Reλ)eiT(Imλ)
= eT(Reλ)(cos[T(Imλ)] + isin[T(Imλ)])
clearly
|eλT| ≥ 1
Therefore
Sustainability 05 00440 i011
implying that λ cannot be an eigenvalue. Q.E.D.
Notice that
h´(x) = βxγ-1e-δx(γ - δx)
and at the steady state
βxγe-δx = αx
implying that
βxγ-1e-δx = α
Therefore
Sustainability 05 00440 i012
so the characteristic equation (11) can be rewritten as
Sustainability 05 00440 i013
We also mention that the condition of Lemma 1 can be rewritten as
Sustainability 05 00440 i014
or equivalently
Sustainability 05 00440 i015
In the special case of γ = 1, this condition has the form
Sustainability 05 00440 i016
which is equivalent to relation
α < β < αe2
In order to give a complete stability analysis, we have to find the possible stability switches. Substituting any stability switch, λ= with ω > 0 into equation (11) yields
Sustainability 05 00440 i017
Separating the real and imaginary parts gives two equations,
Sustainability 05 00440 i018
and
Sustainability 05 00440 i019
Adding the squares of these equations gives
Sustainability 05 00440 i020
so
Sustainability 05 00440 i021
In order to have solution we have to assume now that
Sustainability 05 00440 i022
that is, (14) is violated with strict inequalities. Concerning this assumption, we can give the following interpretation. Let F(x) be y. Then
Sustainability 05 00440 i023
So this means that the absolute value of the elasticity of output with respect to capital is larger than unity. From (16) we have that if h΄( Sustainability 05 00440 i065)>0, then T(n) = T+(n) with
Sustainability 05 00440 i024
and if if h´( Sustainability 05 00440 i065) < 0, then T (n) = T-(n) with
Sustainability 05 00440 i025
and by (13) and (19),
Sustainability 05 00440 i026
so ( Sustainability 05 00440 i065) cannot be zero.
By selecting T as the bifurcation parameter and implicitly differentiating the characteristic equation with respect to T, we have
Sustainability 05 00440 i027
implying that
Sustainability 05 00440 i028
If λ = , then
Sustainability 05 00440 i029
with real part
Sustainability 05 00440 i030
Therefore if a steady state is unstable with T = 0, then it remains unstable for all T > 0, and if a steady state is asymptotically stable at T = 0, then this stability is lost at T = T(0) and cannot be regained later. That is, if the steady state is unstable without delay, then it remains unstable with any delay of positive length. If the steady state is asymptotically stable without delay, then it remains asymptotically stable until the delay reaches a certain threshold, and then becomes unstable and the stability cannot be regained later.
Taking, α = 1, β = 25, γ = 1 and δ = 1, we give an illustrative numerical example in Figure 1. The critical value γ - δ Sustainability 05 00440 i065 is denoted by
Sustainability 05 00440 i031
Introducing the notation z = γ - δ Sustainability 05 00440 i065 transforms the T -(0) curve to
Sustainability 05 00440 i032
and then the corresponding critical value of the delay is
Sustainability 05 00440 i033
In Figure 1(A), the steady state is locally asymptotically stable in the dark-gray region with z > -1 due to Lemma 1. It is also locally asymptotically stable in the light-gray region, which is under the critical curve T = T-(0) and it is unstable in the white region above the curve. Setting z = zc and increasing T along the vertical dotted line in Figure 1(A), we can see that the steady state loses stability at T = Tc. Further increasing T , as observed in Figure 1(B), generates complex dynamics through a quasi period-doubling bifurcation in which T increases from T c- 0.05 to 8.5 with an increment of 0.01 and the local maximum and minimum of the corresponding trajectory are plotted against each value of T.
Figure 1. Dynamics with α = 1, β = 25, γ = 1 and δ = 1.
Figure 1. Dynamics with α = 1, β = 25, γ = 1 and δ = 1.
Sustainability 05 00440 g001

4. Model with Continuously Distributed Delay

Assuming continuously distributed delays in the second term of equation (1) gives the following Volterra-type integro-differential equation:
Sustainability 05 00440 i034
where T > 0 is a positive parameter, the average delay and m ≥ 0 is an integer. The kernel function has the form
Sustainability 05 00440 i035
If m = 0, then the most current value has the highest weight, which is then decreasing exponentially. If m ≥ 1, then the most current value has zero weight which is then increasing until t - s = T, and is decreasing exponentially afterwards. As m increases, the weighting function becomes more peaked around t - s = T and as m → ∞, it converges to the Dirac delta function centered at t - s = T. If T → 0, then it also converges to the Dirac delta function.
A similar model is investigated by Fanti and Manfredi [11] where m = 2 is selected and the stability of the system with a cubic characteristic polynomial is examined based on the Routh-Hurwitz criterion. Stamova and Stamov [12] consider a generalized Solow model with endegenous labor growth and impulsive perturbations. Their stability analysis is based on the Lyapunov-Razumikhin sufficient stability conditions, which is a different approach than ours.
Linearizing equation (22), we have
Sustainability 05 00440 i036
where xδ(t) is the deviation of x (t) from the steady state level Sustainability 05 00440 i065. We are looking for the solution in the usual exponential form
xδ(t) = eλtu
then simple substitution shows that
Sustainability 05 00440 i037
Notice that by introducing the length of the delay as the new variable S = t - s in the integral, we see that
Sustainability 05 00440 i038
and by letting T → ∞, we have the characteristic equation
Sustainability 05 00440 i039
with
Sustainability 05 00440 i040
This equation can be rewritten as
Sustainability 05 00440 i041
Then similarly to the case of fixed delay we can prove the following result:
Lemma 2 Assume that |h´( Sustainability 05 00440 i065)| < α. Then Sustainability 05 00440 i065 is locally asymptotically stable.
Proof. Assume that Reλ≥ 0. Then
Sustainability 05 00440 i042
therefore
Sustainability 05 00440 i043
implying that λ cannot be a solution of equation (23). Q.E.D.
It is well known that the Routh-Hurwitz stability theorem provides necessary and sufficient conditions for all roots of a polynomial equation with real coefficients to have negative real parts. It is also known that it is difficult to locate the eigenvalues with analytic methods in general. However in some special cases, analytic results are still possible to obtain, as it will be next demonstrated.
Case I. T = 0
Assume first that T = 0, which reduces equation (22) with delays to equation (1) without delays. The asymptotic properties of this equation were already discussed earlier.
Case II. T > 0 and m = 0
Assume next that T > 0 and m = 0, when the kernel function becomes exponentially declining. Then characteristic equation (23) becomes quadratic,
Sustainability 05 00440 i044
or
Sustainability 05 00440 i045
If γ ≤ 1, then all coefficients are positive with a positive steady state, which is locally asymptotically stable. Assume next that γ > 1. If
Sustainability 05 00440 i046
then the constant term is zero indicating that one eigenvalue is zero and the other is negative. So Sustainability 05 00440 i065 is marginally stable in the linearized model, so no conclusion can be drawn about its asymptotical behavior in the nonlinear model. If Sustainability 05 00440 i065 < (γ - 1)/δ, then Sustainability 05 00440 i065 is unstable and if Sustainability 05 00440 i065 > (γ - 1)/δ, then Sustainability 05 00440 i065 is locally asymptotically stable.
Case III. T > 0 and m = 1
Assume now that T > 0 and m = 1, when the shape of the kernel function takes a bell-shaped form. Then we have a cubic characteristic polynomial:
Sustainability 05 00440 i047
or
Sustainability 05 00440 i048
If γ ≤ 1, then all coefficients are positive at a positive steady state. If γ > 1, then we can consider three cases. If
Sustainability 05 00440 i070
then zero is an eigenvalue and the other two eigenvalues have negative real parts implying that Sustainability 05 00440 i065 in the linearized system is marginally stable. Therefore no conclusion can be drawn about the stability of Sustainability 05 00440 i065 in the nonlinear system. If Sustainability 05 00440 i065 < (γ - 1)/δ, then the constant term is negative, so Sustainability 05 00440 i065 is unstable. If Sustainability 05 00440 i065 > (γ - 1)/δ then all coefficient of (25) are positive. In this case and when γ ≤ 1 the Routh-Hurwitz criterion implies that the real parts of the eigenvalues are negative if and only if
Sustainability 05 00440 i049
which can be reduced to a quadratic inequality in T:
Sustainability 05 00440 i050
For the sake of simplicity, we re-introduce the notation = γ - δ Sustainability 05 00440 i065. If z ≥ 0, then this inequality holds implying the asymptotical stability of the steady state.So we can assume that z < 0.The discriminant of the left hand side of inequality (26) is
D = z(z + 8)
If z < -8, then D > 0, so the left hand side of (26) has two roots
Sustainability 05 00440 i051
which are positive and Sustainability 05 00440 i067 < Sustainability 05 00440 i068. Notice that Sustainability 05 00440 i067 Sustainability 05 00440 i068 = 1/α2 and (26) holds if and only if T < Sustainability 05 00440 i067 or T > Sustainability 05 00440 i068 when Sustainability 05 00440 i065 is locally asymptotically stable. If Sustainability 05 00440 i067 < T < Sustainability 05 00440 i068, then (26) is violated, so Sustainability 05 00440 i065 is unstable. If z = -8, then D = 0 and there are equal roots
Sustainability 05 00440 i052
so if T ≠ 1/α, then Sustainability 05 00440 i065 is locally asymptotically stable.
If -8 < z < 0, then D < 0, so (26) holds and Sustainability 05 00440 i065 is asymptotically stable. The instability region is shown in Figure 2(A) where z is the horizontal axis and T is the vertical axis. If we start with a very small value of T with any given z < -8, then Sustainability 05 00440 i065 is asymptotically stable. If we gradually increase T then Sustainability 05 00440 i065 remains asymptotically stable until it reaches the critical value Sustainability 05 00440 i067, when the steady state becomes unstable. It remains unstable until Sustainability 05 00440 i068 when stability is regained, and Sustainability 05 00440 i065 remains asymptotically stable for all T > Sustainability 05 00440 i068.
If the steady state is unstable without delay, then it remains unstable with continuous delay with any T and m. If it is asymptotically stable without delay, then either it remains asymptotically stable with all T and m, or loses stability at a certain value of the average delay T and stability is regained with an even larger value of T and the steady state remains asymptotically stable afterwards. In such cases small and large average delays lead to asymptotically stable steady states.
We will next show that at the critical values Sustainability 05 00440 i067 and Sustainability 05 00440 i068, Hopf bifurcation occurs giving the possibility of the birth of limit cycles. We select T as the bifurcation parameter. At the critical values (26) is satisfied with equality, so
Sustainability 05 00440 i053
and the characteristic equation (25) can be rewritten as
Sustainability 05 00440 i054
showing that there is a negative eigenvalue
Sustainability 05 00440 i055
and a pair of pure complex eigenvalues
Sustainability 05 00440 i056
Consider λ as a function of the bifurcation parameter T and differentiate implicitly equation (25) to have
Sustainability 05 00440 i057
By simple calculation we can see that at λ = ± ,
Sustainability 05 00440 i058
with real part
Sustainability 05 00440 i059
Since Sustainability 05 00440 i067 Sustainability 05 00440 i068 = 1/α2, at T = Sustainability 05 00440 i067 the value of d(Reλ)/dT changes from negative to positive showing the loss of stability, and if T = Sustainability 05 00440 i068, then d(Reλ)/dT changes from positive to negative indicating that stability is regained. Since at both critical values d(Reλ)/dT ≠ 0, at both values Hopf bifurcation occurs giving the possibility of the birth of limit cycles.
We perform numerical simulations to illustrate the results obtained above. In Figure 2(A), the steady state is locally asymptotically stable in the dark-gray region with z > -1 due to Lemma 2. It is also locally asymptotically stable in the light gray region and unstable in the white region when z < - 1.The appearance and disappearance of a limit cycle can be observed in Figure 2(B) where we take α = 1, β = e13, γ = 1 and δ = 1 implying z = - 12,
Sustainability 05 00440 i060
Under these specifications, the Volterra-type integro-differential equation (22) can be written as a 3D system of differential equations,
Sustainability 05 00440 i061
where
Sustainability 05 00440 i062
and
Sustainability 05 00440 i063
When T increases from Sustainability 05 00440 i067 - 0.1 to Sustainability 05 00440 i068 + 0.3 with an increment of 0.01, the steady state loses stability at point A and regains stability at point B In Figure 2(B), the local maximum and local minimum of a trajectory generated by the 3D system are depicted against each value of T indicating the birth of a limit cycle for Sustainability 05 00440 i067 < T < Sustainability 05 00440 i068.
Figure 2. Dynamics with T > 0 and m = 1.
Figure 2. Dynamics with T > 0 and m = 1.
Sustainability 05 00440 g002
The cases of m ≥ 2 result in fourth or larger degree polynomial equations. The stability of the steady states can be examined similarly, but the mathematical details become much more complicated. It can be mathematically confirmed that as m → ∞, equation (23) converges to the characteristic equation (11) of the model with fixed delay. In particular, if m → ∞, then expression
Sustainability 05 00440 i064
converges to eλT. For larger values of m dynamics generated by the differential equation with continuously distributed time delay is similar to dynamics generated by the differential equation with fixed time delay.

5. Conclusions

In this paper, a special neoclassical growth model was introduced and examined. A mound-shaped production function for capital growth was assumed in the dynamic equation. Zero is always a steady state, and depending on model parameters there is either no positive steady state, or one, or two positive steady states. A complete steady state analysis was followed by the derivation of stability conditions. By introducing fixed delay we demonstrated that stability can be lost at a certain value of the delay and the equilibrium remains unstable afterwards. In the case of continuously distributed delays it has been shown how stability can be lost at a certain value of the average delay and by further increasing the average delay it can be regained. At the critical values, Hopf bifurcation occurs giving the possibility of the birth of limit cycles. In our further study, more complex kernel functions will be considered and their effect on the asymptotic behavior of the steady state will be examined.
It is now well-known that any discrete-time dynamic model can generate complex dynamics involving chaos. However, in an aggregate model like the neoclassical growth model, it is natural to treat economic phenomena as continuous since the variables to be examined are the outcomes of a great number of decisions taken by different agents at different points of time. This is the main reason why this paper considers the neoclassical growth model in the continuous-time framework. Implication of the analysis is the following. Coexistence of nonlinearity and delay in production can be a source of persistent fluctuations in the capital-labor ration.

Acknowledgments

The authors are grateful to two anonymous referees for their constructive suggestions and helpful comments and highly appreciate financial supports from the Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research (C) 24530202) and Chuo University (Grant for Special Research and Joint Research Grant 0981).

Conflict of Interest

The authors declare no conflict of interest.

References and Notes

  1. Day, R. Irregular growth cycles. Am. Econ. Rev. 1982, 72, 406–414. [Google Scholar]
  2. Day, R. The Emergence of chaos from classical economic growth. Q J. Econ. 1983, 98, 203–213. [Google Scholar]
  3. Day, R. Complex Economic Dynamics: An Introduction to Dynamical Systems and Market Mechanism; MIT Press: Cambridge, MA, USA, 1994. [Google Scholar]
  4. Puu, T. Attractions, Bifurcations and Chaos: Nonlinear Phenomena in Economics, 2nd ed; Springer-Verlag: Berlin, Heidelberg, New York, NY, USA, 2003. [Google Scholar]
  5. Bischi, G.I.; Chiarella, C.; Kopel, M.; Szidarovszky, F. Nonlinear Oligopolies: Stability and Bifurcation; Springer-Verlag: Berlin, Heidelberg, New York, NY, USA, 2010. [Google Scholar]
  6. Li, T.-Y.; Yorke, J.A. Period three implies chaos. Am. Math. Mon. 1975, 82, 985–992. [Google Scholar]
  7. Complexity in Economics I, II, III. In The International Library of Critical Writings in Economic Series; Rosser, J.B. (Ed.) Edward Elgar: Northampton, MA, USA, 2004.
  8. Solow, R. A Contribution to the theory of economic growth. Q. J. Econ. 1956, 70, 65–94. [Google Scholar] [CrossRef]
  9. Swan, T. Economic growth and capital accumulation. Econ. Rec. 1956, 32, 334–361. [Google Scholar] [CrossRef]
  10. Matsumoto, A.; Szidarovszky, F. Delay differential neoclassical growth model. J. Econ. Behav. Organ. 2011, 78, 272–289. [Google Scholar] [CrossRef]
  11. Fanti, L.; Manfredi, P. The Solow’s model with endogenous population: A neoclassical growth cycle model. J. Econ. Dev. 2003, 28, 103–115. [Google Scholar]
  12. Stamova, I.M.; Stamov, A.G. Impulsive control on the asymptotic stability of the solutions of a Solow model with endogenous labor growth. J. Franklin Inst. 2012, 349, 2704–2716. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Matsumoto, A.; Szidarovszky, F. Asymptotic Behavior of a Delay Differential Neoclassical Growth Model . Sustainability 2013, 5, 440-455. https://doi.org/10.3390/su5020440

AMA Style

Matsumoto A, Szidarovszky F. Asymptotic Behavior of a Delay Differential Neoclassical Growth Model . Sustainability. 2013; 5(2):440-455. https://doi.org/10.3390/su5020440

Chicago/Turabian Style

Matsumoto, Akio, and Ferenc Szidarovszky. 2013. "Asymptotic Behavior of a Delay Differential Neoclassical Growth Model " Sustainability 5, no. 2: 440-455. https://doi.org/10.3390/su5020440

APA Style

Matsumoto, A., & Szidarovszky, F. (2013). Asymptotic Behavior of a Delay Differential Neoclassical Growth Model . Sustainability, 5(2), 440-455. https://doi.org/10.3390/su5020440

Article Metrics

Back to TopTop