Abstract:In this paper, an integral equation representation for the early exercise boundary of an American option contract is considered. Thus far, a number of different techniques have been proposed in the literature to obtain a variety of integral equation forms for the early exercise boundary, all starting from the Black-Scholes partial differential equation. We first present a coherent categorization of exiting integral equation methodologies in the American option pricing literature. In the reminder and based on the fact that the early exercise boundary satisfies a fully nonlinear weakly singular non-standard Volterra integral equation, we propose a product integration approach based on linear barycentric rational interpolation to solve the problem. The price of the option will then be computed using the obtained approximation of the early exercise boundary and a barycentric rational quadrature. The convergence of the approximation scheme will also be analyzed. Finally, some numerical experiments based on the introduced method are presented and compared to some exiting approaches.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is the approximate calculation problem of the early exercise boundary in American option pricing. Specifically, the article proposes a product integration method based on linear barycentric rational interpolation to solve this boundary and then calculate the option price.
### Problem Background
American options allow the holder to exercise the option at any time before the expiration date, which makes their pricing problem much more complicated than that of European options. In particular, determining the early exercise boundary is a key issue in American option pricing. The early exercise boundary determines under what circumstances the option should be exercised to maximize the return.
### Existing Methods
Many existing techniques have proposed different forms of integral equations to describe the early exercise boundary. Most of these methods are derived from the Black - Scholes partial differential equation. However, these methods face challenges when dealing with non - standard Volterra integral equations, especially in dealing with weak singularities and numerical stability.
### Paper Contributions
The main contributions of this paper include:
1. **Classification and Review**: First, a systematic classification and review of the integral equation methods for the early exercise boundary in the existing literature are carried out.
2. **Product Integration Method**: A new product integration method based on linear barycentric rational interpolation is proposed to solve the fully nonlinear weakly singular non - standard Volterra integral equation.
3. **Numerical Experiments**: The effectiveness of the proposed method is verified through numerical experiments and compared with existing methods.
### Formula Representation
The key formulas involved in the paper include:
- The weakly singular Volterra integral equation of the early exercise boundary:
\[
K - B(t)=p_E(t, B(t))+\int_{0}^{t}\left[rK e^{-r(t - s)}\aleph(-d_2(B(t), t - s, B(s)))-\delta B(t)e^{-\delta(t - s)}\aleph(-d_1(B(t), t - s, B(s)))\right]ds
\]
where
\[
p_E(t, S)=K e^{-rt}\aleph(-d_2(S, t, K))-B(t)e^{-\delta t}\aleph(-d_1(S, t, K))
\]
\[
d_1(x, t, y)=\frac{\log\left(\frac{x}{y}\right)+(r - \delta+\frac{\sigma^{2}}{2})t}{\sigma\sqrt{t}}, \quad d_2(x, t, y)=d_1(x, t, y)-\sigma\sqrt{t}
\]
- The discretized form of the product integration method:
\[
u_i = g(t_i, u_i)+\sum_{j = 0}^{i}w_{i,j}k_1(t_i, t_j, u_i, u_j)+\sum_{j = 0}^{i}\omega_jk_2(t_i, t_j, u_i, u_j)
\]
where $w_{i,j}=\int_{0}^{t_i}\frac{L_j(s)}{\sqrt{t_i - s}}ds$ and $\omega_j$ are quadrature weights.
Through these methods, the paper aims to provide an efficient and accurate numerical approximation scheme for the early exercise boundary, thereby providing a more reliable tool for American option pricing.