Medical, Pharma, Engineering, Science, Technology and Business

Department of Mathematics and Applied Mathematics, University of Pretoria, Pretoria 0002, Republic of South Africa

- *Corresponding Author:
- Eben Mare

Department of Mathematics and Applied Mathematics

University of Pretoria, Pretoria 0002

Republic of South Africa

**Tel:**+27-12-420-4111

**E-mail:**[email protected]

**Received Date:** July 18, 2017; **Accepted Date:** July 25, 2017; **Published Date:** July 31, 2017

**Citation: **Nteumagn BF, Pindza E, Mare E (2017) Symmetry Analysis of Options
Pricing with Transactions Costs Driven by Fractional Brownian Noises. J Appl
Computat Math 6: 356. doi: 10.4172/2168-9679.1000356

**Copyright:** © 2017 Nteumagn BF, et al. This is an open-access article distributed
under the terms of the Creative Commons Attribution License, which permits
unrestricted use, distribution, and reproduction in any medium, provided the
original author and source are credited.

**Visit for more related articles at** Journal of Applied & Computational Mathematics

We provide a closed-form solution for the European and Asian option pricing models when the source of randomness is a fractional Brownian motion as opposed to the geometric Brownian motion. In addition to the source of randomness, transaction costs are considered to be non-negligible. For the case of the European option, proportional transaction costs hide in the volatility and do not change the form of the model. The construction of the solution is based on the symmetries of the model. The model for Asian options has an additional parameter that makes the volatility time-dependent, which complicates the solution process. However, we are still able to obtain solutions using Lie symmetry methods.

Lie symmetries; European option; Fractional brownian; Transaction cost

In the construction of several models, we make explicit and implicit assumptions. However, quite often, when we solve these models, we overlook these assumptions or they are simply “diluted” in the process [1]. For models encoded as differential equations, assumptions that translate into the invariance of underlying equations are particularly useful in the search of their solutions. This is the main thrust of Lie symmetry analysis of differential equations pioneered by Sophus Lie [2-18]. In such analysis, one algorithmically looks for infinitesimal transformations (i.e., transformations depending on a small parameter and enjoying additional properties [2,9,17] of both dependent and independent variables that does not change the underlying differential equation. These transformations are then systematically used for finding solutions.

An infinitesimal transformation of both dependent and independent variables that leaves a differential equation unchanged up to the first-order in the small parameter is called an infinitesimal symmetry [17] of this equation. For a given differential equation, the set of all its infinitesimal symmetries form an infinitesimal group (Lie algebra in the modern terminology) [19]. Perhaps the most useful property of a symmetry of a differential equation is that it transforms a solution into another one. Thus, for instance, by starting with a trivial solution, one may construct a non-trivial one by acting on the former with a symmetry [2]. Also, for equations belonging to the same family, those that can be invertibly mapped to each other have isomorphic symmetry algebra. We shall employ this last property for transforming a complicated equation into one which is easier to solve [20]. Lie symmetry analysis has a broad range of applications. This includes among others, mathematical biology, mathematical physics and financial mathematics to name just a few. In this paper, we focus our attention on the application of Lie symmetry analysis to financial mathematics, and more precisely to the theory of option pricing.

Option pricing and replication has attracted both the business and the academic audience since the development of the Nobel price winning Black-Scholes differential equation. Significant success has been registered in the attempt to solve the model numerically. However, it is not always possible to find practical numerical algorithms to solve higher dimensional differential equations. Because of its ability to reduce the dimension of the problem, Lie symmetries analysis promises to be a very good instrument to mitigate Bellman’s “curse of dimensionality” [21]. It constitutes a great tool for analytical approximation of the solution even when numerical analysis has not been found to be conclusive. Thus far, it has been applied to the celebrated Black-Scholes equation with success.

(1.1)

Indeed, Gazizov and Ibragimov [22] paved the way in the application of Lie symmetry analysis to solve financial models [23]. Their research was based on the basic BlackScholes equation (1.1). In their work, they applied Lie symmetries to transform the BlackScholes equation into the heat equation and solved it. Pooe et al. [24] deduced a fundamental solution to zero-coupon bonds. Since the introduction of Lie symmetries to finance, many models have been designed to price and hedge options accurately, and models are becoming more sophisticated with the development of new techniques and the evolvement of technology (see for example, [12,13]).

The option pricing theory from inception rested on the arbitrage argument which stipulates that by continuously hedging a portfolio of stocks and risk-free bonds, one can exactly replicate the return on a stock [25]. By doing so, the option value must then be equal to that of the replicating portfolio [5,23,11]. In reality, one cannot completely hedge a portfolio of stocks and bonds without incurring any costs. This is because continuous trading and replication has been proven to be very expensive, even though it is possible to consider very small proportional transaction cost [11,26-28]. The consideration of transaction costs in option pricing is more recent as assumptions of the classical models have minimized these costs in the past. The derivation of the option pricing equation relies heavily on the random process followed by the underlying stock price. In the case of the celebrated Black-Scholes equation, one assumes that the stock follows the classical geometric Brownian motion [5,24]. This process assumes independent increments, which makes it different from the fractional Brownian motion (fBm) where there is a serial correlation in the increments. This correlation allows for predictability in the model. The similarity between the two processes is the fact that they are both geometric.

It has been shown [8] that random changes can be accurately measured by the use of a parameter known as the Hurst parameter. The construction of a dam on the Nile river in the early 1950s saw the birth of the Hurst exponent. The measure of the water levels of increase or decrease was known to follow a random walk. In other words, the next level of increase or decrease is independent of the previous one. However, Edwin Hurst, hydrologist in the Nile dam project, showed that the long term memory of the system indicates an autocorrelation of the time series of water levels, and the rate at which these decrease as the lag between pairs of values increase [8]. The movement of stock price is similar to this model. We therefore exploit the results by Leland [11] for the derivation of our model. Earlier in 2014, we achieved an accurate numerical solution of our model with exponential convergence ([19]) by applying spectral methods to analyze the problem. In this paper, we use Lie symmetries to determine the exact solution, which sets us a step closer to the calibration of our model.

We first give some preliminaries on Lie symmetries of differential equations in Section 2. Secondly, we determine the symmetries of the European and Asian option pricing models and obtain their analytical solutions in Section 3. In Section 4 we apply the Lie symmetry technique to reduce the PDEs and provide solutions. Finally, we provide concluding remarks in Section 5.

Our aim in this section is to provide a brief introduction to Lie symmetry algorithm and define what are equivalence transformations. We shall limit ourselves to scalar second-order PDEs as our models will fall in this category. We shall sometimes sacrifice mathematical rigor for the sake of clarity.

**Lie symmetry algorithm**

Classical methods for solving differential equations are ad-hoc in nature. Lie symmetry analysis offers a unified approach to the solution of differential equations and the reduction of the number of dependent variables. Intuitively, a symmetry of a differential equation is an invertible change of independent and dependent variables that leave the underlying equation unchanged. Finding all the symmetries of a given differential equation is a difficult if not impossible task in many situations. However, as Lie [2,17] discovered, those symmetries that depend on a small parameter and which form a local one-parameter group [2,17] can be systematically calculated in fairly general situations. In this subsection, we specialize Lie’s symmetry algorithm to the class of problems that interest us in this work. Also, we highlight a few properties of Lie symmetries. Consider a second-order partial differential equation (PDE) in one dependent variable u and two independent variables t and x.

E(t, x, u, u_{(1)},u_{(2)})=0 (2.2)

where E is a sufficiently smooth function of its argument, u_{(1)}={ut, ux}, u_{(2)}={utt, utx, uxx}, and the indexes stands for partial differentiations. An invertible transformation of both the independent and dependent variables

(2.3)

(2.4)

(2.5)

is an infinitesimal symmetry or Lie symmetry if

E(t, x, u, u_{(1)},u_{(2)}) =O(∈^{2}). (2.6)

Whenever, eqn. (2.2) is satisfied,

In order to obtain a practical definition of infinitesimal symmetry, we need to transform eqn. (2.6) into an equivalent form by means of Taylor’s expansion in ?. Thus, we must obtain first the Taylor expansions of the derivatives appearing in eqn. (2.6). This is our immediate goal.

We may compute the differential d¯ u in two different ways. The first way comes from the very definition of a differential:

(2.7)

Where we have used eqns. (2.3)-(2.4), Dt and Dx are the operators of total di?erentiations with respect to t and x respectively:

The second way of computing d¯ u is by using eqn. (2.5):

(2.10)

From eqn. (2.7) and eqn. (2.10), we infer the following system of two linear equations in the two unknowns and

(2.11)

(2.12)

By solving the system in eqn. (2.11)-(2.12), we obtain

(2.13)

(2.14)

where

(2.15)

(2.16)

By replacing ¯ u with ¯ u¯ t or ¯ u¯ x in the preceding calculations, we arrive at the following expansions of the second derivatives:

(2.17)

(2.18)

(2.19)

in which are given by

(2.20)

(2.21)

(2.22)

The equations (2.15)-(2.16) are known as first-order prolongation formulas whereas Equation. (2.20)-(2.22) are the second-order prolongation formulas. Note the pattern between first order and second-order prolongation formulas. Such insight may be utilized to infer high order prolongation expressions.

Now let us expand eqn. (2.6) to the second-order in ?:

(2.23)

whenever eqn. (2.2) is satisfied, where

(224)

(2.25)

(2.26)

We shall also refer to X as infinitesimal symmetry of eqn. (2.2). Now, eqn. (2.23) implies that

*X*^{[2]}(E) = 0 (2.27)

Whenever eqn. (2.2) is satisfied. Equation (2.27) together with the phrase “whenever eqn. (2.2) is satisfied” may be written compactly as

(2.28)

Conversely, it can be shown (see for instance [2,3,17]) that if eqn. (2.2) can be locally solved with respect to one of its highest derivative and eqn. (2.28) is satisfied, so is eqn. (2.6). Therefore, modulo the aforementioned local solvability condition, eqn. (2.28) characterizes infinitesimal symmetries of eqn. (2.2). That is why it is called the linearized symmetry condition for infinitesimal symmetries.

After substituting one of the highest derivatives of eqn. (2.2) into the linearized symmetry condition (2.28), we end up with an equation which is a polynomial in the remaining derivatives of u. The latter equation splits into an over-determined system of linear PDEs for the unknowns τ, ξ and η. There are several computer algebra packages for generating automatically this set of linear PDES [12,13] so that one seldom generates them by hand.

Owing the structure of the determining equation (2.28), the infinitesimal symmetries enjoy the following additional properties:

• If X_{1} and X_{2} are two Lie symmetries of eqn. (2.2) and c_{1} and c_{2} are real numbers, then c_{1} X_{1}+c_{2}X_{2} is also Lie symmetries of (2.2).

• If X_{1} and X_{2} are two Lie symmetries of eqn.(2.2), then their Lie bracket

[X_{1},X_{2}]=X_{1}(X_{2})−X_{2}(X_{1})

is also a Lie symmetry of eqn. (2.2).

Thanks to the above two properties, the collection of all the infinitesimal symmetries of eqn. (2.2) form a Lie algebra [2,16,17].

**Equation depending on arbitrary parameters: equivalence transformations**

Here we focus on introducing equivalence transformations and we emphasize some of their applications to problems pertinent to this paper. As, in the previous subsection, we shall limit ourselves to a scenario relevant to this work. Notations that are not introduced here are those of the previous sections. For an in-depth study of the equivalence problem, we refer the reader to the work of Ovsiannikov [18] and former students [2,17]. Consider a family of differential equations

E(t,x,u,u_{(1)},u_{(2)},Λ)=0, (2.29)

which is parameterized by a set Λ made possibly of arbitrary constants and arbitrary functions of t, x, u and u_{(1)}. An equivalence transformation of eqn. (2.29) is an invertible change of dependent and independent variables which preserves its differential structure. The key di?erence between a symmetry and an equivalence transformation is that, an equivalence transformation is required only to preserve the structure of the set of arbitrary parameters Λ: If Λ contains an arbitrary number α, then any equivalent equation’s parameter set, , must contain an arbitrary parameter ¯ α at the same position in the differential structure of the equation. If Λ contains a function f(t,x,u), then the parameter set of the equivalent equation, , must contain a function at the same location in the differential structure of the equation. It is straightforward to see that a symmetry of eqn.(2.29) is always and equivalence transformation of eqn. (2.29). However, the converse is not in general true. Also, two equations of the class in eqn. (2.29) that can be mapped to each other via an equivalence transformation have isomorphic Lie symmetry algebras [2,3,16,17]. The latter results are often used to detect isomorphic equations within the class in eqn. (2.29) and to even construct explicitly the equivalence transformation which realizes this isomorphic.

Equivalence transformations induce an equivalence relation on the family of equations (2.29): Two equations of the family are equivalent if and only if they can be mapped to each other by using an equivalence transformation. Finding all the equivalence classes modulo this equivalence transformation is the group classification problem associated to eqn. (2.29). In the group classification problem, one seeks very simple representative of each equivalence class, which in general, carries fundamental properties of its members. It is generally good to choose a simple enough representative to permit subsequent analysis. Indeed, once the properties of a representative are elucidated, they can always be mapped to those of a given representative via an appropriate equivalence transformation. This is mainly how we are going to exploit equivalence transformations in this work.

The above narrative begs two questions of practical importance: How does one find all the equivalence transformations of the family (2.29)? How does one verify that two equations of the family in eqn. (2.29) are equivalent?

The problem of calculating all the equivalence transformations of the family in eqn. (2.29) is as difficult as finding all its symmetries. In some particular instance, one may find all the equivalence transformations through the use of the chain rule and reasoning about the differential structure of the family [2,3,16,17]. As for Lie symmetries, if we restrict ourselves to continuous groups of equivalence transformations, one arrives at an algorithm [16] which is similar to Lie’s symmetry algorithm in its derivation and application. Concerning the issue of the equivalence of two equations of the family in eqn. (2.29), one uses the fact that two equation are equivalence if and only if their symmetry Lie algebras are isomorphic i.e. they have the same structure constants after an appropriate choice of their bases [2,3,16,17]. We shall take advantage of this result in the sequel.

This section is dedicated to the derivation of both vanilla European and Asian options PDEs under fractional Brownian noises. We only provide a detailed derivation of the former since the latter follows the same logic mutatis mutandis.

**Black-Scholes model under fractional Brownian motion**

We provide the derivation of the European option pricing PDE in fractional Brownian motion with transactions costs.

**Definition:** A fractional Brownian motion (fBm) BH(t) with Hurst exponent H ∈ (0,1) is a Gaussian process with mean zero and covariance [4].

(3.30)

Under this model, the trader predicts a price for the stock at time t as

(3.31)

which implies a stock model

*dS* = *μSdt* + *σdB ^{H}* (

We assume that the European call V (t,S) on the underlying stock S at time t with exercise T and exercise price E (also known as the strike price or just strike) has boundary conditions.

(3.33)

Indeed, if at maturity, the stock price is greater than the strike price, then the option is exercised. If instead the stock price is less than the strike, then the buyer exercises his right not to buy and the option is worthless.

We also considers the risk free bond price D(t) whose evolution is given by

dD(t)=rD(t)dt (3.34)

We build a simple portfolio with value Π, say, made of stocks and risk-free bonds. We suppose that it costs us a proportion k of the turnover to replicate our portfolio (see for example [11]). If we replicate the portfolio with X_{1}(t) units of stocks and X_{2}(t) units of risk free bonds, then the portfolio Π will be valued at

Π=X_{1}(t)S+X_{2}(t)D(t) (3.35)

Assume small changes δΠ, δS and δD in Π, S and D, respectively, caused by a small change in time δt. This change in the portfolio value is obtained at a transaction cost, which is assumed to be proportional to turnover, i.e., . Hence we have from (3.35), we have

(3.36)

Where δD(t) is a small change in risk-free bonds, δX1(t) a small change in the number of stocks held. The replication interval δt is chosen to be minimal. Therefore, by the Taylor theorem [26], expanding V in the δ − neighbourhood of t we can write

(3.37)

and

(3.38)

Now we have

(3.39)

so

(3.40)

Similarly,

(3.41)

for the higher even powers of δS, and since (δt) is chosen to be very small, higher powers are negligible.

Equation (3.38) gives

(3.42)

Similarly, (3.37) and (3.36) yield

(3.43)

If V (t,S) is the value of the option replicated by the above portfolio Π, i.e., V=Π to reduce arbitrage and remain consistent with economic equilibrium [26,27], then by WickIto’s lemmas as in ref. [7]

(3.44)

Since V=Π, δV=δΠ. From (3.35) and (3.36), we have

(3.45)

(3.46)

So we have V=ΔS+B^{H}(t) and thus

(3.47)

Hence, by the fractional Leland [11] approximation

(3.48)

Where . Finally, we get (3.49)

For vanilla when the rate at which the loss occurs accelerates. That is why it is widely assumed that is always positive. Therefore, Equation (3.49) becomes

(3.50)

where

(3.51)

Without loss of generality, we only consider the case of a European call option with terminal condition.

V (S,T)=max{S−E,0} (3.52)

We transform the problem into an initial value problem with the change of variables τ=T−t. This leads to the PDE.

(3.53)

Subject to the initial conditions V (S,0)=max{S−E,0}. The option pricing with transaction cost has a same form as the classical Black- Scholes PDE, except for the transformation τ=T−t and the volatility parameter which is now given by eqn. (3.51). Proportional transaction costs do not a?ect the linearity of the equation when a vanilla option is considered. For the purpose of this work, we will focus on vanilla option pricing. The volatility structure of the fBm model depends on three important parameters: When we vary δt, the optimal hedging interval, we are able to define the minimal volatility for this model as well as the equivalent time interval. The volatility structure is given in equation (3.51) above and **Figure 1** plots the impact of changes in the Hurst parameter H on the volatility of the underlying stock.

(3.54)

Setting this to zero and solving for H leads to meaningless results in view of our model. However if we solve this in terms of δt, we obtain

(3.55)

When. The profile of the minimal volatility is shown in **Figure 2**.

The complete volatility structure is shown in **Figure 3**. Notice that if k=0 and H=0.5, we have .

**Derivation of the geometric average Asian options pricing model under fractional Brownian motion**

In addition to the assumptions of the European option, Asian options posit that the stock price x(t) depends on historical prices. However the steps of the derivation of the equation are similar to those of the European option differential equation provided in the above section. We omit the details of the derivation for this reason.

The model in this context is given by

(3.56)

Where this time the volatility is time-dependent and is given by

(3.57)

the fractional Leland number, while t and S remain the time and price of the stock respectively. here denotes the geometric average of the stock price over the time interval [t_{0},t], V=V (S,t,J) and of course, the Hurst parameter H. We can now proceed to do the symmetry analysis of these models. We start with the European option.

We dedicate this section to solving the equations for pricing European and Asian options in eqn. (3.53) and (3.56). We first consider European option case, and later, the Asian options under fractional Brownian motion.

**European call option under fractional Brownian motion**

The first step is to transform Equation (3.53) into a heat like equation using the Lie symmetry analysis. Thanks to the computer algebra package Symteric [12,13] one does not need to calculate the symmetries by hand, but a line of code success to verify that the symmetries of the PDE in eqn. (3.53) are as in ref. [10].

Equation (3.53) is reduced to the heat equation by applying the transformations

(4.58)

where , and . In fact, the transformations of variables above give

(4.59)

Substituting these variables and their derivatives into Equation (5.53) yields the one dimensional heat equation

(4.60)

together with the initial condition

(4.61)

with solution

(4.62)

Defining new variables d_{1} and d_{2}, then converting the variables back into the variables in our model, we have

(4.63)

The solution of the PDE is obtained by backward substitution from w(y,τ) to V (S,y). We get

V (S,t)=SN(d_{1})−Ee^{r(t−T)}N(d_{2}), (4.64)

where N is the standard normal distribution probability mass function, which satisfies our equation and the terminal condition. Note that when H=0.5 and k=0.002, i.e., 0.2% transaction cost, we obtain the solution to the classical Black-Scholes equation. We vary the Hurst parameter and plot the solution in **Figure 4**.

We now consider the Asian option pricing problem in the same regime.

**Lie symmetry classification of the geometric Asian option pricing model under fractional Brownian motion**

In this subsection, we apply the theory described above to the Asian option pricing problem in fBm. We consider the contextual model derived above in eqn. (3.56) and rewrite it for readability purposes as,

(4.65)

subject to the pay on

V (S,T,J)=max{J −K,0}

where

(4.66)

is the geometric average of the price of the underlying S since inception τ=0 till time τ=t, K is the strike, r is the risk-free rate, which in most emerging markets is nonzero. σ is a constant volatility. The minimum tradable interval (δt) is very important for the practical aspect of our model. Note that if this is zero, it means that we are hedging in continuous time. It has, however, been shown (see for example [14]) that in continuous time we face the problem of arbitrage. In order to avoid the problem of arbitrage, we only consider the case where (δt)≠0. The same condition applies to k, the transaction cost. This is a proportion of turnover that would be charged for replicating the portfolio. If this is allowed to be k=1, we would have more symmetries. However, it is clearly an unrealistic scenario. If the turnover on a transaction is equal to the cost of that transaction, it would never be worth entering such transactions. Despite the above, we find that in eqn. (4.65) in its full glory does not lead to any results. However, we distinguish the below cases:

We have for this case, , which gives the simplified model

(4.67)

which admits the symmetries Lie algebra generators

(4.68)

in addition to the infinite dimensional symmetry f (t, S, J) ∂V, where

Applying the symmetry X_{2} to perform the reduction, we have the quasi-linear PDE

(4.69)

Which gives us the transformation variables

y (S,J)=In S–In J, V(t,S,J)=u(t,y) (4.70)

Where u satisfies the (1+1) PDE

(4.71)

The symmetry Lie algebra of eqn. 4.71 is spanned by the following operators:

(4.72)

(4.73)

Where f (t,y) satisfies the PDE

Reduction of (4.71) is therefore possible by means of the transformation τ=−4σ2/t5, z=y/t2 −rt−σ2t,

(4.74)

to the heat equation

(4.75)

and the solution of the latter is

(4.76)

Substituting back the transformation concludes this case.

Case k=0, (δt)≠ 0

In the case of no transactions cost, we find no symmetry that enables us to perform reduction of the PDE.

This condition yields

(4.77)

with basis for the symmetries Lie algebra given by

(4.78)

as well as the infinite dimensional symmetry. The reduction process using X_{1} yields the transformation

(4.79)

We substitute the above into (4.77) to obtain the 1+1 PDE

(4.80)

Which admits the following basis for the symmetry algebra.

(4.81)

where f satisfies 2(t2H+1 −t2)fξξ=0. It is worth noting that the generators of symmetries above were obtained by transforming those of the heat equation. Using the normal algorithm for generating the symmetries in Section 2, we obtain determining equations, of which the closed form solutions are difficult to construct equation (4.80) can be transformed to the heat equation by the use of the change of independent variables

(4.82)

and the solution is as in case in eqn. 4.2.1. Applying the solution process in as per Case 4.2.1 concludes this section.

European and geometric average Asian options were studied in this work. We successfully constructed the symmetries of the European option pricing model and provided the mathematical tools for mapping this model to a known PDE. Despite the transaction cost and the source of randomness, the model keeps similar algebraic structure to the classical Black-Scholes equation. In other words the model with transaction cost, when considering a problem of a vanilla option, keeps its classification. Some suggestions were given to relax the assumption on the sign of Γ, however, the problem becomes ill-posed if Γ<0.

For geometric average Asian options, we have provided a classification of the PDE and closed form solutions were provided whenever it was possible. We notice that though there has been some progress in technology, it is still relevant and important to keep the rigor of the algorithm. This is due to the fact that though some computer algebras are capable to handle generic systems, it is di?cult to obtain symmetries of some PDEs. As an example, with the parametric time dependency of the PDE for the Asian option, we were unable to apply Sym to obtain the symmetries. We had to fall back on Ya Lie, in order to regain control of the solution process and in the end, we realized that not all the determining equations were solvable using the normal solution process in Mathematical. It is in this regard that we opted to exploit a well known use of symmetries to transform the symmetries of the heat equation.

Methods that have been used so far to obtain the solutions of the above problems are ad hoc in nature. We acknowledge that the application of Lie symmetries analysis provided here is systematic and can be automated. Applying equivalence transformations to the reduced model, we transformed these to the heat equation and provided a solution thereof. We have limit ourselves to European option pricing with constant volatility. We know, however, that in real market conditions, the volatility is not constant. Considering stochastic volatility in fBm would be an extension of this work. Additionally, it will be interesting to now calibrate these models and estimate H for a specific asset class in a given period of time. A question of interest may be: Could non-market conditions dictate the level of the Hurst parameter? In other words, what would have been the e?ect of the 2008 crisis to H for a specific product type? Can H capture a similar market reaction to another crisis in future? Answers to these question could be an interesting extension to the building of these models.

This work was partially financed by Brad Welch at Ridge Cape Capital. This financial assistance enabled BFN to present this work at AIMS in Orlando. BFN would like to express his sincere appreciations to Prof. C. Wafo Soh for his feedback and advice during the preparation of this manuscript.

The authors declare that there is no conflict of interest regarding the publication of this paper.

- Zhang Y, Pan D, Zhou SW, Han M (2014) Asian option pricing with transaction costs and dividend’s under the fractional Brownian motion model. Journal of Applied Mathematics Hindawi Publishing Corporation, p: 8.
- Bluman GW, Kumei S (1989) Symmetries and di?erential equations. Applied Mathematical Sciences Springer-Verlag 81: 229-237.
- Sophocleous C, OHara JG, Leach PGL (2011) Symmetry analysis of a model of stochastic volatility with time-dependent parameters. Journal of Computational and Applied Mathematics 235: 4158-4164
- Nualart D (2006) Fractional Brownian motion: Stochastic calculus and applications. Proceedings of the International Congress of Mathematicians 3: 1541-1560.
- Hull JC (2012) Options, futures and other derivatives. Pearson, 8th edition.
- Mandelbrot BB, Van Ness JW (1968) Fractional Brownian motion, fractional noises and applications. SIAM Review 10: 422-437.
- Dai W, Heyde CC (1996) It’s formula with respect to fractional Brownian motion and its application. Journal of Applied Mathematics and Stochastic Analysis 9: 439-448.
- Hurst HE (1987) Long-term storage capacity in reservoirs, transactions of the American Society of Engineers116: 770-799.
- Ibragimov NH, Kara AH, Mahomed FM (1997) Lie-Backlund¨ and Noether symmetries with applications. Nonlinear Dynamics 15: 115-136.
- Ibragimov NH, Gazizov RK (1998) Lie symmetry analysis of di?erential equations in finance. Nonlinear Dynamics 17: 387-407.
- Leland HE (1985) Option pricing and replication with transaction costs. Journal of Finance 40: 1283-1301.
- Dimas S, Tsoubelis D (2004) SYM: A new symmetry-finding package for Mathematica. Proceedings of the 10th International Conference in Modern Group Analysis, p: 24-30.
- Dimas S, Tsoubelis D (2006) A new Mathematica-based program for solving over determined systems of PDE. Proceedings of the 8th International Mathematica Symposium, p: 19-23.
- Cheridito P(2003) Arbitrage in fractional Brownian motion models. Finance and Stochastics 7: 533-553.
- Wilmott P, Howison S, Dewynne J (1995) The Mathematics of financial derivatives. Cambridge University Press.
- Mahomed FM (2008) Complete invariant characterization of scalar linear (1+1) parabolic equations. Nonlinear Mathematical Physics 1: 112-123 .
- Olver PJ (1986) Applications of Lie groups to di?erential equations. Graduate texts in Mathematics.
- Ovsiannikov LV (1982) Group analysis of di?erential equations.Academic Press, p: 482.
- Nteumagne BF, Pindza E, Mare E (2014) Applying the barycentric Jacobi spectral method to price options with transaction costs in a fractional Black-Scholes framework. Journal of Mathematical Finance 4: 35-46.
- Ablamovitz A, Stagun N (1970) Handbook of Mathematical functions with formulas, graphs and Mathematical tables.
- Houle ME, Kreigel HP, Kroger P, Schubert E, Zimek A (2001) Quality of similarity rankings in time series. Advances in spacial and temporal database 6849: 422.
- Gazizov RK, Ibragimov NH (1998) Lie symmetry analysis of di?erential equations in finance. International Journal of Nonlinear Dynamics 17: 387-407.
- Goard J (2000) New solutions to the bond-pricing equation via Lie’s classical method. Mathematical and Computer Modelling 32: 299-313.
- Pooe CA, Mahomed FM, Wafo Soh C (2004) Fundamental solutions for zero-coupon bond-pricing models. Nonlinear Dynamics 36: 69-76.
- Paulian AD, Zaitsev VF (2002) Handbook of exact solutions for ordinary di?erential equations. Chapman & HALL CRC.
- Wang XT (2010) Scaling and long-range dependence in option pricing I: Pricing European option with transaction costs under fractional Black-Scholes model Physics A 389: 438-444.
- Wang XT, Zhu EH, Tang MM, Yan HG (2010) Scaling and long-range dependence in option pricing II: Pricing European option with transaction costs under the mixed Brownian-fractional. Brownian model Physics A 389: 445-451.
- Wang XT (2010) Scaling and long-range dependence in option pricing I: Pricing European option with transaction costs under the multifractional. Black-Scholes model Physics A 389: 789-796.

Select your language of interest to view the total content in your interested language

- Adomian Decomposition Method
- Algebraic Geometry
- Analytical Geometry
- Applied Mathematics
- Axioms
- Balance Law
- Behaviometrics
- Big Data Analytics
- Binary and Non-normal Continuous Data
- Binomial Regression
- Biometrics
- Biostatistics methods
- Clinical Trail
- Complex Analysis
- Computational Model
- Convection Diffusion Equations
- Cross-Covariance and Cross-Correlation
- Differential Equations
- Differential Transform Method
- Fourier Analysis
- Fuzzy Boundary Value
- Fuzzy Environments
- Fuzzy Quasi-Metric Space
- Genetic Linkage
- Hamilton Mechanics
- Hypothesis Testing
- Integrated Analysis
- Integration
- Large-scale Survey Data
- Matrix
- Microarray Studies
- Mixed Initial-boundary Value
- Molecular Modelling
- Multivariate-Normal Model
- Noether's theorem
- Non rigid Image Registration
- Nonlinear Differential Equations
- Number Theory
- Numerical Solutions
- Physical Mathematics
- Quantum Mechanics
- Quantum electrodynamics
- Quasilinear Hyperbolic Systems
- Regressions
- Relativity
- Riemannian Geometry
- Robust Method
- Semi Analytical-Solution
- Sensitivity Analysis
- Smooth Complexities
- Soft biometrics
- Spatial Gaussian Markov Random Fields
- Statistical Methods
- Theoretical Physics
- Theory of Mathematical Modeling
- Three Dimensional Steady State
- Topology
- mirror symmetry
- vector bundle

- Total views:
**399** - [From(publication date):

September-2017 - May 22, 2018] - Breakdown by view type
- HTML page views :
**352** - PDF downloads :
**47**

Peer Reviewed Journals

International Conferences 2018-19