Journal of Nonlinear Mathematical Physics

In Press, Corrected Proof, Available Online: 21 January 2021

On Ramsey Dynamical Model and Closed-Form Solutions

Gülden Gün Polat1, Teoman Özer2, *
1İstanbul Kültür University, Faculty of Science and Letters, Mathematics and Computer Sciences, 34158 Bakırköy, İstanbul-Turkey
2İstanbul Technical University, Division of Mechanics, Faculty of Civil Engineering, 34469 Maslak, İstanbul-Turkey
*Corresponding author. Email:
Corresponding Author
Teoman Özer
Received 2 July 2020, Accepted 2 January 2021, Available Online 21 January 2021.
DOI to use a DOI?
Ramsey dynamical model, economic growth models, Lie point symmetries, Prelle-Singer approach, Jacobi last multiplier, Hamiltonian dynamics, closed-form solutions

This study focuses on the analysis of Ramsey dynamical model with current Hamiltonian defining an optimal control problem in a neoclassical growth model by utilizing Lie group theory. Lie point symmetries of coupled nonlinear first-order ordinary differential equations corresponding to first-order conditions of maximum principle are analyzed and then first integrals and corresponding closed-form (analytical) solutions are determined by using Lie point symmetries in conjunction with Prelle-Singer and Jacobi last multiplier methods. Additionally, associated λ-symmetries, adjoint symmetries, Darboux polynomials, and the properties of the model are represented.

© 2021 The Authors. Published by Atlantis Press B.V.
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (


The concept of the development of an economic growth model is one of the most important research areas in the theory of economics. The main purpose of an economic growth model as an applied dynamical system is to investigate the active forces to determine the growth and economic development. The classical models are based on the fact that economic growth depends on investments and improving productive capacity in an economy. However, in the case of the neoclassical economic growth models, economic growth depends on land, labor, and capital. In the literature, there are well-known economic growth models. For example, in the Solow model [27], agents and the planner follow up a straightforward linear rule for consumption and investments. However, in the case of the Ramsey model, the agents and the planner follow up consumption and investments to maximize their welfare. As early as 1928 a sophisticated model of a society’s optimal saving was published by the British mathematician Frank Ramsey [24]. In the Ramsey model, there is a finite number of agents with an infinite time horizon; further, these agents are completely alike and thus it is a representative agent model. The Ramsey framework can be formulated in discrete time as well as in continuous time. In this study, the continuous-time version as did Ramsey’s original contribution is considered.

The essence of the optimal control problem is to choose, for the entire time interval, a path c that maximizes υ(k(t), c(t), t) subject to the equation of the motion governing k(t). These optimization problems are commonly defined as:

maxV=0Tυ(k(t),c(t),t)dt,k˙(t)=g(k(t),c(t),t),t[0,T],k(T)e-R(T)T0. (1.1)

In this formulation, the objective function is the integral over a predetermined interval [0, T] of the payoff function υ. This payoff function depends, at each time instance, on the control variable c(t), which the planner can directly control, the state variable k(t), which the planner cannot directly control but is affected by c(t) via its equation of motion, and the time t. R(s) denotes the average discount rate between time zero and time t. For example, in an optimal savings problem, the control and state variables can be considered as the consumption and the level of assets, respectively. In general, the control problem involves several control and state variables. The constraint function expressing the change in the state variable may depend on the state variable, control variable, and time. The constraints state that: (1) At each moment, the change in state variable depends on the state variable itself and/or control variable and/or time; (2) the initial level of the state variable is given; (3) the discounted value of the state variable at the end of the planning horizon has to be weakly positive. If the planning horizon T is infinite, the last constraint can equivalently be stated as requiring that k(T) ≥ 0.

The dynamical optimization problems can be analyzed by utilizing Hamiltonian, which is set in two different ways. The first one can be written with the discount factor ahead and then a multiplier to every constraint is assigned. It is known as the present value Hamiltonian since it includes the discount of the future. Furthermore, the multipliers are called costate variables since it is associated with the constraint of the accumulation of the state variable. The second one called current value Hamiltonian can be set without taking into account the discount factor. In this study, the analysis of the Ramsey model with current value Hamiltonian is carried out by considering the first-order conditions (FOCs).

In the literature, there are some qualitative and quantitative methods to investigate closed-form and numerical solutions for the optimal control models including current value Hamiltonian. Hamiltonian dynamical systems in economics were introduced and its potential usefulness to economic theory was presented in the study [3]. A method for solving a family of four-dimensional nonlinear modified Hamiltonian dynamic systems in closed form was proposed in the studies [9,25,26]. A closed-form solution to the Ramsey model with logistic population growth was investigated in [8]. Furthermore, a closed-form solution to the transitional dynamics of the Lucas-Uzawa model was given in [5]. However, it can be said that there is no general solution procedure for these kinds of optimal control problems.

It is a fact that the theory of Lie groups is one of the most important methods for the investigation of invariant solutions of systems for differential equations [1,2,14,16,17,1923]. This approach can also be used effectively for the analysis of optimal controls problem of dynamical systems in the economic growth models. In the literature, some economic growth models are analyzed by the partial Hamiltonian and partial Lagrangian approaches [15,16]. The integrability and linearization properties of the coupled ODEs are evaluated based on the Prelle-Singer procedure in the study [4]. Besides, the JLM approach was built up by Jacobi [10] for the application of the multipliers to classical mechanics problems and it has a direct relationship with Lie point symmetries. It is important to mention that all these methods have important properties for the investigation of solutions of partial or ordinary differential equations [18,19]. Here, we deal with FOCs of the Ramsey model representing by the coupled nonlinear ordinary differential equations (ODEs) by using the aforementioned methods, first, by determining its Lie point symmetry properties of the Ramsey model. However, it is a fact that the determination of Lie point symmetries of the model requires some tough mathematical computations because of its nonlinearity and order of the system. In this study, we represent the theory of Lie group as a general analysis and solution procedure for the optimal control problems with economic growth models by using its relations with Prelle-Singer (PS) and Jacobi last multiplier (JLM) methods.

The rest of the paper is arranged as follows: In section 2, we introduce the preliminaries on definitions for current value Hamiltonian, as well as mathematical methods such as Lie point symmetries, Prelle-Singer, Jacobi last multiplier methods, and related mathematical relations. Section 3 deals with the analytical analysis and shows the results for Ramsey’s economic growth model with current value Hamiltonian including first integrals and closed-form solutions as well as λ-symmetries and Darboux polynomials. From the mathematical point of view, analytical solutions of the coupled system of nonlinear ordinary differential equations (ODEs) are investigated. The relationships between PS procedure and some analytical methods such as Lie point symmetries, Jacobi last multiplier (JLM), λ-symmetries, and Darboux polynomials [12] are indicated. In section 3, some conclusions and discussions on the analysis of the model are represented.


In general, an optimal control problem is defined as

max𝒡=0TF(t,qi,ci)e-ρtdt (2.1)
subjecttoq˙i=fi(t,qi,ci),i=1,,n, (2.2)
where q˙i=dqi/dt , ci = (c1, ... , cm) is the control vector, mn or m > n with the proper boundary conditions, t is the independent variable representing time and the pair (qi, pi) = (q1, ... , qn, p1, ... , pn) are the space coordinates, where q1, ... , qn are called the state variables and p1, ... , pn are called the costate variables. The integrand in the integral (2.1) contains the discount factor eρt. Furthermore, the investigation of the optimal control problems in economic growth models can be considered by taking into account Hamiltonian function, namely current value Hamiltonian.

Definition 2.1.

The current value Hamiltonian is defined by

H(t,qi,ci)=F(t,qi,ci)+pifi(t,qi,ci). (2.3)

Remark 1.

The necessary conditions for optimal control based on the Pontryagin maximum principle for current value Hamiltonian H are

Hci=0, (2.4)
q˙i=Hpi, (2.5)
p˙i=-Hqi+Γi,i=1,,n, (2.6)
where Γi generally is taken as a nonzero function of variables t, pi, and qi.

The conditions (2.42.6) derive optimality with all variables starting out from any given initial position. In this way, the optimal control problems and associated optimal initial values become much simpler for determining the optimality.

In the literature, Ramsey’s neoclassical growth model based on a dynamical system [24] is one of the principal models in the case of the theory of dynamic macroeconomics, the investigation of the related optimal control problem is considered by taking into account Hamiltonian. The extended types of this growth model are considered to analyze some important macroeconomic problems. In general, the Ramsey neoclassical growth model as a dynamical system dealing with representative consumer’s utility maximization problem is defined as

Max0e-rtc1-σdt,σ0,1,r0 (2.7)
subject to the capital accumulation equation
k˙(t)=kβ(t)-δk(t)-c(t),k(0)=k0,0<β<1,δ0 (2.8)
where c = c(t) is consumption per person, k(t) is capital-labor ratio, β is capital share, δ is depreciation rate, and r is time preferences. The intertemporal elasticity of substitution is given by 1 and k0 is the capital stock. In the economic growth models with a dynamical system, the maximum principle deals with the system of coupled and nonlinear first-order ordinary differential equations in terms of state and costate variables. In this approach, the main aim is to maximize the current value Hamiltonian function via control variables.

The current value Hamiltonian function for the model (2.7) by using Definition 2.1 from the equation (2.3) can be written as

H(t,c,k,p)=c1-σ+p(kβ-δk-c) (2.9)
in which p = p(t) is the costate variable. Secondly, the necessary first-order conditions for the optimal control problem are obtained from equations (2.42.6) of the form
p(t)=(1-σ)c(t)-σ,k˙(t)=k(t)β-δk(t)-c(t),p˙(t)=-p(t)(βk(t)β-1-δ)+Γ, (2.10)
where the “overdot” represents the first-order derivatives concerning time parameter t and it can be checked that Γ = rp(t). Besides, we can transform the nonlinear coupled system (2.10) into the following coupled first-order ODEs by eliminating control variable c(t) by using the application of the implicit function theorem based on the fact that the second derivative of the Hamiltonian function with respect to variable c(t) is an invertible map
p˙(t)=-p(t)(βk(t)β-1-δ)+rp(t),k˙(t)=k(t)β-δk(t)+(-1)1+1σ(σ-1)1σp(t)-1σ,σ0,1. (2.11)

Remark 2.

Besides, it is important to mention about some comments on the difference between the Ramsey model and the economic growth model with logarithmic utility function analyzed in the study [22] In general, especially, the definition of Ramsey model includes a parameter σ, which means “the intertemporal elasticity of substitution is given by 1/σ (σ ≠ 0, 1)” in which the economic growth model with logarithmic utility function has not a similar term. Hence, because of this parameter, the economic growth model with logarithmic utility function and Ramsey dynamical system are quite different not only from mathematical point of view but from economic growth models point of view.

Remark 3.

Let us suppose that a system of two coupled first-order ODEs of the form

p˙=ϕ1(t,p,q),q˙=ϕ2(t,p,q), (2.12)
in which ϕi, i = 1, 2 are analytic functions. The PS method is based on the solution of the following set of equations
St+ϕ1Sp+ϕ2Sq=-ϕ2q+(ϕ2q-ϕ1p)S+ϕ1qS2, (2.13)
Kt+ϕ1Kp+ϕ2Kq=-K(Sϕ1q+ϕ2q), (2.14)
Kp=SKq+KSq. (2.15)
where R(t, p, q) and K(t, p, q) are integrating factors of the system and S(t, p, q) is called a null form. Furthermore, the associated first integral is obtained as
I=r1+r2-[K+ddq(r1+r2)], (2.16)
r1=(Rϕ1+Kϕ2)dtandr2=-(R+ddp(r1))dp. (2.17)

Remark 4.

Assume that a partial differential equation (PDE) is given by in the form

Af=i=1na1(x1,,xn)fxi=0, (2.18)
and its equivalent Lagrange’s system is
dx1a1=dx2a2==dxnan. (2.19)

The system of first-order differential equation of the form

x˙=f(x), (2.20)
is given. Jacobi last multiplier of the above system, M, can be determined by the solution of the following partial differential equation
ddtlogM+divf=0 (2.21)

Based on the original formulation of the Jacobi last multiplier method, complete knowledge of the system (2.18) or (2.19) should have been available. Because of the existence of a solution/first integral is consequent upon the existence of symmetry, an alternative formulation in terms of symmetries was presented by Sophus Lie [11]. If n − 1 symmetries of (2.18)/(2.19) are known

Γi=j=1nξij(x1,,xn)xj,i=1,,n-1, (2.22)
and JLM, that is equal to M = Δ−1, is calculated with the following determinant
Δ=det[a1anξ1,1ξ1,nξn-1,1ξn-1,n], (2.23)
in which Δ ≠ 0 [17].

Remark 5.

Let ξ(t, p, q), η1(t, p, q), and η2(t, p, q) be Lie symmetries of (2.12) and then the characteristic notation Q [19] can be written as

Q1=η1-p˙ξ=η1-ϕ1ξ,Q2=η2-q˙ξ=η2-ϕ2ξ. (2.24)

Thus, a relation with S null function (2.13) is given by

S=Q1Q2. (2.25)

If the Lie symmetries of the equation (2.12) are known, and then it is clear that (2.13) is satisfied by using the relation (2.25). Hence, via the equations (2.14) and (2.15), the process of PS method is achieved. Now let us suppose that Λ1 and Λ2 are adjoint symmetries of (2.12), then they must satisfy the adjoint equation for the linearized symmetry conditions

D[Λ1]=-(Λ1ϕ1p+Λ2ϕ2p),D[Λ2]=-(Λ1ϕ1q+Λ2ϕ2q). (2.26)

The function K an integrating factor satisfying the second equation given by (2.14) can be determined by using the relations between Jacobi last multipliers, Lie point symmetries (2.33), and PS method. However, it is a fact that there are no standard techniques to determine first integrals of ODEs, yet there is several special methods widely discussed in the literature. One of them is called Darboux polynomial method. There is a simple relation between Darboux polynomials and Jacobi’s last multipliers. The determining equation for the Darboux polynomial F [6] of equation (2.12) is defined by

D[F]=αF, (2.27)
in which α = (ϕ2q + ϕ1p) represents the corresponding cofactor. The relation between Jacobi last multiplier and Darboux polynomial [12] is given by
F=M-1. (2.28)

Remark 6.

It is clear that, given a system of first-order differential equations like (2.20) (allowing f to depend also on t) the determining equations are the solution of the following system of first order PDEs:

ηit+i=1Nηkxkfk+fi(ξit+i=1Nξkxkfk)=fit+i=1Nfixkηk,i=1,,N, (2.29)
where X=ξt+k=1Nηkxk is the infinitesimal generator of the Lie symmetry. For every N ≥ 1 the solution of this system will depend on arbitrary functions, so that a system of first order differential equations admits infinite-dimensional Lie algebra.

Proposition 2.1.

The full symmetry group of (2.11) is always infinite-dimensional. However, the symmetry X=t presents one-dimensional Lie subalgebra if there are no any constraints on the parameters β, σ , δ, and r. On the other hand, the Lie subalgebra L4 arises only for the given combination of parameters such as σ =−1, β=12 , δ=-2r3 .

Proof. In order to analyze the solutions of the system (2.11), firstly Lie point symmetries of the coupled-system are examined. In order to determine the Lie point symmetries of the system (2.11), let us consider the vector field X=ξ(t,p,k)t+η1(t,p,k)p+η2(t,p,k)k . The first prolongation of this vector field written according to the prolongation formula [19] is applied separately to the equations in the system (2.11). Thus, the overdetermined systems of PDEs, called determining equations, are obtained in terms of infinitesimal functions of the infinitesimal generator as follow

(-1)1+1σkp1σ(σ-1)1ση1+pσ(-p2/σ(kββ-kδ)η2-k(p1σ(kδ-kβ)+(-1)1σ(σ-1)1σ)2ξk+p1σ(p(-kββ+k(r+δ))(p1σ(kδ-kβ)+(-1)1σ(σ-1)1σ)ξp+(-1)1σk(σ-1)1σ(ξt-η2k)+p1σ(k(kδ-kβ)ξt+k(kβ-kδ)η2k-kβpβη2p+k(p(r+δ)η2p+η2t))))=0, (2.30)
-kp1σ(-kββ+k(r+δ))η1+kβp1+1σ(β-1)βη2-k2+βp1+1σrξk+k1+2βp1+1σβξk-k2+βp1+1σδξk+k3p1+1σrδξk-k2+βp1+1σβδξk+k3p1+1σδ2ξk-k2p2+1σr2ξp+2k1+βp2+1σrβξp-k2βp2+1σβ2ξp-2k2p2+1σrδξp+2k1+βp2+1σβδξp-k2p2+1σδ2ξp-k2p1+1σrξt+k1+βp1+1σβξt-k2p1+1σδξt+k2+βp1ση1k-k3p1ση1k+(-1)1σk(σ-1)1σ(p(-kββ+k(r+δ))ξ1k-kη1k)+k2p1+1σrη1p-k1+βp1+1σβη1p+k2p1+1σδη1p+k2p1ση1t. (2.31)

Based on the above Remark 6, the solution of this system of first-order differential equations as determining equations depends on an arbitrary function and then it admits an infinite-dimensional Lie algebra. On the other hand, in order to obtain a finite-dimensional subalgabra, from the direct integration of the determining equations (2.30) and (2.31), the following algebraic equations including the model parameters, namely β, δ, σ , r and the group parameters, namely s1, s2, and s3

s1(-2+β)βδ2(-r2(-1+β)-r(2-2β+β2)δ+(1-β+β2)δ2)=0,s2(-1)((1/β))r3(1/(-1+β))(-4+1/β)(-1+2β)=0,s3(-1+β)δ(r+δ-βδσ)=0,δ0,β1, (2.32)
are obtained. If there are no any constraints on the parameters and then we have s1 = 0, s2 = 0, and s3 = 0 in (2.32) and then it can be shown that the analytical integration of the associated determining equations yields a finite subalgebra L1 with basis operator of the form X=t .

In addition, the solution of (2.32) gives σ =−1, β=12 , δ=-2r3 for s1 ≠ 0, s2 ≠ 0, s3 ≠ 0 (for a nontrivial symmetry group). Consequently, under these constraints, the solution of determining equations (2.30) and (2.31) yields a four-dimensional infinite subalgabra given by

X1=-32rt,X2=ert3t+r3ert3pp+2r3ert3kk,X3=ert3(3+2rk)3kp+ert3k,X4=e2rt3kp+e2rt3k. (2.33)

Remark 7.

It is clear that the use of Lie point symmetry X=t with Prelle-Singer procedure and Jacobi last multipliers method does not yield functions S and K. Therefore, Lie point symmetry X=t does not give rise to an invariant by using Prelle-Singer procedure for the system (2.11) without any constraints on the model parameters.

Remark 8.

It is important to mention that the relation between Jacobi last multiplier and Lie symmetries enable to us to determine first integrals. Since this case has four Lie symmetries, then six Jacobi last multipliers of this system can be obtained. Hence, by using the ratio of two multipliers, thirty first integrals of the system are achieved, it is clear that some of these first integrals should be equivalent. Via the first integrals for this case, K and S functions can be found based on Remark 4. As a result, six different S null forms by using first integrals that obtained by JLM method are determined. Four of them overlap the null forms S1, ... , S4 that are derived by Lie symmetries and other two of them S5 and S6, which cannot be obtained directly from Lie symmetries are given at the Table 1.

Generators Null Forms Integration Factors λ-symmetries and First Integrals
X1 S1(t,p,k)=k(3p-6k-4kr)p(2rk-3) K1(t,p,k)=Λ2=-(4kr2-9+pr(3-2rk))3p(p-4k)2rk λ1(t,p,k)=r3+3(p-4k)12k-8rk3/2
R1(t,p,k)=Λ1=4(pr-3-2rk)(3p-6k-4kr)3rp2(p-4k)2 V1=p-(k(3p-6k-4kr)p(2rk-3))k
X2 S2(t,p,k)=2kp-k K2(t,p,k)=Λ2=-8e2rt/3p(p-4k)2k λ2(t,p,k)=p4k-1k+r3
R2(t,p,k)=Λ1=8e2rt/3(p-2k)p2(p-4k)2 V2=p-(2kp-k)k
X3 S3(t,p,k)=3k3+2rk K3(t,p,k)=Λ2=-12r2e2rt/3(3+2rk)k((3+2rk)2-6pr)2 λ3(t,p,k)=112(-6k+4r+9p3k+2rk3/2)
R3(t,p,k)=Λ1=36e2rt/3r2(3+2rk)2-6pr)2 V3=p-(3k3+2rk)k
X4 S4(t,p,k)=-k K4(t,p,k)=Λ2=-r2e2rt/3k(3-pr+2rk)2 λ4(t,p,k)=3p-6k+4kr12k
R4(t,p,k)=Λ1=r2e2rt/3(3-pr+2rk)2 V4=p+kk
S5(t,p,k)=-k(-6p+12k+(p2-4pk+8k)r)p(pr-6) K5(t,p,k)=Λ2=-2ert3(pr-6)pr(p-4k)2k λ5(t,p,k)=72k+p(-18+3r(p-6k)+4r2k)12k(pr-6)
R5(t,p,k)=Λ1=2ert3(-6p+12k+r(p2-4pk+8k))p2r(p-4k)2 V5=p+(k(-6p+12k+(p2-4pk+8k)r)p(pr-6))k
S6(t,p,k)=k(4kr2-9)4kr2+9-4kr(pr-3) K6(t,p,k)=Λ2=ert34kr2+9-4kr(pr-3)k((2rk+3)2-6pr)2 λ6(t,p,k)=16(2r-3k)+p(9-4kr2)4k(9+4kr2-4rk(pr-3))
R6(t,p,k)=Λ1=3ert3r(4kr2-9)((2rk+3)2-6pr)2 V6(t,p,k)=p-(k(4kr2-9)4kr2+9-4kr(pr-3))k
Table 1

Prelle-Singer Procedure of the system (2.11) for σ =−1, β=12 , δ=-2r3

It is clear that due to the presence of symmetries (2.33), such an invariant can be derived in a purely algebraic way from JLM. For instance, from JLMs M13=6re-rt/36pr-(3+2rk)2 and M24 = ert corresponding symmetries X1, X3 and X2, X4 respectively, the first integral

I=6re2rt/36pr-(3+2rk)2, (2.34)
is determined and the analytical solution for the coupled system (2.11) using the invariant (2.34)
p(t)=3r+23c2rert/3+4c2ert/3+9c1-6re2rt/3c1r2,k(t)=c1ert/3+9c1-6re2rt/34c1r2, (2.35)
can be obtained, where c1 and c2 are arbitrary constants. Furthermore, Darboux polynomials of the system for this case (2.11) are
F1(t,p,k)=(p2-4pk)ert/34,F2(t,p,k)=ert/3(6pr-(3+2rk)2)6r,F3(t,p,k)=e2rt/3(-3+pr-2rk)2r,F4(t,p,k)=e2rt/3(3-pr+2rk)3,F5(t,p,k)=ert,F6(t,p,k)=2rert3, (2.36)
which satisfy the condition (2.27).

Proposition 2.2.

From the analytical integration of the determining equations (2.30) and (2.31), the following finite subalgebra corresponding to the model without constraints

X˜1=t. (2.37)
can be found and for the constraint σ=r+δβδ , the system (2.11) admits two-dimensional finite subalgebra L2.

Proof. First of all, similarly, from the direct integration of the determining equations (2.30) and (2.31), we have the following relation for the model parameters

s(β-1)δ(r+δ-βδσ)=0, (2.38)
where s is a group parameter. It can be shown that s = 0 gives the trivial symmetry group as given in (2.37) and s ≠ 0 yields a nontrivial symmetry group and for the conditions δ ≠ 0, and β ≠ 1, thus we have a constraint σ=r+δβδ from (2.38). By similar calculations in the proof of the Proposition 2.1 for the constraint σ=r+δβδ , the analytical integration of the determining equations (2.30) and (2.31) yields Lie subalgebra with basis operators of the form
X˜2=etδ(β-1)δ(β-1)t+etδ(β-1)p(r+δ)δ(β-1)p+etδ(β-1)k(β-1)k. (2.39)

Remark 9.

For the condition σ=r+δβδ , Lie point symmetry X˜1 does not produce an analytical solution for the system (2.11) by using PS aproach.

For the condition σ=r+δβδ and Lie point symmetry X˜1=t , the null function

S˜1(t,p,k)=(-1)1σp-1σ(σ-1)1σ+kδ-kβp(r-βkβ-1+δ), (2.40)
is determined, which clearly satisfies the equation in (2.13). The determining equation (2.14) requires some demanding calculations for obtaining the solution of the function K. For this purpose, let us consider the ansatz of the form
K˜(t,p,k)=Sd(A(t,p)+B(t,p)k)z, (2.41)
where Sd is the denominator of S˜1 , and A(t, p), B(t, p) are functions of t and p, and z is a constant. Since the function K is in a rational form, the differentiation or integration the form of the denominator remains the same. On the other hand, the power of the denominator decreases or increases by a unit order from that of the initial one [4]. Using (2.41) for function K by considering the expressions
ϕ1(t,p,k)=-p(βkβ-1-δ)+pr,ϕ2(t,p,k)=kβ-δk+(-1)1+1σ(σ-1)1σp-1σ, (2.42)
it can be checked that the equation (2.13) is satisfied and the solution of the equation (2.14) in terms of the function K using the ansatz form (2.41) is
K˜(t,p,k)=α1-ze-rtp(r-βkβ-1+δ), (2.43)
where α1 is a constant. In addition, the compatibility condition (2.15) gives α 1zertr=0 . However, for the time preference parameter r, the condition r ≠ 0 means α1 = 0(z < 0), that is K = 0. Thus, we conclude that the symmetry X˜1 does not produce a nontrivial solution for the system (2.11). Furthermore, λ-function and λ-symmetry are determined as
λ˜1(t,p,k)=r-βkβ-1+δ+kβ-2(β-1)β((-1)1σp-1σ(σ-1)1σ+kδ-kβ)r-βkβ-1+δ,V˜1=p-((-1)1σp-1σ(σ-1)1σ+kδ-kβp(r-βkβ-1+δ))k. (2.44)

Remark 10.

For the condition σ=r+δβδ , Lie point symmetry X˜2 produces a nontrivial analytical solution for system (2.11) via PS approach.

To show this, let us consider the Lie point symmetry X˜2=etδ(β-1)δ(β-1)t+etδ(β-1)p(r+δ)δ(β-1)p+etδ(β-1)k1-βk with the condition σ=r+δβδ and investigate two different cases as below.

Case I: Based on the applications of the equations (2.25) for Lie point symmetry X˜2 , then the related null form

S˜2(t,p,k)=k(1+(-1)1+1σp-1σk-β(σ-1)1σ)βp, (2.45)
is obtained, indeed a particular solution of determining equation (2.13), and then the function K in the form (2.41) associated with the function S˜2 can be determined as
K˜21(t,p,k)=Λ2=e-t(r+δ-βδ)pkβ-1. (2.46)

The integrating factor R is

R˜21(t,p,k)=Λ1=e-t(r+δ-βδ)kββ(1-(-1)1σp-1σk-β(σ-1)1σ), (2.47)
where Λ1 and Λ2 are the adjoint symmetries of (2.11) and it is clear that they satisfy the adjoint invariance condition Λ1k = Λ2p. Additionally, the first integral of the system (2.11)
I˜21(t,p,k)=-e-t(r+δ-βδ)p1-1σ(p1σkβ(r+δ-βδ)-(-1)1σ(r+δ)(σ-1)1σ)β(r+δ-βδ), (2.48)
is obtained. As a result, by using first integral above, an analytical solution to the Ramsey neoclassical economic growth model (2.7) corresponding to the condition σ=r+δβδ can be determined as
p˜(t)=(c1et(β-1)(r+δ)βσ+(-1)β-1βσβ1+β(σ-1)βσδ1-ββσ(r+δ)-1/β(r+δ-βδ)-1+β+σ-1βσ)βσβ-1,k˜(t)=(-1)1βσβ-1βσδ-1βσ(r+δ)1/β(r+δ-βδ)-1β+1βδ×(c1et(β-1)(r+δ)βσ+(-1)β-1βσβ1+β(σ-1)βσδ1-ββσ(r+δ)-1/β(r+δ-βδ)-1+β+σ-1βσ)11-β,c˜(t)=(-1)1σ(σ-1)1σ(c1et(β-1)δ+(-1)(β-1)δr+δβr+2δ-βδr+δδδ-βδr+δ×(r+δ)-1/β(r+δ-βδ)(1-β)(r+δ-βδ)β(r+δ))β1-β, (2.49)
where c1 is an arbitrary constant.

Case II: Now let us consider the other solution in the following form for the function K

K˜22(t,p,k)=Λ2=et(1-σ)kβ-1(σ-1)(-1)p(σ-1)-1δ, (2.50)

which is the first integrating factor for the system. The second integrating factor R of (2.11) can be found as

R˜22(t,p,k)=Λ1=1-et(1-σ)kβ(σ-1)(-1)p(σ-1)-1(1+(-1)1+1σp-1σk-β(σ-1)1σ)δβp, (2.51)
where Λ1 and Λ2 functions are the adjoint symmetries of (2.11) and they clearly satisfy the adjoint invariance condition Λ1k = Λ2p. Thus, the first integral of system (2.11) is
I˜22(t,p,k)=et(1-σ)p-1σ(-1)p(σ-1)-1(p1σkβ(σ-1)+(-σ)1+1σ(r+δ)(σ-1)1σ)βδ, (2.52)

Another form of an analytical solution of the model (2.7) for the same condition σ=r+δβδ can be obtained by using the first integral (2.52) I22 as follows

p˜(t)=(c1et(β-1)(r+δ)βδ+(-1)β-1βσβ(r+δ)-1/β(σ-1)(β-1)(1-σ)βδ)βδβ-1,k˜(t)=(-1)1βσ(r+δ)1β(c1et(β-1)(r+δ)βδ+(-1)β-1βσβ(r+δ)-1/β(σ-1)(β-1)(1σ-1)β)11-β(σ-1)1-σβσ,c˜(t)=(-1)1σ(c1et(β-1)(r+δ)βδ+(-1)β-1βσβ(r+δ)-1/β(σ-1)(1-β)(σ-1)βσ)β1-β(σ-1)1σ, (2.53)
where c1 is an arbitrary constant. For the system (2.11), its λ-function and λ-symmetry are determined as below
λ˜2(t,p,k)=r+δ+(-1)1+1σp-1σ(β-1)(σ-1)1σ-kβk,V˜2=p-(k(1+(-1)1+1σp-1σk-β(σ-1)1σ)βp)k. (2.54)

In addition, Darboux polynomial for this case is

F˜1(t,p,k)=p(kβ(r+δ-βδ)+(-1)1+1σp-1σ(r+δ)(σ-1)1σ)etδ(β-1)(β-1)δ, (2.55)
which satisfies the condition (2.27).


This study introduces an analytical approach for analyzing Ramsey’s dynamical model defining a neoclassical growth model in the theory of economics via the theory of Lie groups. In the literature, there are different qualitative and quantitative methods to deal with the optimal control problems related to dynamic economic growth models in terms of Hamiltonians. However, the crucial problem for nonlinear economic growth models is the fact that there is no general analytical solution procedure for the analytical investigation for the related problems. On the other hand, with this study, we represent a unified approach including Lie point symmetry, Jacobi last multiplier (JLM), and Prelle-Singer (PS) methods for Ramsey economic growth model with current value Hamiltonian.

Firstly, we deal with the associated coupled-nonlinear first-order ordinary differential equations (ODEs) corresponding to first-order conditions of maximum principle for current value Hamiltonian form. From the Lie group point of view, it is a fact that the symmetry group analysis of first-order coupled nonlinear ODEs requires troublesome calculations since it is not possible to determine the Lie point symmetries by using the software programs known in the literature. From the theory of economics point of view, we consider the Ramsey economic growth model in terms of the current value Hamiltonian. The problem is that, in general, the Lie symmetries of the system (2.11) cannot be determined in closed form. Indeed, solving such problem is equivalent to the original problem. However, also in the particular cases the full symmetry group of (2.11) is infinite-dimensional. It is pointed out that Lie point symmetries for the Ramsey growth model do not present enough information in order to determine the first integrals and the exact closed-forms solution since the model has a restricted number of point symmetries. For this purpose, the relation between Lie point symmetries and JLM is considered. We also show that the JLM approach can be used as an effective method for analysis of the first-order nonlinear coupled system by considering its Lie point symmetries together.

For the case of current value Hamiltonian, we show that the system has a finite subalgebra such as X=t and prove that there is no analytical solution for the Ramsey model by using the PS method if there is no any restriction on the model parameters, which corresponds to the general Ramsey model. Additionally, based on the parametric restrictions σ = −1, β=12 , δ=-2r3 , it is shown that the system admits four-dimensional Lie algebra L4 and six different Jacobi last multipliers. The associated closed-form solutions for each subcase are found. As the other case, for the condition σ=r+δβδ , we present that the system admits two-dimensional finite Lie subalgebra L2 and only symmetry X˜2 gives a nontrivial closed-form solution by using PS method. In this case, two different analytical solutions to the Ramsey dynamical system are represented.

We also take into account the relation between Lie point symmetries and extended PS method to determine integrating factors and first integrals since the JLM approach does not always yield a first integral thus an analytical solution due to the symmetry restriction problem in which the first-order coupled and highly nonlinear ODEs may have. For the determination of the function K called the integrating factor, it is clear that the determination of the solution for the related system of equations is not straightforward and then we consider an ansatz form for the function K. However, it is shown that the compatibility equation (2.15) is not always satisfied by functions K and S, null function, obtained from (2.13) and (2.14), which then gives a trivial solution for function K. For the other cases, the associated first integrals and nontrivial closed-form solutions are obtained and they are represented for the different cases. Additionally, the related λ-symmetries, adjoint symmetries, Darboux polynomials, their properties, and relations for all cases are determined.

This study can be seen as an application of the theory of Lie groups to the Ramsey dynamical system. Based on the theoretical Lie group analysis, the novel closed-form solutions to the Ramsey economic growth model are reported for the first time in the literature.


The authors declare they have no conflicts of interest.


[6]G Darboux, Mémoire sur les équations différentielles algébriques du premier ordre et du premier degré, Bull. Sci. Math., Vol. 2, 1878, pp. 151-200.
[11]S Lie, Veralgemeinerung und neue Verwerthung der Jacobischen Multiplicator-Theorie, Christ. Forh., 1874, pp. 198-226.
[19]PJ Olver, Applications of Lie Groups to Differential Equations, Springer-Verlag, 1993.
[22]GG Polat and T Özer, On group analysis of optimal control problems in the theory of economic growth, Discrete Continuous Dyn. Syst. Ser., Vol. S 13, 2020, pp. 2853-2876.
Journal of Nonlinear Mathematical Physics
Publication Date
ISSN (Online)
ISSN (Print)
DOI to use a DOI?
© 2021 The Authors. Published by Atlantis Press B.V.
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (

Cite this article

AU  - Gülden Gün Polat
AU  - Teoman Özer
PY  - 2021
DA  - 2021/01
TI  - On Ramsey Dynamical Model and Closed-Form Solutions
JO  - Journal of Nonlinear Mathematical Physics
SN  - 1776-0852
UR  -
DO  -
ID  - Polat2021
ER  -