An important problem in mathematics of finance is the mathematical modelling of optimal investment-consumption decisions under uncertainty conditions. The investment-consumption problem  has been extensively investigated in many works with various modifications and extensions. Cox  and Cox and Ross  have derived the well-known constant elasticity of variance (CEV) option pricing model and Schroder  afterwards broadened the model by stating the CEV option pricing formula in respect of noncentral Chi-square distribution. The CEV model is mostly used to investigate the option and asset pricing formula, as was investigated by Beckers , Davydov and Linetsky , Emanuel and Macbeth  and recently by Hsu et al. . Here, we reconsider the CEV model . Literature survey witness that there have been few studies [10,11,12,13,14,15,16,17] recently reported where its solutions are presented. Nonetheless, to the best of our knowledge, no work has been published on the closed-form solutions of the CEV model, which is the aim of the present work.
The classical Lie symmetry theory was discovered by the Norwegian mathematician Marius Sophus Lie (1842-1899) in the nineteenth century. This theory systematically unites the widely known ad hoc methods to find exact solutions for differential equations. After many years of discovery, the Lie’s theory was popularized by Ovsiannikov in Novosibirsk, Russia in the middle of twentieth century and by Birkhoff and Olver in the West. Lie’s theory is one of the most effective tools to find exact analytical solutions of nonlinear partial differential equations and is established on the analysis of the invariance under one-parameter group of point transformations. See for example, the references [18,19,20,21,22,23,24].
Recently Lie’s theory has been applied to partial differential equations (PDEs) of mathematical finance. One of the earliest study on the subject is , where the classical Black-Scholes equation was discussed. The bond-pricing equation via Lie group approach was investigated in [26, 27]. The invariant analysis of some well-known financial mathematics models was presented [28,29,30,31] . More recently the Lie’s theory has been applied to various PDEs of financial mathematics. See for example, Motsepa et al. , Lekalakala et al. , Nteumagne and Moitsheki , Bozhkov and Dimas , Caister et al. [36, 37], Naicker et al. , Lo , Taylor and Glasgow , Wang et al.  and Pooe et al. , are few important studies to mention.
In this paper, we discuss the optimal investment-consumption problem under the CEV model
with the terminal condition
from the viewpoint of Lie symmetry approach. The application of Lie’s theory, in general, reduces the partial differential equation in two independent variables to an ordinary differential equation and provides us with the group-invariant solution. Some nontrivial conservation laws are also constructed for the CEV model by employing a general theorem proved in .
2 Lie symmetries of (1)
In this section, we compute the Lie symmetries admitted by the model equation (1) and utilised them to obtain the closed-form group-invariant solution of the PDE (1) satisfying the terminal condition (2).
For detailed description on Lie symmetry method and its applications to various disciplines the reader is referred to Refs [18,19,20,21,22,23,24]. However, in this section we give detailed calculations on finding Lie point symmetries of (1).
In order to facilitate the calculations when employing the Lie group approach, we rewrite the PDE (1) in the form
is a Lie symmetry of PDE (3) if and only if
where X denotes the second prolongation of X which is defined by
and the total derivative operators
which are linear PDEs in τ, ξ and η. Solving the above system gives
where and C1, …, C4 are arbitrary constants. Since τ, ξ and η contain four constants, we conclude that the Lie algebra of infinitesimal symmetries of the PDE (3) is spanned by the four vector fields
3 Group invariant solution
Many researchers have developed various analytical methods for solving partial differential equations, such as inverse scattering transform , the Bäcklund transformation , the Hirota bilinear method , the Painlevé analysis method , the Bell polynomials , the homoclinic breather limit method .
We now obtain the closed-form group-invariant solution for the PDE (3) by making use of the Lie symmetry algebra calculated in the previous section. Firstly we calculate symmetry Lie algebra admitted by (3) that satisfies the terminal condition (2) [26, 28].
We consider the linear combination of the Lie point symmetries, namely
and use the terminal conditions
The condition (9) yields
while condition (10) gives
Splitting Eq. (12) on powers of x yields
Thus, the group-invariant solution of the PDE (3) is given by J2 = G(J1), which yields
Solving the reduced ODE (17) for G(z), we obtain
Substituting the value of G into Eq. (16), the solution for F(x, t) is written as
Finally, making use of the terminal condition F(x, t) = 1 in Eq. (18), we obtain
Therefore, the solution of (1+1) evolution PDE (3) satisfying the terminal condition is given by
4 Conservation laws
In this section we derive conservation laws for the (1+1) evolution partial differential equation (3). In classical physics, conservations laws are physical quantities which describe the conservation of energy, mass, linear momentum, angular momentum, and electric charge. One particularly important result concerning conservation laws is the celebrated Noether theorem, which gives us a sophisticated and useful way of constructing conservation laws when a Noether point symmetry connected to a Lagrangian is known for the corresponding Euler-Lagrange equation.
Recently, a theorem due to Ibragimov was proved, which gives a method to construct conservation laws irrespective of the existence of a Lagrangian. We start here by stating the definition of an adjoint equation.
be a second-order PDE with x, t as independent variables and f a dependent variable. Then, its adjoint equation is 
denoting the Euler–Lagrange operator and g is a new dependent variable. Here
are the total derivative operators.
The 2-tuple vector T = (Tt, Tx), is a conserved vector of (19) if Tt and Tx satisfy
We now state the following theorem.
 Any Lie point, Lie–Bäcklund or non-local symmetry
where w and 𝓛 are determined as follows:
We now apply the above theorem to our problem. The adjoint equation to the equation (3) is
We now apply the above theorem to each Lie point symmetry of Eq. (3). We start with X1 = ∂/∂t. Corresponding to symmetry X1 we obtain the conserved vector with components
The symmetry X2 provides us with the conserved vector whose components are
where . Likewise, the conserved vectors associated with the symmetries X3 and X4 are given by
5 Concluding remarks
The evolution (1+1) PDE (1) describing the optimal investment-consumption problem under the CEV model  satisfied the classical Black–Scholes–Merton equation with boundary condition which differ from those often used in the most common cases. It is well-known that the evolution (1+1) PDE (1) is related to the heat equation via the equivalence transformations and thus its general solution can be obtained. However, in this paper, for the first time, we have solved the PDE (1) subject to the terminal condition (2) by utilizing the Lie group method. This demonstrates the usefulness of Lie’s theory. We found a four-dimensional Lie symmetry algebra for evolution PDE (1). Using the nontrivial Lie point symmetry operator, we have shown that the governing PDE can be transformed into a second-order variable coefficient ODE. The reduced ODE is solved to obtain a new exact closed-form solution of the CEV model which also satisfy the terminal condition. Thus for the first time with the application of Lie’s theory closed-form solution of (1) is derived. Finally, we constructed conservation laws corresponding to the four Lie point symmetries by employing a general theorem on conservation laws. This is the first time that the evolution PDE (1) for optimal investment-consumption problem has been considered from the view point of group theoretical approach and the conservation laws have been derived in the literature.
Cox J.C., The constant elasticity of variance option pricing model, J. Portfolio Manage., 1996, 22, 16-17. Google Scholar
Emanuel D., Macbeth J., Further results on the constant elasticity of variance call option pricing model, J. Financial Quan. Anal., 1982, 17, 53-54. Google Scholar
Chang H., Rong X.M., Zhao H., Zhang C.B., Optimal investment and consumption decisions under the constant elasticity of variance model, Math. Prob. Eng., 2013, 974098 Google Scholar
Xiao J., Zhai H., Qin C., The constant elasticity of variance (CEV) model and the Legendre transform-dual solution for annuity contracts, Insurance: Mathematics and Eco., 2007, 40, 302 - 310. Google Scholar
Gao J., An extended CEV model and the Legendre transform-dual-asymptotic solutions for annuity contracts, Insurance: Mathematics and Eco., 2010, 46, 511-530. Google Scholar
Ballestra L.V., Cecere L., Pricing American options under the constant elasticity of variance model: An extension of the method by Barone-Adesi and Whaley, Finance Research Lett., 2015, 14, 45-55. CrossrefGoogle Scholar
Rong X., Zhao H., Time-consistent reinsurance-investment strategy for an insurer and a reinsurer with mean-variance criterion under the CEV model, J. Comput. Appl. Math., 2015, 283, 142-162. CrossrefGoogle Scholar
Ovsiannikov L.V., Group Analysis of Differential Equations, Academic Press, New York, 1982. Google Scholar
Olver P.J., Applications of Lie Groups to Differential Equations, Springer, New York, 1986. Google Scholar
Bluman G.W., Kumei S., Symmetries and Differential Equations, Springer, New York,1989. Google Scholar
Ibragimov N.H., Elementary Lie Group Analysis and Ordinary Differential Equations, John Wiley & Sons, Chichester, 1999. Google Scholar
Ma P.L., Tian S.F., Zhang T.T., On symmetry-preserving difference scheme to a generalized Benjamin equation and third order Burgers equation, Appl. Math. Lett., 2015, 50, 146-152. CrossrefGoogle Scholar
Tu J.M., Tian S.F., Xu M.J., Zhang T.T., On Lie symmetries, optimal systems and explicit solutions to the Kudryashov-Sinelshchikov equation, Appl. Math. Comp., 2016, 275, 345-352 CrossrefGoogle Scholar
Sinkala W., Leach P.G.L., O’Hara J.G., Zero-coupon bond prices in the Vasicek and CIR models: Their computation as group invariant solutions, Math. Meth. Appl. Sci., 2008, 31, 665-678. CrossrefGoogle Scholar
Sinkala W., Leach P.G.L., O’Hara J.G., An optimal system and group-invariant solutions of the Cox-Ingersoll-Ross pricing equation, Appl. Math. Comput., 2008, 201, 95-107. Google Scholar
Sinkala W., Leach P.G.L., O’Hara J.G., Embedding the Vasicek model into the Cox-Ingersoll-Ross model, Math. Meth. Appl. Sci., 2001, 34, 152-159. Google Scholar
Motsepa T., Khalique C.M., Molati M., Group classification of a general bond-option pricing equation of mathematical finance, Abstract Appl. Anal. Volume 2014, Article ID 709871, 10 pages. Google Scholar
Bozhkov Y., Dimas S., Group classification of a generalization of the Heath equation, Appl.Math. Comput., 2014, 243, 121-131. Google Scholar
Caister N.C., Govinder K.S., O’Hara J.G., Solving a nonlinear pde that prices real options using utility based pricing methods, Nonlinear Anal.: Real World Appl., 2011, 12, 2408-2415. CrossrefGoogle Scholar
Naicker V., O’Hara J.G., Leach P.G.L., A note on the integrability of the classical portfolio selection model, Appl. Math. Lett., 2010, 23, 1114-1119. Google Scholar
Lo C.F., Lie algebraic approach for pricing zero-coupon bonds in single-factor interest rate models, J. Appl. Math., Volume(2013), Article ID 276238, 9 pages. Google Scholar
Wang Z., Wang L., Wang D.S., Jin Y., Optimal system, symmetry reductions and new closed form solutions for the geomet- ric average Asian options. Appl. Math. Comput., 2014, 226, 598-605. Google Scholar
Ablowitz M.J., Clarkson P.A., Solitons; Nonlinear Evolution Equations and Inverse Scattering, Cambridge Univ. Press, 1991. Google Scholar
Xu M.J., Tian S.F., Tu J.M., Zhang T.T., Bäcklund transformation, infinite conservation laws and periodic wave solutions to a generalized (2+1)-dimensional Boussinesq equation, Nonlinear Analysis: Real World Appl,.2016, 31, 388-408. CrossrefGoogle Scholar
Hirota R., Direct Methods in Soliton Theory, Springer, 2004. Google Scholar
Tu J.M., Tian S.F., Xu M.J., Ma P.L., Zhang T.T., On periodic wave solutions with asymptotic behaviors to a (3+1)-dimensional generalized B-type Kadomtsev-Petviashvili equation in fluid dynamics, Comput. Math App., 2016, 72, 2486-2504. Google Scholar
Wang X.B., Tian S.F., Qin C.Y., Zhang T.T., Dynamics of the breathers, rogue waves and solitary waves in the (2+1)- dimensional Ito equation, Appl. Math. Lett., 2017, 68, 40-47. CrossrefGoogle Scholar
Wolfram Research, Inc., Mathematica, Version 10.0, Champaign, IL (2014). Google Scholar
Abramowitz M., Stegun I.A., Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 9th printing. New York: Dover, 1972, 260 Google Scholar
About the article
Published Online: 2018-03-08
Conflict of interestConflict of Interests: The authors declare that there is no conflict of interests regarding the publication of this paper.