Anbo Le

A robust and accurate finite difference method for a generalized Black-Scholes equation

In this paper we present a numerical method for a generalized Black-Scholes equation, which is used for option pricing. The method is based on a central difference spatial discretization on a piecewise uniform mesh and an implicit time stepping technique. Our scheme is stable for arbitrary volatility and arbitrary interest rate, and is second-order convergent with respect to the spatial variable. Furthermore, the present paper efficiently treats the singularities of the non-smooth payoff function. Numerical results support the theoretical results.

A second-order hybrid finite difference scheme for a system of singularly perturbed initial value problems

A system of coupled singularly perturbed initial value problems with two small parameters is considered. The leading term of each equation is multiplied by a small positive parameter, but these parameters may have different magnitudes. The solution of the system has boundary layers that overlap and interact. The structure of these layers is analyzed, and this leads to the construction of a piecewise-uniform mesh that is a variant of the usual Shishkin mesh. On this mesh a hybrid finite difference scheme is proved to be almost second-order accurate, uniformly in both small parameters. Numerical results supporting the theory are presented.

Finite difference scheme with a moving mesh for pricing Asian options

In this paper we propose a stable numerical method for pricing Asian call options, which is based on a central difference scheme with a moving mesh in the spatial discretization and the Rannacher time stepping scheme in the time discretization. At each time mesh point we make a piecewise uniform mesh to discretise the space interval, which ensures that the matrix associated with the discrete operator is an M-matrix. Hence the spatial discretization scheme is maximum-norm stable for arbitrary volatility and arbitrary interest rate. We show that the scheme is second-order convergent with respect to both time and spatial variables. Numerical results support the theoretical results.

Exponential Time Integration and Second-Order Difference Scheme for a Generalized Black-Scholes Equation

We apply an exponential time integration scheme combined with a central difference scheme on a piecewise uniform mesh with respect to the spatial variable to evaluate a generalized Black-Scholes equation. We show that the scheme is second-order convergent for both time and spatial variables. It is proved that the scheme is unconditionally stable. Numerical results support the theoretical results.

A posteriori error analysis for a fractional differential equation

ABSTRACT Numerical treatment for a fractional differential equation (FDE) is proposed and analysed. The solution of the FDE may be singular near certain domain boundaries, which leads to numerical difficulty. We apply the upwind finite difference method to the FDE. The stability properties and a posteriori error analysis for the discrete scheme are given. Then, a posteriori adapted mesh based on a posteriori error analysis is established by equidistributing arc-length monitor function. Numerical experiments illustrate that the upwind finite difference method on a posteriori adapted mesh is more accurate than the method on uniform mesh.

A hybrid finite difference scheme for pricing Asian options

In this paper we apply a hybrid finite difference scheme to evaluate the prices of Asian call options with fixed strike price. We use the Crank-Nicolson method to discretize the time variable and a hybrid finite difference scheme to discretize the spatial variable. The hybrid difference scheme uses the central difference approximation whenever the mesh points are sufficiently away from the left-hand side of the domain to ensure the stability of the scheme; otherwise a midpoint upwind difference scheme is used. The matrix associated with the discrete operator is an M-matrix, which ensures that the spatial discretization scheme is maximum-norm stable. It is proved that the scheme is second-order convergent with respect to both time and spatial variables. Numerical experiments support these theoretical results.

A Finite Difference Scheme for Pricing American Put Options under Kou's Jump-Diffusion Model

We present a stable finite difference scheme on a piecewise uniform mesh along with a penalty method for pricing American put options under Kous jump-diffusion model. By adding a penalty term, the partial integrodifferential complementarity problem arising from pricing American put options under Kous jump-diffusion model is transformed into a nonlinear parabolic integro-differential equation. Then a finite difference scheme is proposed to solve the penalized integrodifferential equation, which combines a central difference scheme on a piecewise uniform mesh with respect to the spatial variable with an implicit-explicit time stepping technique. This leads to the solution of problems with a tridiagonal M-matrix. It is proved that the difference scheme satisfies the early exercise constraint. Furthermore, it is proved that the scheme is oscillation-free and is second-order convergent with respect to the spatial variable. The numerical results support the theoretical results.

An Alternating-Direction Implicit Difference Scheme for Pricing Asian Options

We propose a fast and stable numerical method to evaluate two-dimensional partial differential equation (PDE) for pricing arithmetic average Asian options. The numerical method is deduced by combining an alternating-direction technique and the central difference scheme on a piecewise uniform mesh. The numerical scheme is stable in the maximum norm, which is true for arbitrary volatility and arbitrary interest rate. It is proved that the scheme is second-order convergent with respect to the asset price. Numerical results support the theoretical results.

A robust numerical method for a fractional differential equation

Abstract This paper is devoted to giving a rigorous numerical analysis for a fractional differential equation with order α  ∈ (0, 1). First the fractional differential equation is transformed into an equivalent Volterra integral equation of the second kind with a weakly singular kernel. Based on the a priori information about the exact solution, an integral discretization scheme on an a priori chosen adapted mesh is proposed. By applying the truncation error estimate techniques and a discrete analogue of Gronwall’s inequality, it is proved that the numerical method is first-order convergent in the discrete maximum norm. Numerical results indicate that this method is more accurate and robust than finite difference methods when α is close to 0.

A robust upwind difference scheme for pricing perpetual American put options under stochastic volatility

In this paper, we present an upwind difference scheme for the valuation of perpetual American put options, using Hestons stochastic volatility model. The matrix associated with the discrete operator is an M-matrix, which ensure that the scheme is stable. We apply the maximum principle to the discrete linear complementarity problem in two mesh sets and derive the error estimates. Numerical results support the theoretical results.