Abstract:Computational efficiency is essential for enhancing the accuracy and practicality of pricing complex financial derivatives. In this paper, we discuss Isogeometric Analysis (IGA) for valuing financial derivatives, modeled by two nonlinear Black-Scholes PDEs: the Leland model for European call with transaction costs and the AFV model for convertible bonds with default options. We compare the solutions of IGA with finite difference methods (FDM) and finite element methods (FEM). In particular, very accurate solutions can be numerically calculated on far less mesh (knots) than FDM or FEM, by using non-uniform knots and weighted cubic NURBS, which in turn reduces the computational time significantly.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to improve the computational efficiency and accuracy of financial derivatives pricing. Specifically, the paper explores the use of Isogeometric Analysis (IGA) to evaluate financial derivatives modeled by two nonlinear Black - Scholes partial differential equations (PDEs): the Leland model for European call options with transaction costs, and the AFV model for convertible bonds with default options. By comparing with the solutions of the Finite Difference Method (FDM) and the Finite Element Method (FEM), the paper shows that IGA can calculate numerical solutions with higher accuracy while reducing the number of grid points (i.e., the number of nodes), thereby significantly reducing the computation time.
### Main Problem Analysis
1. **Computational Efficiency and Accuracy**
- The paper points out that in the complex financial derivatives pricing, computational efficiency is crucial for improving the accuracy and practicality of pricing. Although traditional FDM and FEM methods perform well in some aspects, they have limitations when dealing with nonlinear PDEs. For example, the complexity of FDM increases when dealing with boundary conditions, and although FEM can provide higher - order approximation and effective boundary condition management, the use of Lagrange polynomial basis functions may lead to highly oscillating solutions.
- IGA can produce higher - order smooth solutions within the computational domain by using Non - Uniform Rational B - splines (NURBS) as basis functions, and NURBS basis functions are non - negative and non - oscillating, and can approximate non - smooth functions, such as the payoff function in option pricing, more accurately.
2. **Application of Nonlinear Models**
- **Leland Model**: This model takes into account the impact of transaction costs on European call options. The PDE model in the paper is:
\[
\frac{\partial V}{\partial t}+\frac{1}{2}\sigma^{2}S^{2}\left(1 + L_{e}\text{sgn}\left(\frac{\partial^{2}V}{\partial S^{2}}\right)\right)\frac{\partial^{2}V}{\partial S^{2}}+rS\frac{\partial V}{\partial S}-rV = 0
\]
where \( V = V(S,t) \) is the option price, \( S \) is the value of the underlying asset, \( r \) is the risk - free interest rate, \( \sigma \) is the volatility, and \( L_{e} \) is the Leland number, representing the impact of transaction costs.
- **AFV Model**: This model is used for pricing convertible bonds with default strategies. The PDE system in the paper is:
\[
\begin{aligned}
\frac{\partial U}{\partial t}+\frac{1}{2}\sigma^{2}S^{2}\frac{\partial^{2}U}{\partial S^{2}}+(r + p\eta)S\frac{\partial U}{\partial S}-(r + p)U + p\max(kS(1-\eta),\hat{R}B)&= 0,\\
\frac{\partial B}{\partial t}+\frac{1}{2}\sigma^{2}S^{2}\frac{\partial^{2}B}{\partial S^{2}}+(r + p\eta)S\frac{\partial B}{\partial S}-(r + p)B+\hat{R}pB&= 0,\\
\frac{\partial C}{\partial t}+\frac{1}{2}\sigma^{2}S^{2}\frac{\partial^{2}C}{\partial S^{2}}+(r + p\e