Application of Elzaki Transform Method to Market Volatility Using the Black-Scholes Model ()
1. Introduction
Volatility has to do with the degree of variation in prices of financial instruments over a period of time. The market expectation of future price fluctuations can be measured and expressed as a financial model. The Black-Scholes model is a partial differential equation responsible for the mathematical analysis of price evolution in a European call ( [1] [2] ). Basically, the Black-Scholes Model (B-SM) as the foundation for option pricing can be used to explore a variety of options in a European call. It is efficient to compute the price of an option correctly [3] . It can also be used to reveal different observable factors such as underlying price options, determine the range of price movements, rates and other factors not simultaneously observable and underlying fear factors according to Hicks [4] . For a European call depending on stock paying with no dividends, the B-SM is given as [5]
(1)
with
as
. (2)
where;
· V is the price of the option;
·
is the volatility of the stock;
· S is the stock price; and
· r is the rate-free interest rate.
Biazar et al. [6] carried out a detailed study of Equation (1) and suggested a straightforward version which is rewritten as
, (3)
for proper financial interpretation. Equation (3) is referred to as the standard B-SM.
FOKKER-PLANCK EQUATION ( [7] [8] )
Let us consider the general form of a nonlinear FPE
(4)
With initial conditions as
,
. Such that
is the unknown distribution function in Equation (4). Also, a special case of equation (4) is the FPE derived from Plasma Physics having the form
(5)
in Equation (5) is the probability density function. Expressing
and substituting for it, then Equation (5) is reduced to a linear Brownian motion of FPE given as
with initial conditions
Or as
(6)
using the Mellin transformation [9] .
S is the the exact price that undergoes Brownian motion,
is the call option price, V is the volatility and r is the interest rate. Equation (6) is more useful if the problem is to determine the call option. However, B-SM when solved provides a strategy for determining not just the prices of underlying security, current time and time of expiration, and return of the free-risk assets but also eliminating the risk involved in the whole process.
Assumptions of the Black-Scholes Model
· The model B-SM computes the present value of a call option without dividends.
· The model requires the value of the standard deviation. It can be calculated from the variance by taking the square root.
· The hedge ratio is the ratio of the expected stock price at expiration to the current stock price at the onset.
· Free risk assets, underlying prices, strike prices and expiration can be used to hedge positions on the option.
Equation (1) can both be solved analytically and numerically. The few available analytic methods are complex and difficult to handle requiring linearization, perturbation, or weak assumptions. Therefore, several numerical methods have been developed and implemented over the years for seeking the solution of the Black-Scholes equation. The Finite Difference Method (FDM) has often been the most privileged numerical scheme used by many researchers in recent times to seek the approximate solutions of the B-SM. For instance, Andallah and Anwar [5] used the FDM for the numerical solution of B-SM. Dura and Mosneagu [10] also considered the numerical approximation of the B-SM via the FDM. Duffy [11] employed the FDM in financial engineering through the B-SM. A unitary transformation has been applied to the classical B-SM to obtain the quantum effects associated with the market fear factor especially that causing the increase in volatility rate [12] . Characterizing and determining volatility rate has also been the focus of many researchers [13] . There is also the non-commutative B-SM presented in Accardi and Boukas [14] that was transformed into an integral equation. Wei [15] listed the three benchmarks of a stochastic system for B-SM to test the accuracy of Fokker-Planck Equation (FPE) via time-dependent methods of FDM, Finite Volume and Finite Element. The B-SM of option pricing can also be regarded as a reaction-diffusion equation that is entirely based on stochastic analysis. Here, the B-SM is first transformed into a simple FPE and then analyzed. Several research has been carried out in the computational sense for the numerical solution of this derived FPE ( [9] [15] [16] ).
The empirical performance of B-SM has been examined in [17] . However, the link between B-SM and FPE has been investigated and analyzed by many researchers. These include [18] ; they even introduced Fractional FPE in their work by deriving the fractional FPE order governing the dynamics of the equation and then determine the Black-Scholes differential equation that involves the stock asset and fair prices of the European option. Before this, others have used fractional to analyze stock exchange market dynamics ( [19] [20] [21] ). Consequently, once the market prices have been observed for the option, the B-SM can be inverted to determine the volatility. In addition, the B-SM when solved provides a strategy for eliminating the risk and determining the prices of underlying security ( [22] [23] ).
In recent years, the focus on B-SM has been on reviews in terms of redefinition ( [24] [25] ), Computations [26] , Applications ( [27] [28] ) and solution methods [29] and sensitivity analysis [30] .
This paper is motivated to seek the numerical solution of the B-SM using the Elzaki transform method (ETM). This method is necessitated by its simplicity, less computational rigor, and rapid rate of convergence as will be observed in section 3.
2. Elzaki Transform Method
The set [31]
, (7)
define by
, (8)
is called the Elzaki transform of the function
.
2.1. Properties of Elzaki Transforms Method
1)
.
2)
.
3)
, m is the order of the highest derivative.
4)
.
5)
.
2.2. Elzaki Transform Method for Black-Scholes Model
Here, we apply the ETM to solve B-S equation of the form (1).
Applying the Elzaki transform on both sides of Equation (2), we have
(9)
Using 2.1(iii), we have that
(10)
Having the inverse on both sides if Equation (10), we have
(11)
Thus the approximate solution is given as
Thus, Equation (11) becomes
(12)
By comparing both sides of Equation (12), we obtain the recurrence relations
(13)
(14)
Hence,
for
are calculated using the Equations (13) and (14).
Finally, the required approximate solution of the Black-Scholes equation becomes
(15)
Equation (15) is the determining index for cash flow, risk profile for buying and selling underlying asset and predicting future price movement.
3. Numerical Illustrations
Here, we implement the ETM on the B-SM using Maple 18 software. Results obtained is presented as a power series.
Thus, using the ETM scheme (13) and (14) on
,
with the parameters
,
,
and
, (Andallah and Anwar, 2018), we obtain the following approximations:
Thus, the required computed solution is
(16)
Using the values of the call option on Equation (16) over the range 0 ≤ S ≤ 100 for t = 0 to 1, the striking price is illustrated in Figure 1 produced by MAPLE 18. The approximate values of the model t = 0 to 1
Example 2
Let
(17)
be an approximate solution of (6) with
being constant parameters, and
, Mamadu-Njoseh polynomials.
Figure 1. Numerical simulation for European call at different time steps with
,
,
and
.
Now, substituting (17) into (6), we get
(18)
Evaluating (18) at the conditions
we obtain the following system of equations:
(19)
Thus, solving the above equations for
, using the estimates
,
,
,
,
, we obtain Figure 2 below.
Figure 2. Numerical simulation for Equation (6) at different time steps with
,
,
,
and
.
Substituting the above into (17) yields the approximate solution to (6). Further simulation gives.
4. Conclusion
We have considered the numerical solution of the B-SM using the Elzaki transform method. Results obtained showed that the ETM is very reliable for simulating the call depending on stock paying with no dividends. The numerical illustration shows that the method is accurate for solving the FPE, hence for B-S equation. Also, the results agreed with those found in the literature ( [5] [31] ) for comparison. B-SM in its numerous modifications and expansion [31] can be used for hedging and risk mitigation, put and call options, Forex options, cash call, assets with continuous yielding dividends, stock options etc (DF, 2022) [32] .