All submissions of the EM system will be redirected to Online Manuscript Submission System. Authors are requested to submit articles directly to Online Manuscript Submission System of respective journal.

Numerical Solution of a Linear Black-Scholes Model: A Comparative Overview

Md. Kazi Salah Uddin, Md. Noor-A-Alam Siddiki*, Md. Anowar Hossain

Department of Natural Science, Stamford University Bangladesh, Dhaka-1209, Bangladesh

*Corresponding Author:
Md. Noor-A-Alam Siddiki
Department of Natural Science, Stamford University Bangladesh, Dhaka-1209, Bangladesh
E-mail: siddiki@stamforduniversity.edu.bd

Received: 03/08/2015 Accepted: 02/09/2015 Published: 12/09/2015

Visit for more related articles at Research & Reviews: Journal of Statistics and Mathematical Sciences

Abstract

Black-Scholes equation is a well-known partial differential equation in financial mathematics. In this paper we try to solve the European options (Call and Put) using different numerical methods as well as analytical methods. We approximate the model using a Finite Element Method (FEM) followed by weighted average method using different weights for numerical approximations. We present the numerical result of semidiscrete and full discrete schemes for European Call option and Put option by Finite Difference Method and Finite Element Method. We also present the difference of these two methods. Finally, we investigate some linear algebra solvers to verify the superiority of the solvers.

Keywords

Black-Scholes model; Call and put options; Exact solution; Finite difference schemes, Finite Element Methods

Introduction

A powerful tool for valuation of equity options is the Black-Scholes mode [1,2]. This model is used for finding the prices of stocks. Company et al. [3] solved the modified Black-Scholes equation pricing option with discrete dividend. A delta-defining sequence of generalized Dirac-Delta function and the Mellin transformation are used toobtain an integral formula. Finally numerical quadrature approximation is used to approximate the solution. In some papers like [4] Mellin transformation is used. They were required neither variable transformation nor solving diffusion equation. Jodar et al. [4] found the solution of BS equation with a wide class of payoff functions that contains not only the Dirac delta type functions but also the ordinary payoff functions with discontinuities of their derivatives. Julia Ankudiova and Matthias Ehrhardt [5] solved nonlinear Black-Scholes equations numerically. They focused on various models relevant with the Black-Scholes equations with volatility depending on several factors. They also worked on the European Call option and American Call option analytically using transformation into convection diffusion equation with nonlinear term and the free boundary problem respectively. In our previous paper [6] we discussed about the analytical solution of Black-Scholes equation using Fourier Transformation method for European options. We formulated the Finite Difference Scheme and find the solutions of them. In this paper we discuss the solution with Finite Element Method and compare the result with the result obtained by Finite Difference Schemes.

Model Equation

The linear Black-Scholes equation [1,2] developed by Fischer Black and Myron Scholes in 1973 is

image (1)

where

V =V (S,t ), the pay − off function

S = S (t ), the stock price, with S = S (t ) ≥ 0,

t = time,

r = Risk − Free interest rate,

σ =Volatility Condition

And also t ∈ (0,T )

where T is time of maturity.

The terminal and boundary conditions [7] for both the European Call and Put options stated below.

European Call Option [7]

The solution to the Black-Scholes equation (1) is the value V (S,t ) of the European Call option on $ 0 ≤ S < ∞ , 0 ≤ t ≤ T . The boundary and terminal conditions are as follows

V (0,t) = 0 for 0 ≤ t ≤ T,

image (2)

image

European Put Option [7]

European Put option is the reciprocal of the European Call option and the boundary and terminal conditions are

image

V (S,t )→ 0 as S →∞, (3)

image

Transformation

The model problem stated in (1) is a backward type. This type is little bit difficult to solve. To solve the problem in (1) with the conditions stated in (2) and (3) we need to make the model in forward type. In this regard, we have the following transformations.

Let

image

image

And

image

image

image

image

image

inserting these derivatives in equation (1) we have

image

implies

image (4)

Let

image

(4) implies

image

Now let

image

so that

image

image

image

image

image

image

inserting these into equation(5) and dividing by image we get

image

Implies

image

image

image

And the initial an boundary conditions for the European Call and Put options are respectively

image (7)

and

image

image(8)

image

Thus the Black-Scholes equation reduced to a heat diffusion equation

Numerical Approximation of Transformed Linear Black-Scholes Model

Now we solve the problems numerically. We use the Finite Element Method (FEM) to solve the problems related to the differential equation (6). Finally back substitution of the coordinate transformation gives the solution of the problems related to the differential equation (1).

We have the model

image (9)

and the initial and boundary conditions for call option are

image

image (10)

image

The weak form of the governing equation is

image (11)

image

Discretizing image spatially, we have

image (12)

where imageare given shape functions, and image are unknown, and n is the ordinal number of nodes.

Substituting (12) into (11), we get the weak semidiscretized equation

image(13)

Let imagedenote the so-called mass and stiffness matrices, respectively, defined by:

image (14)

image (15)

Then (13) can be expressed as:

image (16)

whereimage is a vector function with the componentsimage

After performing the integral in (14) and (15) for the linear shape functions, the mass and the stiffness matrices have the following form

image

where h is the length of the spatial approximation.

Now we would like to discrete the equation (16) with respect to time. One may start with a simple scheme.

One of the trivial choice is to use the forward Euler scheme. Firstly we discrete (16) explicitly and we have

image

image

image

image

image

image

The difficulty of using the scheme is that it needs very little step size to converge , as a result the scheme is a slow one, and is not of interest in this advance study.

We want a fast and efficient scheme, so we want larger time stepping, and interested in using implicit techniques. We dis crete (16) implicitly and have

image

image

image

image (18)

which is a system of linear equations with unknowns image . The advantage of using (18)

is that the scheme is unconditionally stable . Equation (18) accuracy of order image . It is faster than the explicit Euler scheme since (18) allows us to use large time steps.

We use an weighted average method to discrete (16) with weight and we have

image

image

image

image(18)

which is a system of linear equations with unknowns image. The advantage of using (18)

is that the scheme is unconditionally stable . Equation (18) accuracy of order O(k ) . It is faster than the explicit Euler scheme since (18) allows us to use large time steps.

We use an weighted average method to discrete (16) with weight and imagewe have

image

image

image

image (19)

This system is also a linear one with unknowns image where image varies from 0 to 1. This method turns to the explicit method when imagei.e., equations (17) and (19) are same and implicit method when image,i.e., equations (18) and (19)are same. For image the scheme is conditionally stable and unconditionally stable for image

The order of the accuracy of the scheme is image

Resultas, Discusión and Conclusión

In this section we have presented the results by various methods. We have solved the model analytically [8] by the method of Fourier Transformation. In Figure 1 we placed the analytic solution of two options (Call Option and Put Option). To solve the model numerically we have applied [8] Finite difference methods (FDM) and have shown the result of the two options in Figure 2. Our interest in this paper was in the methods of Finite Elements (FEM) [9] .Firstly, we have discretized the model (6) spatially in the section (4). Then we have used various one step Euler’s time integrations to discretize the system of linear equations obtained by semi-discretization. The results have been presented in Figure 3. We have tried to show comparison between the methods (FDM and FEM) in Figure 4.

statistics-and-mathematical-sciences-Analytic-solutions

Figure 1: Analytic solutions

statistics-and-mathematical-sciences-Numerical

Figure 2: Numerical Solutions by Finite Difference Method.

statistics-and-mathematical-sciences-Solutions

Figure 3: Numerical Solutions by Finite Element Method.

statistics-and-mathematical-sciences-Comparison

Figure 4: Comparison of Finite Difference Method and Finite Element Method.

The system of linear equations (19) generated by the discretization of the Black-Scholes model can be solved by many conventional processes. For a large scale linear system, scientists rarely use direct methods as they are computationally costly. Here, in this section, it is our motivation to solve the system of equation (19) using various iterative techniques. Here we first investigate which linear solver converges swiftly. To that end, we consider Jacobi iterative method, Gauss-Seidel iterative method and successive over relaxation method to start with. In terms of matrices, the Jacobi method can be expressed as

image

Gauss-Seidel method

image

and the SOR algorithm can be written as

image

where in each case the matrices D, − L , and −U represent the diagonal, strictly lower triangular, and strictly upper triangular parts of A, respectively.

We investigate Preconditioned Conjugate Gradient (PCG) Method and Generalized Minimal Residual (GMRES) Method with a diagonal preconditioning [10]. Here for all computations we consider image. The results are presented with different weights δ. Observing Figure 5, we notice that Preconditioned Conjugate Gradient (PCG) Method performs the best.

statistics-and-mathematical-sciences-linear

Figure 5: Time comparison of different linear algebra solvers.

References