Analytical Solutions of the Black–scholes Pricing Model for European Option Valuation via a Projected Differential Transformation Method

In this paper, a proposed computational method referred to as Projected Differential Transformation Method (PDTM) resulting from the modification of the classical Differential Transformation Method (DTM) is applied, for the first time, to the Black–Scholes Equation for European Option Valuation. The results obtained converge faster to their associated exact solution form; these easily computed results represent the analytical values of the associated European call options, and the same algorithm can be followed for European put options. It is shown that PDTM is more efficient, reliable and better than the classical DTM and other semi-analytical methods since less computational work is involved. Hence, it is strongly recommended for both linear and nonlinear stochastic differential equations (SDEs) encountered in financial mathematics.


Introduction
Pricing of options is a key aspect of financial mathematics and financial engineering.In 1973 Black and Scholes derived the most famous and significant valuation model, known as Black-Scholes Model for options [1].The model is used for European or American options-be it a call option or a put option.The model is based on some assumptions, among such are the no-arbitrage opportunities, no inclusion of transaction costs associated with hedging, the asset price is lognormally distributed, the drift and the volatility rates are assumed constants, trading of all securities and derivatives are assumed continuous [2].
The assumption of the volatility as a constant function has really posed a challenge in option valuation using the Black-Scholes Model, since it is not the case in reality.
In a bid to address parts of the challenges posed by the aforementioned assumptions, many researchers have resorted to different approaches and modified models.Among these are the inclusion of jumps or stochastic parameters such as volatility in the price processes of option, Levy processes [3,4], the derivative of a no arbitrage determinant theorem for Liu's stock model in uncertain markets [5], models driven by uncertain processes for option pricing [6][7][8] and so on.
We remark here that the market value of a call option is a function of the underlying asset price, the exercise price, interest rate, expiration time, and the stock volatility [ ] . Despite these shortcomings, the Black-Scholes Model remains the hallmark of option pricing models for derivative security, and still proves very useful and vital both empirically and theoretically for the following reasons: the price of the option does not explicitly rely on the preferences of investors-risk -neutral valuation relationship (RNVR.Therefore, there is a greater need for better semi-analytical methods for solutions of such models resulting from stochastic differential equations (SDEs) and uncertain differential equations (UDEs).
In what follows, we will consider the classical celebrated Black-Scholes option pricing model: , we have: where For the derivation of Equation ( 1), Theorem 1 is applied (see [11] and the references therein for other necessary details).We remark here that Equation (1) holds for options whose underlying stock do not pay dividends provided that [ ] , and upon the satisfaction of the assumptions associated with the Black-Scholes model.
Many researchers and authors have attempted to obtain the solution of Equation (1) analytically and/or numerically, thereby adopting and using various direct and iterative methods, respectively.The classical Black-Scholes model is notable for its explicit closed form solution of European-style options (call and put options).On the other-hand, this is not generally true for non-European-style options whose closed form solutions do not exist, and even if they do exist, the techniques and approaches are complicated and even not easy to obtain using the conventional approaches, or methods, as such, Smeureanu and Fanache in [12], by means of several processors, via the finite difference method consider numerical solution of the Black-Scholes equation.
Cen and Le in [13] consider a numerical method based on central difference spatial discretization on a piecewise uniform mesh and an implicit time stepping technique for generalized Black-Scholes equation.In [14], Mosneagu and Dura apply numerical methods based on finite differences for solving Black-Scholes equation.Their intention is to create a general numerical scheme for different types of options.
Uddin, Ahmed and Bhowmik in [15], consider solution methods for the Black Scholes model with European options, by studying a weighted average method using different weights numerical approximations, and as such approximate the model using finite difference scheme.Algliardi, Popivanov and Slavova [16,17] consider the solution of the Black-Scholes equation via a Mellin transformation approach.Qiu and Lorenz in [18] study a modification of the Black-Scholes equation with regard to existence and uniqueness of solution to the Cauchy problem.
For the solution of fractional type Black-Scholes equation, Elbeleze, Kilicman and Taib combine the homotopy perturbation method (HPM), Sumudu transform, and He's polynomial [19], Kumar, Kumar and Singh in [20] apply a numerical algorithm [HPM], Ahmad et al. in [21] apply the Variational Iteration Method (VIM), while Kumar et al. [22], provide an analytical solution for the fractional Black-Scholes option pricing equation by homotopy perturbation method with coupling of the Laplace transform.
In solving both ordinary differential equations and partial differential equations of various forms including integro-differential equations encountered in finance [26], the fractional type ordinary differential equation in [27]; the relatively new semi-analytical method known as differential transform method (DTM) is shown to be effective, reliable and easier in application when compared to other semi-analytical methods [28,29], even when the results agree.
In this work, a modification of the DTM referred to as projected differential transform method (PDTM) is adopted and presented for the first, in solving the famous Black-Scholes equation in option valuation.
The remaining part of the paper is structured as follows: in Section 2, we give a brief introduction to DTM, PDTM and their fundamental properties; in Section 3, the PDTM is applied to solve some examples of the Black-Scholes equation, while in Section 4, we compare our results with regards to graphical interpretation, and in Section 5, we give concluding remarks.

The Differential Transformation Method (DTM) and Its Modification
In this section, we give a brief introduction to DTM, PDTM and their fundamental properties.
The differential transformation method (DTM) has been studied by many researchers and showed to be easier in terms of application when solving both linear and nonlinear differential equations as it converts said problems into their equivalents in algebraic recursive form.This is unlike other semi-analytical methods: ADM, VIM, HAM and so on, that require the determination of a successive term only by integrating a previous component [32,33].
Despite the many advantages of the DTM over other semi-analytical methods, some level of difficulty is still met when dealing mainly with the nonlinearity of differential equations and differential equations with variable coefficients.This again gives room for modification of the DTM in various forms by many authors and researchers [34,35].
In this work, a relatively new version of the modification referred to as the projected differential transform method (PDTM) will be applied to Black-Scholes equations for analytical and numerical solutions.

The Overview of the PDTM
Suppose ( ) , w x t is analytic at ( ) * * , x t in a domain D , then in considering the Taylor series of ( ) , w x t , regard is given to some variables v s t = instead of all the variables as seen in the classical DTM.Thus, the projected DTM of ( ) , w x t with respect to t at * t is defined and denoted by: and as such: where Equation ( 8) is referred to as the projected differential inverse transform (PDIT) of ( ) , W x h with respect to t .

Applications and Illustrative Examples
In this section, we solve some examples of the Black-Scholes equations with the proposed modified version of the DTM.

Discussion of Results
We present in this section Figures 1-4 to discuss our obtained results in comparison with their associated exact forms.

Concluding Remarks
In this study, a proposed computational method known as Projected Differential Transformation Method (PDTM) has been successfully applied, for the first time, to the Black-Scholes Equation for European Option Valuation.We solved some illustrative and numerical examples to test the efficiency of the proposed method.The results obtained converge faster to their associated exact solutions, even with less computation, without linearization or perturbation; showing that the method can also be used easily for approximate solutions in a direct form; these easily computed results represent the analytical values of the associated European call options, the same algorithm can be followed for European put options.
Finally, we remarked that the PDTM is very efficient, reliable; and faster in application (even without giving up accuracy) when compared with the classical DTM [20], ADM [17], HAM and HPM [12]; though the results via these methods are in strong agreement.Hence, it is strongly recommended for both linear and nonlinear stochastic differential equations (SDEs) encountered in financial mathematics and other areas of applied sciences.

2. 1 .
Analysis of a Two-Dimensional DTM Suppose ( ) , r x y a two-variable function is analytic at ( ) * * , x y in the Domain, D then, the differential transform of ( ) , r x y is defined and denoted as: