Next Article in Journal
Singularity Formation in the Inviscid Burgers Equation
Next Article in Special Issue
Boundary Value Problems of Hadamard Fractional Differential Equations of Variable Order
Previous Article in Journal
Constraints on Newtonian Interplanetary Point-Mass Interactions in Multicomponent Systems from the Symmetry of Their Cycles
Previous Article in Special Issue
New Hermite–Hadamard Inequalities in Fuzzy-Interval Fractional Calculus and Related Inequalities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Invariant Solutions of Black–Scholes Equation with Ornstein–Uhlenbeck Process

by
Maba Boniface Matadi
*,† and
Phumlani Lawrence Zondi
Department of Mathematical Sciences, Faculty of Science, Agriculture and Engineering, University of Zululand, Private Bag X1001, KwaDlangezwa 3886, South Africa
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2021, 13(5), 847; https://doi.org/10.3390/sym13050847
Submission received: 3 February 2021 / Revised: 23 March 2021 / Accepted: 24 March 2021 / Published: 11 May 2021
(This article belongs to the Special Issue Functional Equations and Analytic Inequalities)

Abstract

:
This paper analyses the model of Black–Scholes option pricing from the point of view of the group theoretic approach. The study identified new independent variables that lead to the transformation of the Black–Scholes equation. Furthermore, corresponding determining equations were constructed and new symmetries were found. As a result, the findings of the study demonstrate of the integrability of the model to present an invariant solution for the Ornstein–Uhlenbeck stochastic process.

1. Introduction

In 1973, Fischer Black and Myron Scholes formulated a mathematical model that describes option pricing. Their model forms the cornerstone for modern financial theory; it is not often that one can talk about modern finance without ever mentioning the revolutionary Black–Scholes (BS) model. The central goal of the BS model is to find fair prices of options by combining the price variation of the stock, the time value of money, the option’s stock price, and the time to the option’s expiry [1]. In 1976, Robert Merton provided several extensions to the BS model, the most prominent of these extensions being the ability to account for dividend yields [2]. The importance of the BS model was recognised in 1997 by Royal Swedish Academy of Sciences, where Scholes and Merton were awarded the Nobel Prize in Economics. Unfortunately, having died two years earlier, Fischer Black was not present for the award.
In this paper, we demonstrate how the theory of Lie symmetry analysis can be used to find invariant solutions of the Black–Scholes option pricing model under stochastic volatility. The stochastic volatility case was handled using the method introduced by Palianthanasis et al. [3], where the volatility was assumed to follow the Ornstein–Uhlenbeck process. The model is governed by the given partial differential equation [3]:
V t + 1 2 f 2 ( y ) S 2 V S S + ρ β S f ( y ) V S y + 1 2 β 2 V y y + r S V S + [ α ( m y ) β ρ μ r f ( y ) ] V y r V = 0 ,
where S is the price of an underlying asset (the stock); ρ is the correlation coefficient; β is the drift; r is the riskless rate; α is the rate of mean reversion; m is the long-run mean of Y; σ is the volatility of the stock. To derive the model for stochastic volatility with the Ornstein–Uhlenbeck process, it is convenient to regard the volatility of the stock price as being a function of Y t , i.e., σ = f ( Y t ) , where Y t is the stochastic process that contains the Ornstein–Uhlenbeck term [3].
This paper is organised as follows. In Section 3, a heuristic background of the concepts underlying Lie symmetry analysis are introduced. In Section 4, Lie symmetry analysis is presented for the case where the volatility follows a stochastic process. The obtained symmetries were used to reduce Equation (1) into a one dimensional linear ordinary differential equation of order two. In Section 5, the modified local one-parameter transformations are presented and used to calculate invariant solutions of Equation (1). Furthermore, Lie symmetry analysis was used to examine a stochastic Heston model. The numerical solution that represents the evolution of the solutions is presented graphically in Section 6. The conclusion is provided in Section 7.

2. Fundamental Definitions and Theorems

In this section, a comprehensive review of a group theoretic approach to the solution of differential equations is given. The theory entails the tools necessary for subsequent employment throughout the paper. To start with, the mathematical idea of a symmetry is explained, and then the general properties of groups are explained; the properties are then extended to the Lie groups. Numerous textbooks are accessible [4,5,6,7,8], as well as research papers extensively published on the symmetry analysis of ordinary differential equations (ODEs) and partial differential equations (PDEs). In [5], Bluman and Kumei described the significance of the theory of Lie symmetry for a PDE. The authors pointed out that invariant functions can be constructed to reduce the order of a differential equation or the number of dependent variables [9].
According to Lie theory, the kth-order partial differential equation [10,11]:
u t F ( t , x , u , u ( 1 ) , u ( 2 ) , , u ( k ) ) = 0 ,
admits the given Lie group of transformations of one-parameter:
t ^ t + a ξ 0 ( t , x , u , u ( 1 ) , u ( 2 ) , , u ( k ) ) , x ^ i x i + a ξ i ( t , x , u , u ( 1 ) , u ( 2 ) , , u ( k ) ) , u ^ i α u i α + a η i α ( t , x , u , u ( 1 ) , u ( 2 ) , , u ( k ) ) ,
with infinitesimal Lie generators [12,13]:
X = ξ 0 t + ξ i x i + η α u α ,
if:
u ^ t F ( t ^ , x ^ , u ^ , u ^ ( 1 ) , u ^ ( 2 ) , , u ^ ( k ) ) = 0 .
the group transformations t ^ , x ^ , and u ^ are obtained by solving the following Lie equations [14]:
d t ^ d a = ξ 0 ( t ^ , x ^ , u ^ , u ^ ( 1 ) , u ^ ( 2 ) , , u ^ ( k ) ) , d x ^ i d a = ξ i ( t ^ , x ^ , u ^ , u ^ ( 1 ) , u ^ ( 2 ) , , u ^ ( k ) ) , d u ^ i α d a = η i α ( t ^ , x ^ , u ^ , u ^ ( 1 ) , u ^ ( 2 ) , , u ^ ( k ) ) ,
with the initial conditions:
t ^ a = 0 = t , x ^ i a = 0 = x i , u ^ i α a = 0 = u i α .
the infinitesimal form of u ^ t ¯ , u ^ ( 1 ) , u ^ ( 2 ) , . . . , u ^ ( k ) are found by the given formulas [15,16]:
u ^ i α u i α + a η i α ( x , u , u 1 ) , u ^ i j α u i j α + a η i j α ( x , u , u 1 , u 2 ) , u ^ i 1 i k α u i 1 i k α + a η i 1 i k α ( x , u , u 1 , , u k ) .
the functions η i α ( x , u , u 1 ) , η i j α ( x , u , u 1 , u 2 ) , and η i 1 i k α ( x , u , u 1 , . . . , u k ) are obtained from the following prolongation formulas [6]:
η i α = D i ( η α ) u j α D i ( ξ j ) , η i j α = D j ( η i α ) u i l α D j ( ξ l ) , η i 1 i k α = D i k ( η i 1 . . . i k 1 α ) u i 1 i k 1 l α D i k ( ξ l ) ,
where D i denotes the operator of total differentiation with respect to ( x 1 , x 2 . . . x n ) , then:
D i = x i + u i α u α + u i j α u j α .
the transformed derivatives u ^ ( 1 ) , u ^ ( 2 ) , . . . , u ^ ( k ) can be computed from the formulas:
D i = D i ( f i ) D ^ j .
the generators are therefore given by:
X [ 1 ] = X + η i α ( x , u , u 1 ) u i α , X [ k ] = X [ 1 ] + + η i 1 i k α ( x , u , , u k ) u i 1 . . . i k α .
Theorem 1.
A function F ( x , u , , u k ) is invariant under the prolonged group G, if and only if [17]
X [ k ] F = 0 ,
where X [ k ] is the generator of G.
Theorem 2.
Every one-parameter group of transformations ( x ^ = f ( x , y , ε ) , y ^ = g ( x , y , ε ) ) is reduced to a group of translations t ^ = t + ε , u ^ = u with the generator [17]:
X = t ,
by suitable change of variables
t = t ( x , y ) , u = u ( x , y ) .
Considering the Lie groups of point transformations related to a given differential equation E involve n independent variables x = ( x 1 , x 2 , . . . , x n ) R n and m dependent variables u = ( u 1 , u 2 , . . . , u m ) R m [4,18], let:
x * = X ( x , u ; a ) , u * = U ( x , u ; a )
be a group of transformations in the space R n + m of the variables ( x , u ) [4,18]. Furthermore, for:
u = Θ ( x ) ( Θ 1 ( x ) , Θ 2 ( x ) , . . . , Θ m ( x ) ) ,
be a solution of the equation E . A Lie group of transformations of the form (12) admitted by E has the following two corresponding properties [4,18]:
1
A transformation of the group maps any solution of E into another solution of E ;
2
A transformation of the group leaves E invariant, say, E reads the same in terms of the variables ( x , u ) and in terms of the transformed variables ( x * , u * ) .
Definition 1.
[18]-The function u = Θ ( x ) with components u A = Θ A ( x ) ( A = 1 , 2 , , m ) is said to be an invariant solution of (2) if u A = Θ A ( x ) is an invariant surface of (12) and is a solution of (2).

3. Symmetry Analysis for Black–Scholes with the Ornstein–Uhlenbeck Process

We assume that Equation (1) admits a Lie symmetry satisfying the generator:
X [ 1 ] = ξ 1 t + ξ 2 S + ξ 3 y + η V .
the prolonged symbol is then:
X [ 2 ] = X [ 1 ] + η t V t + η S V S + η y V y + η S y V S y + η S S V S S + η y y V y y .
using Equation (7), we can obtain the following exact expressions for the prolonged coefficients:
η t = η t + V t η V ξ t 1 V t V t V t ξ V 1 ξ t 2 V S V t V S ξ V 2 ξ t 3 V y V t V y ξ V 3 η S = η S + V S η V ξ S 1 V t V t V S ξ V 1 ξ S 2 V S V S V S ξ V 2 ξ S 3 V y V S V y ξ V 3 η y = η y + V y η V ξ y 1 V t ξ V 2 V y V S V t V y ξ V 1 ξ y 2 V S ξ y 3 V y y V y V y ξ V 3 η S S = η S S + 2 η S V V S + η V V S S + η V V V S V S ξ S S 2 V S ξ S S 3 V y ξ S S 1 V t 2 ξ S 2 V S S 2 ξ S 3 V S y ξ S 1 V S t 2 ξ S V 2 V S V S 2 ξ S V 3 V S V y 2 ξ S V 1 V S V t ξ V 2 V S V S S 2 ξ V 2 V S V S S ξ V 3 V y V S S 2 ξ V 3 V S y ξ V 1 V t V S S 2 ξ V 1 V S V S t ξ V V 2 V S V S V S ξ V V 3 V S V S V y ξ V V 1 V S V S V t η y y = η y y + 2 η y V V y + η V V y y + η V V V y V y ξ y y 2 V S ξ y y 3 V y ξ y y 1 V t 2 ξ y 2 V y S 2 ξ y 3 V y y ξ y 1 V y t 2 ξ y V 2 V y V S 2 ξ y V 3 V y V y 2 ξ y V V y V t ξ V 2 V S V y y 2 ξ V 2 V y V y S ξ V 3 V y V y y ξ V V 2 V y V y V S ξ V V 3 V y V y V y ξ V V 1 v y V y V t η S y = η S y + η y V V S ξ S y 2 V S + η S V V y ξ S y 3 V y ξ S y 1 V t + η V V V S V y ξ y V 3 V S V y ξ S V 2 V S V y ξ y V 1 V S V t ξ y V 1 V y V t ξ y V 1 V S V S ξ S V 3 V y V y ξ y 2 V S S ξ S 3 V y y + η V V S y ξ y 3 V S y ξ S 2 V S y ξ y 1 V S t ξ S 1 V y t ξ V 2 V S V S y ξ V 1 V S S V y ξ V 3 V S V y y ξ V V 3 V S V y V y .
the symbols η t , η S , η y , η S y , η S S , and η y y are substituted into the determining equation. Using the SYM Mathematica package conceived by Dimas and Tsoubelis in [19], we obtained the following two Lie operators, which are inline in [20]:
X 1 = t ,
X 2 = S S .
in as much as Equation (1) is a linear equation, it therefore constantly admits the respective linear and infinite symmetries given by:
X V = V V ,
X b = b ( t , S , y ) V ,
in constructing invariant solutions, the following three facts must be taken into consideration: First, the infinite symmetry cannot be used for the reduction of the differential equation; therefore, they are not considered here. Second, all those solutions in which V does not depend on one of the independent variables are discarded. Third, linear combinations of symmetries are also symmetries. These three facts mean that a reduction of Equation (1) can be performed using the following symmetry generators, which are linear combinations of the two symmetries obtained:
Y 1 = X 1 + k 1 X u ,
Y 2 = X 2 + k 2 X u ,
Y 12 = X 1 + c X 2 + k 3 X u .

4. Invariant Solution through Lie Operators

4.1. Invariant Solution through Symmetry Y 1

The symmetry Y 1 is given by:
Y 1 = t + k 1 V V .
the characteristic system is therefore given by:
k 1 d t 1 = d V V
integrating both sides, we can obtain:
k 1 t = ln V + ln ν ( S , y ) ,
which results in:
V ( t , S , y ) = ν ( S , y ) e k 1 t .
the computation of the partial derivatives of V gives:
V t = k 1 ν ( S , y ) e k 1 t , V y = ν y e k 1 t , V S = ν S e k 1 t , V S S = ν S S e k 1 t , V S y = ν S y e k 1 t , V y y = ν y y e k 1 t .
the substitution of Equation (27) into Equation (1) gives the following one-dimensional second-order partial differential equation of function ν ( S , y ) :
0 = 1 2 f 2 ( y ) S 2 ν S S + ρ β S f ( y ) ν S y + 1 2 β 2 ν y y +
r S ν S + α ( m y ) β ρ μ r f ( y ) ν y ( r k 1 ) ν .
using the SYM [19], it can be found that Equation (28) admits one symmetry X 2 = S S . Applying the symmetry S S + k 2 V d V , which is equivalent to the symmetry vector Y 2 , to Equation (28), a second-order differential equation can be obtained. This is shown below in detail.

4.2. Invariant Solution through Symmetry X 2

The symmetry Y 2 is given by:
Y 2 = S S + k 2 ν ν ,
which gives the following characteristic system:
k 2 d S S = d ν ν .
the integration of Equation (31) gives:
k 2 ln S + ln w ( y ) = ln ν ,
i.e.:
ln S k 2 + ln w ( y ) = ln ν .
hence
ν ( S , y ) = S k 2 w ( y ) .
substituting Equation (34) into Equation (26), we can obtain:
V ( t , S , y ) = e k 1 t S k 2 w ( y ) .
one can use Equation (34) to reduce Equation (28) by computing partial derivatives of ν with respect to S and y as follows:
ν y = S k 2 w y , ν S = k 2 S k 2 1 w , ν S y = k 2 S k 2 1 w y , ν S S = k 2 ( k 2 1 ) S k 2 2 w , ν y y = S k 2 w y y .
substituting these partial derivatives into Equation (28) gives:
0 = β 2 w y y + 2 α ( m y ) + 2 ρ β f ( y ) ( k 2 f 2 ( y ) μ + r ) w y + [ ( k 2 2 k 2 ) f 2 ( y ) + 2 ( r k 2 r + k 1 ) ] w .

4.3. Invariant Solution through Symmetry Y 12

The symmetry Y 12 is given by:
Y 12 = t + c S S + k 3 V V .
the characteristic is given by:
d t 1 = d S c S = d V k 3 V .
the first characteristic is solved by integrating:
d t 1 = d S c S ,
to get:
t c = ln S + ln z .
hence:
S = z e c t ,
and so:
z = S e c t .
the second characteristic is solved by integrating:
d t 1 = d V k 3 V ,
to get:
k 3 t = ln V + ln v ( z , y ) ,
which then gives:
V ( t , S , y ) = e k 3 t v ( z , y ) .
now, derivatives are obtained as follows:
V t = k 3 e k 3 t v c S e k 3 t e c t , V S = e k 3 t e c t v z , V y = e k 3 t v y , V y y = e k 3 t v y y , V S S = e k 3 t e ( 2 c t ) v z z , V S y = e k 3 t e c t v z y .
substituting these derivatives into Equation (1), we can find:
0 = 1 2 f 2 ( y ) S 2 e k 3 t e 2 c t v z z + ρ β S f ( y ) e k 3 t e c t v z y + + 1 2 β 2 e k 3 t v y y + r S e k 3 t e c t v z + α ( m y ) ρ β μ r f ( y ) e k 3 t v y r e k 3 t v + e k 3 t v c S e c t v z .
multiplying throughout by 2, cancelling out the term e k 3 t , and imposing the equation z = S e c t leads to:
0 = f 2 z 2 ( y ) v z z + 2 ρ β z f ( y ) v z y + β 2 v y y + 2 ( r c ) z v z + 2 α ( m y ) ρ β μ r f ( y ) v y 2 ( r k 3 ) v .
using the SYM, it can be found that Equation (49) admits the symmetry z z , which is a reduced symmetry [3]. Hence, one can apply the zeroth-order invariant of the symmetry vector z z + w w to Equation (49) to get a reduced form, as follows: Taking the characteristic of z z + k 4 v v , we get:
d z z = d v k 4 v .
integrating gives:
k 4 ln z + ln u ( y ) = ln v ,
and hence:
v ( z , y ) = z k 4 u ( y ) .
the substitution of Equation (52) into (46) gives
V ( t , S , y ) = e k 3 t z k 4 u ( y ) .
now, from Equation (52), we can compute the derivatives of v with respect to z and y as follows:
v z = k 4 z k 4 1 u , v z z = k 4 ( k 4 1 ) z ( k 4 2 ) u , v z y = k 4 z k 4 1 u y , v y = z k 4 u y , v y y = z k 4 u y y .
substituting these into Equation (49) gives the following second-order differential equation:
0 = β 2 u y y + 2 α ( m y ) + 2 ρ β f ( y ) ( k 2 f 2 ( y ) μ + r ) u y + [ ( k 2 2 k 2 ) f 2 ( y ) + 2 ( r k 2 r + k 1 ) ] u .
the solution of differential Equation (55) can be obtained using Maple software. This results in a very complicated solution that contains hypergeometric functions. However, this solution can then be substituted into Equation (53) to get the complete solution V ( t , S , y ) for the original partial differential equation.

5. Symmetry Analysis of the Heston Model

This section is dedicated to performing Lie symmetry analysis for the Heston model. The Black–Scholes equation for the Heston model is given by [3]:
0 = 1 2 Y S 2 V S S + ρ δ Y S V S Y + 1 2 δ 2 Y V Y Y + r S V S + ( α ( m Y ) λ Y ) V Y r V + V t .
in order to ease the process of calculating symmetries for Equation (56), the following change of variables was employed:
Y = y 2 β = δ 2 c 1 = ( α λ ) c 2 = α m β 2
which provided the new differential equation below:
0 = 1 2 y 2 S 2 V S S + β ρ y S V S y + 1 2 β 2 V y y + r S V S + 1 2 c 1 y + c 2 y V y r V + V t .
using the SYM [19], we obtained the given Lie operators:
X 1 = t , X 2 = S S , X V = V V , X b = b ( t , S , y ) V ,
where X V and X b are linear symmetry and infinite symmetry, respectively. Next, the linear combinations of these symmetries were considered to obtain:
Y 1 = X 1 + k 1 X v ,
Y 2 = X 2 + k 2 X V ,
Y 12 = X 1 + c X 2 + k 3 X V .

5.1. Invariant Solution through Symmetry Y 1

The symmetry Y 1 is given by:
Y 1 = t + k 1 V V .
the characteristic system is therefore given by:
k 1 d t 1 = d V V .
integrating both sides, we can obtain:
k 1 t = ln V + ln Φ ( S , y ) ,
which results in:
V ( t , S , y ) = Φ ( S , y ) e k 1 t .
the computation of partial derivatives of V with respect to its parameters gives:
V t = k 1 Φ ( S , y ) e k 1 t , V y = Φ y e k 1 t , V S = Φ S e k 1 t , V S S = Φ S S e k 1 t , V S y = Φ S y e k 1 t , V y y = Φ y y e k 1 t .
substituting these in Equation (58), the function Φ ( S , y ) satisfies the following equation:
0 = 1 2 y 2 S 2 Φ S S + β ρ y S Φ S y + 1 2 β 2 Φ y y + 1 2 c 1 y + c 2 y Φ y r Φ + k 1 Φ .
using the SYM, it can be found that Equation (67) admits one symmetry X 2 = S S . This means that if the symmetry S S + k 2 V d V , which is equivalent to the symmetry vector Y 2 applied to Equation (46) being able to reduce the equation into a second-order ordinary differential equation. This is shown below in detail.

5.2. Invariant Solution through Symmetry X 2

The symmetry Y 2 is given by:
Y 2 = S S + k 2 Φ Φ ,
which gives the following characteristic system:
k 2 d S S = d Φ Φ .
integrating both sides gives:
k 2 ln S + ln W ( y ) = ln Φ
i.e.:
ln S k 2 + ln W ( y ) = ln Φ
and therefore:
Φ ( S , y ) = S k 2 W ( y ) .
substituting Equation (68) where there is Φ into Equation (66), we can find:
V ( t , S , y ) = e k 1 t S k 2 W ( y ) .
one can use Equation (68) to reduce Equation (46) by computing partial derivatives of Φ with respect to S and y as follows:
Φ y = S k 2 W y , Φ S = k 2 S k 2 1 W , Φ S y = k 2 S k 2 1 W y , Φ S S = k 2 ( k 2 1 ) S k 2 2 W , Φ y y = S k 2 W y y .
substituting these partial derivatives into Equation (46) gives:
0 = β 2 W y y + 2 β ρ y k 2 + c 1 y + c 2 y W y + 2 k 1 2 r ( 1 k 2 ) + y 2 ( k 2 2 k 2 ) .
as explained in Section 4, Maple software can be used to obtain asolution of differential Equation (70). The complete solution V ( t , S , y ) of the original partial differential equation can then be found by substituting the obtained solution into Equation (69).

6. Numerical Solutions

The numerical solutions were computed using Maple software, and they are depicted in Figure 1 and Figure 2 below. The parameters in Figure 1 are chosen as ρ = 0.5 , β = 0.7 , k 1 = 1 , k 2 = 0.5 , r = 0.5 , c 1 = 0.01 , c 2 = 0.05 , for the dashed line; ρ = 0.5 , β = 0.7 , k 1 = 1 , k 2 = 0.5 , r = 0.5 , c 1 = 0.01 , c 2 = 0.003 , for solid blue line and ρ = 0.5 , β = 0.7 , k 1 = 1 , k 2 = 0.5 , r = 0.5 , c 1 = 0.01 c 2 = 0.01 , for red dotted line. However, parameters in Figure 2 were chosen as: ρ = 0.5 , β = 0.7 , k 1 = 1 , k 2 = 0.5 , r = 0.5 , c 1 = 0.01 , c 2 = 0.05 , for the dashed line; ρ = 0.5 , β = 0.7 , k 1 = 1 , k 2 = 0.5 , r = 0.5 , c 1 = 0.01 , c 2 = 0.003 , for solid blue line; ρ = 0.5 , β = 0.7 , k 1 = 1 , k 2 = 0.5 , r = 0.5 , c 1 = 0.01 , and c 2 = 0.01 , for the red dotted line.

7. Conclusions

This study looked at the evolution of the solution of the Black–Scholes model for stochastic volatility using the technique known as the modified local one-parameter transformation, and symmetries were obtained and then used to obtain an invariant solution. The model was assumed to follow the Ornstein–Uhlenbeck process, and the Lie symmetry analysis reduced the model to a second-order ordinary differential equation. This process was then applied to the Heston model. The future work in this regard will be to incorporate the dividend yield and observe how the solutions evolve. Another possible extension of the model is to consider an interest rate that is not constant, as an instance interest rate can be considered to be a function of time or a stochastic process.

Author Contributions

These authors contributed equally to this work. All authors have read and approved the final version of this manuscript.

Funding

This research was financially supported by the research office of the University of Zululand.

Acknowledgments

This research was financially supported by the research office of the University of Zululand. The Authors thank the reviewers for their comprehensive comments and recommendations, which significantly improved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Miller, R.M. Option Valuation, Economic and Financial Modelling with Mathematica; Varian, H.R., Ed.; Springer: Berlin/Heidelberg, Germany, 1993; pp. 32–58. [Google Scholar]
  2. Merton, R.C. Theory of Rational Option Pricing. Bell J. Econ. Manag. Sci. 1973, 41, 141–183. [Google Scholar] [CrossRef] [Green Version]
  3. Paliathanasis, A.; Krishnakumar, K.; Tamizhmani, K.M.; Leach, P.G.L. Lie Symmetry Analysis of the Black-Scholes-Merton Model for European Options with Stochastic Volatility. Mathematics 2015, 4, 28. [Google Scholar] [CrossRef] [Green Version]
  4. Bluman, G.W.; Anco, S.C. Symmetry and Integration Methods for Differential Equations; Springer: New York, NY, USA, 2002. [Google Scholar]
  5. Bluman, G.W.; Kumei, S. Symmetries of Differential Equations; Springer: New York, NY, USA, 1989. [Google Scholar]
  6. Hydon, P. Symmetry Methods for Differential Equations: A Beginner’s Guide; Cambridge Texts in Applied Mathematics; Cambridge University Press: Cambridge, UK, 2000; pp. i–vi. [Google Scholar]
  7. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer: New York, NY, USA, 1986. [Google Scholar]
  8. Sedov, L.I. Similarity and Dimensional Methods in Mechanics; Academic Press: New York, NY, USA, 1959. [Google Scholar]
  9. Paliathanasis, A.; Leach, P.G.L. Symmetry analysis for a fourth-order noise-reduction partial differential equation. arXiv 2020, arXiv:2008.06323. [Google Scholar]
  10. Matadi, M.B.; Zondi, P.L. Lie Theoretic Perspective of Black-Scholes Equation under Stochastic Heston Model. Int. J. Appl. Math. 2020, 33, 753–764. [Google Scholar] [CrossRef]
  11. Matadi, M.B. Singularity and Lie group Analyses for Tuberculosis with Exogenous Reinfection. Int. J. Biomath. 2015, 8, 1–12. [Google Scholar] [CrossRef]
  12. Matadi, M.B. Invariant solutions and conservation laws for a pre-cancerous cell population model. J. Interdiscip. Math. 2020, 23, 1121–1140. [Google Scholar] [CrossRef]
  13. Matadi, M.B. Lie Symmetry Analysis of Early Carcinogenesis Model. Appl. Math. E-Notes 2018, 18, 238–249. [Google Scholar]
  14. Matadi, M.B. Symmetry and conservation laws for tuberculosis model. Int. J. Biomath. 2017, 10, 1750042. [Google Scholar] [CrossRef]
  15. Matadi, M.B. The Conservative Form of Tuberculosis Model with Demography. Far East J. Math. Sci. 2017, 102, 2403–2416. [Google Scholar] [CrossRef]
  16. Oke, S.I.; Matadi, M.B.; Xulu, S.S. Optimal control analysis of a mathematical model for breast cancer. Math. Comput. Appl. 2018, 23, 1–28. [Google Scholar]
  17. Masebe, T.P. A Lie Symmetry Analysis of the Black-Scholes Merton Finance Model Through Modified Local One-Parameter Transformations; University of South Africa Press: Unisa, South Africa, 2014; pp. 37–80. [Google Scholar]
  18. Oliveri, F. Lie Symmetries of Differential Equations: Classical results and recent contributions. Symmetry 2010, 2, 658–706. [Google Scholar] [CrossRef] [Green Version]
  19. Dimas, S.; Tsoubelis, D. SYM: A new symmetry-finding package for Mathematica. In Proceedings of the 10th the International Conference in Modern Group Analysis, Larnaca, Cyprus, 24–31 October 2004. [Google Scholar]
  20. Gazizov, R.K.; Ibragimov, N.H. Lie symmetry analysis of differential equations in finance. Nonlinear Dyn. 1997, 17, 387–407. [Google Scholar] [CrossRef]
Figure 1. Numerical solution of the invariant solution for the Heston model Equation (52); the parameters were chosen as follows:
Figure 1. Numerical solution of the invariant solution for the Heston model Equation (52); the parameters were chosen as follows:
Symmetry 13 00847 g001
Figure 2. Numerical solution of the invariant solution for the Heston model Equation (52); the parameters were chosen as follows:
Figure 2. Numerical solution of the invariant solution for the Heston model Equation (52); the parameters were chosen as follows:
Symmetry 13 00847 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Matadi, M.B.; Zondi, P.L. Invariant Solutions of Black–Scholes Equation with Ornstein–Uhlenbeck Process. Symmetry 2021, 13, 847. https://doi.org/10.3390/sym13050847

AMA Style

Matadi MB, Zondi PL. Invariant Solutions of Black–Scholes Equation with Ornstein–Uhlenbeck Process. Symmetry. 2021; 13(5):847. https://doi.org/10.3390/sym13050847

Chicago/Turabian Style

Matadi, Maba Boniface, and Phumlani Lawrence Zondi. 2021. "Invariant Solutions of Black–Scholes Equation with Ornstein–Uhlenbeck Process" Symmetry 13, no. 5: 847. https://doi.org/10.3390/sym13050847

APA Style

Matadi, M. B., & Zondi, P. L. (2021). Invariant Solutions of Black–Scholes Equation with Ornstein–Uhlenbeck Process. Symmetry, 13(5), 847. https://doi.org/10.3390/sym13050847

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop