Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Physics

formerly Central European Journal of Physics

Editor-in-Chief: Seidel, Sally

Managing Editor: Lesna-Szreter, Paulina

1 Issue per year

IMPACT FACTOR 2016 (Open Physics): 0.745
IMPACT FACTOR 2016 (Central European Journal of Physics): 0.765

CiteScore 2016: 0.82

SCImago Journal Rank (SJR) 2015: 0.458
Source Normalized Impact per Paper (SNIP) 2015: 1.142

Open Access
See all formats and pricing
More options …
Volume 16, Issue 1


Volume 13 (2015)

Algebraic aspects of evolution partial differential equation arising in the study of constant elasticity of variance model from financial mathematics

Tanki Motsepa
  • International Institute for Symmetry Analysis and Mathematical Modeling, Department of Mathematical Sciences, North-West University, Mafikeng Campus, Private Bag X 2046, Mmabatho, 2735, South Africa
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Taha Aziz
  • Corresponding author
  • School of Computer, Statistical and Mathematical Sciences, North-West University, Potchefstroom Campus, Private Bag X 6001, Potchefstroom, 2531, South Africa
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Aeeman Fatima
  • International Institute for Symmetry Analysis and Mathematical Modeling, Department of Mathematical Sciences, North-West University, Mafikeng Campus, Private Bag X 2046, Mmabatho, 2735, South Africa
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Chaudry Masood Khalique
  • International Institute for Symmetry Analysis and Mathematical Modeling, Department of Mathematical Sciences, North-West University, Mafikeng Campus, Private Bag X 2046, Mmabatho, 2735, South Africa
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
Published Online: 2018-03-08 | DOI: https://doi.org/10.1515/phys-2018-0006


The optimal investment-consumption problem under the constant elasticity of variance (CEV) model is investigated from the perspective of Lie group analysis. The Lie symmetry group of the evolution partial differential equation describing the CEV model is derived. The Lie point symmetries are then used to obtain an exact solution of the governing model satisfying a standard terminal condition. Finally, we construct conservation laws of the underlying equation using the general theorem on conservation laws.

Keywords: CEV model; Terminal condition; Lie point symmetries; Group-invariant solution; Conservation laws

PACS: 02.20.Sv; 02.20.Tw; 02.30.Jr; 02.60.Lj

1 Introduction

An important problem in mathematics of finance is the mathematical modelling of optimal investment-consumption decisions under uncertainty conditions. The investment-consumption problem [1] has been extensively investigated in many works with various modifications and extensions. Cox [2] and Cox and Ross [3] have derived the well-known constant elasticity of variance (CEV) option pricing model and Schroder [4] 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 [5], Davydov and Linetsky [6], Emanuel and Macbeth [7] and recently by Hsu et al. [8]. Here, we reconsider the CEV model [9]. 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 [25], 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. [32], Lekalakala et al. [33], Nteumagne and Moitsheki [34], Bozhkov and Dimas [35], Caister et al. [36, 37], Naicker et al. [38], Lo [39], Taylor and Glasgow [40], Wang et al. [41] and Pooe et al. [42], 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 [43].

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




The operator


is a Lie symmetry of PDE (3) if and only if


where X[2] denotes the second prolongation of X which is defined by




and the total derivative operators


Expanding the symmetry condition (4) with the help of (5)-(7) and then splitting on the derivatives of F leads to


which are linear PDEs in τ, ξ and η. Solving the above system gives


where θ=c24be 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 [44], the Bäcklund transformation [45], the Hirota bilinear method [46], the Painlevé analysis method [47], the Bell polynomials [48], the homoclinic breather limit method [49].

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



Solving equations (11), (13) and (14), we obtain


Substituting the values of a1, a3 and a4 from Eq. (15) into Eq. (8) and solving X(F) = 0, we obtain the two invariants


Thus, the group-invariant solution of the PDE (3) is given by J2 = G(J1), which yields


where z = J1. Now substituting Eq. (16) into PDE (3) yields a second-order ODE for G(z), viz.,


Solving the reduced ODE (17) for G(z), we obtain


where C1 and C2 are constants and Γ [.,.] is the incomplete Gamma function [50, 51]. The terminal condition dictates that we take C2 = 0, hence the above solution for G(z) takes the form


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.

Definition 1



be a second-order PDE with x, t as independent variables and f a dependent variable. Then, its adjoint equation is [43]






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.

Theorem 1

[43] Any Lie point, Lie–Bäcklund or non-local symmetry


of Eq. (19) gives a conservation law Di(Ti) = 0 for the system (19)(20). The conserved vector is given by


where w and 𝓛 are determined as follows:


We now apply the above theorem to our problem. The adjoint equation to the equation (3) is


and hence the second-order Lagrangian of the equation (3) and its adjoint equation (22) is written as


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 θ=c24be.. 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 [9] 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.


  • [1]

    Merton R.C., Optimum consumption and portfolio rules in a continuous time model, J. Economic Theory, 1971, 3, 373-413. CrossrefGoogle Scholar

  • [2]

    Cox J.C., The constant elasticity of variance option pricing model, J. Portfolio Manage., 1996, 22, 16-17. Google Scholar

  • [3]

    Cox J.C, Ross S.A., The valuation of options for alternative stochastic processes, J. Financial Economics, 1976, 3, 145-166. CrossrefGoogle Scholar

  • [4]

    Schroder M., Computing the constant elasticity of variance option pricing formula, J. Finan., 1989, 44, 211-219. CrossrefGoogle Scholar

  • [5]

    Beckers S., The constant elasticity of variance model and its implications for option pricing, J. Finan., 1980, 35, 661-673. CrossrefGoogle Scholar

  • [6]

    Davydov D., Linetsky V., The valuation and hedging of barrier and lookback option under the CEV process, Management Sci., 2001, 47, 949-965. CrossrefGoogle Scholar

  • [7]

    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

  • [8]

    Hsu Y.L., Lin T.I., Lee C.F., Constant elasticity of variance (CEV) option pricing model: Integration and detailed derivation, Math. Comput. Simul., 2008, 79, 60-71. CrossrefGoogle Scholar

  • [9]

    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

  • [10]

    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

  • [11]

    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

  • [12]

    Lo C.F., Yuen P.H., Hui C.H., Constant elasticity of variance option pricing model with time-dependent parameters, Int. J. Theor. Appl. Fin., 2000, 3, 661-674. CrossrefGoogle Scholar

  • [13]

    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

  • [14]

    Ballestra L.V., Pacelli G., The constant elasticity of variance model: calibration, test and evidence from the Italian equity market, Appl. Financial Eco., 2011, 21, 1479-1487. CrossrefGoogle Scholar

  • [15]

    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

  • [16]

    Tsai Wei-Che, Improved method for static replication under the CEV model, Finance Research Lett., 2014, 11, 194-202. CrossrefGoogle Scholar

  • [17]

    Zhao H., Rong X., On the constant elasticity of variance model for the utility maximization problem with multiple risky assets, IMA J. Management Math., 2017, 28, 299-320. CrossrefGoogle Scholar

  • [18]

    Ovsiannikov L.V., Group Analysis of Differential Equations, Academic Press, New York, 1982. Google Scholar

  • [19]

    Olver P.J., Applications of Lie Groups to Differential Equations, Springer, New York, 1986. Google Scholar

  • [20]

    Bluman G.W., Kumei S., Symmetries and Differential Equations, Springer, New York,1989. Google Scholar

  • [21]

    Ibragimov N.H., Elementary Lie Group Analysis and Ordinary Differential Equations, John Wiley & Sons, Chichester, 1999. Google Scholar

  • [22]

    Motsepa T., Khalique C.M., Conservation laws and solutions of a generalized coupled (2+1)-dimensional burgers system, Comput. Math. Appl., 2017, 74, 1333-1339. CrossrefGoogle Scholar

  • [23]

    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

  • [24]

    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

  • [25]

    Gazizov R.K., Ibragimov N.H., Lie symmetry analysis of differential equations in finance, Nonlinear Dynam., 1998, 17, 387-407. CrossrefGoogle Scholar

  • [26]

    Goard J., New solutions to the bond pricing equation via Lie’s classical method, Math. Comput. Model., 2000, 32, 299-313. CrossrefGoogle Scholar

  • [27]

    Goard J., P. Broadbridge and G. Raina, Tractable forms of the bond pricing equation, Math. Comput. Model., 2004, 40, 151-172. CrossrefGoogle Scholar

  • [28]

    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

  • [29]

    Sinkala W., Leach P.G.L., O’Hara J.G., Invariance properties of a general bond-pricing equation, J. Differential Equ., 2008, 244, 2820-2835. CrossrefGoogle Scholar

  • [30]

    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

  • [31]

    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

  • [32]

    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

  • [33]

    Lekalakala S.L., Motsepa T., Khalique C.M., Lie Symmetry Reductions and Exact Solutions of an Option-Pricing Equation for Large Agents, Mediterr. J. Math., 2016, 13, . CrossrefGoogle Scholar

  • [34]

    Nteumagne B.F., Moitsheki R.J., Optimal systems and group in variant solutions for a model arising in financial mathematics, Math. Model. Anal., 2009, 14, 495-502. CrossrefGoogle Scholar

  • [35]

    Bozhkov Y., Dimas S., Group classification of a generalization of the Heath equation, Appl.Math. Comput., 2014, 243, 121-131. Google Scholar

  • [36]

    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

  • [37]

    Caister N.C., Govinder K.S., O’Hara J.G., Optimal system of Lie group invariant solutions for the Asian option PDE, Math. Meth. Appl. Sci., 2011, 34, 1353-1365. CrossrefGoogle Scholar

  • [38]

    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

  • [39]

    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

  • [40]

    Taylor S.M., Glasgow S.A., A novel reduction of the simple Asian option and Lie-group invariant solutions, Int. J. Theoret. Appl. Fin., 2009, 12, 1197. CrossrefGoogle Scholar

  • [41]

    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

  • [42]

    Pooe C.A., Mahomed F.M., Wafo-Soh C., Fundamental solutions for zero-coupon bond pricing models, Nonlinear Dyn., 2004, 36, 69-76. CrossrefGoogle Scholar

  • [43]

    Ibragimov N.H., A new conservation theorem, J. Math. Anal. Appl., 2007, 333, 311-328. CrossrefGoogle Scholar

  • [44]

    Ablowitz M.J., Clarkson P.A., Solitons; Nonlinear Evolution Equations and Inverse Scattering, Cambridge Univ. Press, 1991. Google Scholar

  • [45]

    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

  • [46]

    Hirota R., Direct Methods in Soliton Theory, Springer, 2004. Google Scholar

  • [47]

    Weiss J., Tabor M., Carnevale G., The Painlevé property for partial differential equations, J. Math. Phys., 1983, 24, 522-526. CrossrefGoogle Scholar

  • [48]

    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

  • [49]

    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

  • [50]

    Wolfram Research, Inc., Mathematica, Version 10.0, Champaign, IL (2014). Google Scholar

  • [51]

    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

Tel.: +27-18-2852494

Received: 2017-08-07

Accepted: 2017-10-28

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.

Citation Information: Open Physics, Volume 16, Issue 1, Pages 31–36, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2018-0006.

Export Citation

© 2018 T. Motsepa et al.. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0

Citing Articles

Here you can find all Crossref-listed publications in which this article is cited. If you would like to receive automatic email messages as soon as this article is cited in other publications, simply activate the “Citation Alert” on the top of this page.

Chaudry Masood Khalique and Tanki Motsepa
Physica A: Statistical Mechanics and its Applications, 2018

Comments (0)

Please log in or register to comment.
Log in