
This study presents a numerical method for accurately computing the option values and Greeks of equity-linked securities (ELS) near early redemption dates. The Black–Scholes (BS) equation is solved using the finite difference method (FDM), and a Dirichlet boundary condition is applied at strike prices instead of directly replacing option values above the strike price with predefined option prices. This approach improves the accuracy of option pricing, particularly in the presence of early redemption structures. The proposed method is demonstrated to be effective in computing Greeks, which are crucial for risk management and hedging strategies in ELS markets. The computational tests validate the reliability of the method in capturing the sensitivities of ELS prices to various market factors.
Citation: Yunjae Nam, Changwoo Yoo, Hyundong Kim, Jaewon Hong, Minjoon Bang, Junseok Kim. Accurate computation of Greeks for equity-linked security (ELS) near early redemption dates[J]. Quantitative Finance and Economics, 2025, 9(2): 300-316. doi: 10.3934/QFE.2025010
[1] | Mustafa Tevfik Kartal, Özer Depren, Serpil Kılıç Depren . The determinants of main stock exchange index changes in emerging countries: evidence from Turkey in COVID-19 pandemic age. Quantitative Finance and Economics, 2020, 4(4): 526-541. doi: 10.3934/QFE.2020025 |
[2] | Dumisani Pamba, Sophia Mukorera, Peter Moores-Pitt . The asymmetric effects of cross-border equity flow volatility on equity market returns in SANEK countries. Quantitative Finance and Economics, 2025, 9(1): 40-75. doi: 10.3934/QFE.2025002 |
[3] | Muath Asmar, Susanne Trimbath . Regulatory reform and trade settlement failures in USA equity markets: Does regulatory reform matter?. Quantitative Finance and Economics, 2022, 6(4): 537-552. doi: 10.3934/QFE.2022023 |
[4] | Mehmet F. Dicle . US Implied Volatility as A predictor of International Returns. Quantitative Finance and Economics, 2017, 1(4): 388-402. doi: 10.3934/QFE.2017.4.388 |
[5] | Manuela Larguinho, José Carlos Dias, Carlos A. Braumann . Pricing and hedging bond options and sinking-fund bonds under the CIR model. Quantitative Finance and Economics, 2022, 6(1): 1-34. doi: 10.3934/QFE.2022001 |
[6] | Andreas Tsalas, Platon Monokroussos . PIIGS in and out of sync: the changing face of financial business cycle synchronization in Europe. Quantitative Finance and Economics, 2018, 2(1): 261-278. doi: 10.3934/QFE.2018.1.261 |
[7] | Anwen Yin . Equity premium prediction: keep it sophisticatedly simple. Quantitative Finance and Economics, 2021, 5(2): 264-286. doi: 10.3934/QFE.2021012 |
[8] | Yonca Erdem Demirtaş, Neslihan Fidan Keçeci . The efficiency of private pension companies using dynamic data envelopment analysis. Quantitative Finance and Economics, 2020, 4(2): 204-219. doi: 10.3934/QFE.2020009 |
[9] | Abdulsalam Abidemi Sikiru, Afees A. Salisu . Hedging with financial innovations in the Asia-Pacific markets during the COVID-19 pandemic: the role of precious metals. Quantitative Finance and Economics, 2021, 5(2): 352-372. doi: 10.3934/QFE.2021016 |
[10] | Takashi Kanamura . Supply-side perspective for carbon pricing. Quantitative Finance and Economics, 2019, 3(1): 109-123. doi: 10.3934/QFE.2019.1.109 |
This study presents a numerical method for accurately computing the option values and Greeks of equity-linked securities (ELS) near early redemption dates. The Black–Scholes (BS) equation is solved using the finite difference method (FDM), and a Dirichlet boundary condition is applied at strike prices instead of directly replacing option values above the strike price with predefined option prices. This approach improves the accuracy of option pricing, particularly in the presence of early redemption structures. The proposed method is demonstrated to be effective in computing Greeks, which are crucial for risk management and hedging strategies in ELS markets. The computational tests validate the reliability of the method in capturing the sensitivities of ELS prices to various market factors.
The step-down equity-linked security (ELS) is a structured investment instrument where the coupon or interest rate paid to investors increases over time based on predefined conditions, frequently tied to the performance of a specific equity index or stock (Jo and Kim, 2013).
Since the fourth quarter of 2022, Korean financial regulators have been directing Hang Seng China Enterprises Index (HSCEI, H-index)-based ELS sellers (5 banks, 7 securities companies) to establish a customer response system to prepare for potential investor losses following the sharp decline in the H-index. The problematic product is mainly ELS issued since early 2021, when the H-index was at its peak, and given the recent trend of the H-index, investor losses are likely to occur from early 2024, when the ELS will mature.
It is necessary to compute the Greeks of ELS and proceed with the overall interpretation so that we can prepare for these risks. Greeks computations, which refer to the measures of risk derived from options pricing models, are important in step-down ELS derivatives. Through Greek values, we can interpret the sensitivity of step-down ELS to fluctuations in factors like the underlying asset price, maturity, volatility, and interest rate. When designing step-down ELS, issuers use Greek computations to determine the optimal structure of the product to meet the desired risk-return profile. By incorporating Greeks into pricing models, market participants can estimate the fair value of the product under different market scenarios and assess its attractiveness relative to alternative investments.
Meanwhile, ongoing numerical investigations have been conducted and continue to be actively pursued to price the fair value of an ELS derivative with a step-down structure composed of multiple underlying assets. Wu et al. (2023) studied step-down ELS options under the fractional BS model, which incorporates long-term stock correlations. They revealed that higher Hurst exponents increase repayment returns in multi-asset ELS. Cui et al. (2024) proposed the step-down ELS evaluation by applying a unified pricing framework with continuous-time Markov chain approximation and employed stochastic gradient descent to minimize the stochastic local volatility of hedging loss under various asset models. Lee et al. (2024) studied evaluating step-down ELS by applying logistic regression with the Monte Carlo method from Hwang et al. (2023) to predict exit probabilities in a multi-dimensional Brownian bridge. This approach demonstrates the efficiency of numerically estimating the exit and co-exit probabilities, which enables effective pricing of complex structured products.
In this study, we introduce a computational method for the accurate computation of the option prices and Greeks for ELS near early redemption dates. We consider the price of an ELS option with one asset. The parameters are set as follows: maturity time T, the computational domain Ω=(0,300), strike prices [K1,K2,K3,K4,K5,K6], coupon rates [c1,c2,c3,c4,c5,c6], knock-in barrier D, dummy rate d=F(1+c6) with face value F, and stock price S0. The strike prices are arranged as K1≥K2≥K3≥K4≥K5≥K6, and the coupon rates set as c1≤c2≤c3≤c4≤c5≤c6 for each early redemption date T1≤T2≤T3≤T4≤T5≤T6.
Two option prices v(S,t) and u(S,t) can be defined through the equation solved numerically, corresponding to whether the knock-in barrier is hit or not. If the knock-in barrier is touched, the payoff of ELS is as follows:
v(S,Ti)={F(1+c6),if S≥K6,FS/S0,otherwise. | (1) |
In more detail, the payoff scenario when the knock-in barrier is touched is shown in Figure 1 (①, ②, ③, ④, ⑤, ⑥, ⑧, ⑨). In contrast, without the knock-in barrier, the payoff of ELS is given by
u(S,Ti)={F(1+c6),if S≥K6,F(1+d),if D<S<K6,FS/S0,otherwise. | (2) |
A detailed illustration of the payoff scenario when the knock-in barrier is not breached is presented in Figure 1 (①, ②, ③, ④, ⑤, ⑥, ⑦, ⑧).
To evaluate the price of derivatives, the no-arbitrage theory in continuous time provides a fundamental framework. Under this framework, the price of the underlying asset St is assumed to follow a geometric Brownian motion (GBM), which is modeled by the following stochastic differential equation:
dSt=rStdt+σStdWt, |
where r represents the drift rate, σ is the volatility, and Wt is a standard Wiener process. This assumption implies that the logarithm of the asset price follows a GBM with drift, which serves as the foundation for the derivation of the BS equation. Black and Scholes (1973) formulated a partial differential equation for option pricing based on the no-arbitrage principle and the concept of a self-financing portfolio. By constructing a risk-free portfolio composed of the option and the underlying asset, the BS equation is obtained as follows:
∂v(S,t)∂t+(σS)22∂2v(S,t)∂S2+rS∂v(S,t)∂S−rv(S,t)=0, | (3) |
where v(S,t) represents the option price at the underlying asset price S and time t and r is the risk-free interest rate. The BS equation is solved by imposing a terminal condition at maturity T, v(S,T), along with appropriate boundary conditions.
The BS equation is a widely used method for evaluating derivatives, and various numerical methods have been developed to evaluate increasingly complex products. Liu et al. (2024) proposed a novel numerical method that employs the radial basis function-finite difference solver to approximate derivatives and solve multi dimensional option pricing problems. Roul and Goura (2020) introduced a novel numerical scheme for option pricing, utilizing the Crank–Nicolson and sextic B-spline collocation methods to solve the generalized BS equation. Tao et al. (2023) conducted a study that proposed an Asian option pricing formula based on sub fractional Brownian motion and the sub fractional Vasicek model, analyzing its impact using the BS equation and a novel numerical scheme.
The ELS price and its Greeks are calculated through the finite difference method (FDM). The proposed algorithm solves the BS equation at early redemption dates using the Dirichlet boundary condition at strike prices instead of overwriting the option values on greater than or equal to strike prices by the given option prices.
This article is structured as follows. Section 2 describes the computational solution with FDM. Section 3 provides the values and surfaces of Greeks. In Section 4, we numerically evaluate the prices and Greeks of ELS products for different payoffs. Finally, Section 5 concludes the paper.
Let τ=T−t. Then, Equation (3) becomes
∂v(S,τ)∂τ=(σS)22∂2v(S,τ)∂S2+rS∂v(S,τ)∂S−rv(S,τ) | (4) |
with an initial condition v(S,0). To numerically solve Equation (4), we restrict the infinite underlying domain for S into a finite domain, [0,Smax]. Equation (4) is discretized using an implicit FDM with a nonuniform grid (Lee et al., 2023):
vn+1i−vniΔτ=(σSi)22(2vn+1i−1hi−1(hi−1+hi)−2vn+1ihi−1hi+2vn+1i+1hi(hi−1+hi))+rSi(−hivn+1i−1hi−1(hi−1+hi)+(hi−hi−1)vn+1ihi−1hi+hi−1vn+1i+1hi(hi−1+hi))−rvn+1i, | (5) |
where vni≡v(Si,τn), Si+1=Si+hi, S1=0, SNx=Smax, τn+1=τn+Δτ, and τ1=0. Here, hi and Δτ are the step sizes of the underlying asset and time, respectively (Liu et al., 2024). We use a linear boundary condition; therefore, it is convenient to define SNx+1=2SNx−SNx−1 so that the numerical second derivative at SNx is always zero and the numerical solution is linear at the right boundary. Figure 2 displays a schematic illustration of the underlying asset S.
Equation (5) can be rewritten as follows:
αivn+1i−1+βivn+1i+γivn+1i+1=fi,i=2,…,Nx, | (6) |
where
αi=−(σSi)2+rSihihi−1(hi−1+hi),βi=1Δτ+(σSi)2−rSi(hi−hi−1)hi−1hi+r, | (7) |
γi=−(σSi)2−rSihi−1hi(hi−1+hi),fi=vniΔτ. | (8) |
When early redemption is not possible, the following boundary conditions are satisfied. From the zero Dirichlet boundary condition at i=1, we have vn+11=0. At i=Nx, we apply the linear boundary condition, vn+1Nx+1=2vn+1Nx−vn+1Nx−1 and get
(αNx−γNx)vn+1Nx−1+(βNx+2γNx)vn+1Nx=fNx. | (9) |
Then, the solution vector vn+12:Nx=[vn+12vn+13⋯vn+1Nx]T can be obtained by solving the tridiagonal system
Avn+12:Nx=f2:Nx, | (10) |
where A is a tridiagonal matrix formed from Equation (6) with the zero Dirichlet at i=1 and linear boundary conditions at i=Nx, That is,
A2:Nx=(β2γ200⋯000α3β3γ30⋯0000α4β4γ4⋯000⋮⋮⋮⋮⋱⋮⋮⋮0000⋯αNx−1βNx−1γNx−10000⋯0αNx−γNxβNx+2γNx) |
and
f2:Nx=(f2f3f4⋮fNx−1fNx). |
Let SKj be the price of early repayment and cj be the coupon rate of early repayment for j=5,⋯,1. When early redemption occurs at the evaluation time, v has the Dirichlet boundary conditions at i=1, we have vn+11=0 and at i=Kj, vn+1Kj=F(1+cj). Then, the solution vector vn+12:Kj−1=[vn+12vn+13⋯vn+1Kj−1]T can be determined by solving the tridiagonal system
Avn+12:Kj−1=f2:Kj−1, | (11) |
where A is a tridiagonal matrix constructed from Equation (6) with the zero Dirichlet boundary condition at i=1 and the Dirichlet boundary condition at i=Kj, That is,
A2:Kj−1=(β2γ200⋯000α3β3γ30⋯0000α4β4γ4⋯000⋮⋮⋮⋮⋱⋮⋮⋮0000⋯αKj−2βKj−2γKj−20000⋯0αKj−1βKj−1) |
and
f2:Kj−1=(f2f3f4⋮fKj−1fKj−γKj−1vn+1Kj). |
We utilize the Thomas algorithm (Thomas) to solve Equations (10) and (11), which inverts the tridiagonal matrix directly.
First, we calculate v(S,τ). Subsequently, Figure 3 illustrates the option values v at both the pre-redemption and the redemption date. These values are then used to compute u(S,τ) by applying Dirichlet boundary conditions at the knock-in barrier. From the Dirichlet boundary condition at i=D and i=Kj, we obtain un+1D=vn+1D and un+1Kj=F(1+cj), as shown in Figure 4.
βD+1un+1D+1+γD+1un+1D+2=gD+1−αD+1un+1D, | (12) |
αKj−1un+1Kj−2+βKj−1un+1Kj−1=gKj−1−γKj−1un+1Kj. | (13) |
Then, the solution vector un+1D+1:Kj−1=[un+1D+1un+1D+2⋯un+1Kj−1]T can be obtained by resolving the tridiagonal system
Bun+1D+1:Kj−1=gD+1:Kj−1, | (14) |
where A is a tridiagonal matrix formed from similar to Equation (6) with the Dirichlet at i=D and i=Kj, That is,
AD+1:Kj−1=(βD+1γD+100⋯000αD+2βD+2γD+20⋯0000αD+3βD+3γD+3⋯000⋮⋮⋮⋮⋱⋮⋮⋮0000⋯αKj−2βKj−2γKj−20000⋯0αKj−1βKj−1) and fD+1:Kj−1=(fD+1−αD+1un+1DfD+2fD+2⋮fKj−2fKj−1−γKj−1un+1Kj). |
We also utilize the Thomas algorithm (Thomas) to solve Equation (11), which inverts the tridiagonal matrix (Kwak et al., 2023). Then, the option value u at the redemption date can be obtained, as illustrated in Figure 5.
In this section, we compute the values of the Greeks for ELS and visualize their surfaces, which include delta (Δ), gamma (Γ), theta (Θ), vega (ν), and rho (ρ). Unless otherwise stated, this section uses the parameters summarized in Table 1.
Parameter | Value |
Maturity time, T | 3 |
Repayment times, [T1,T2,T3,T4,T5,T6] | [0.5,1,1.5,2,2.5,3] |
Time step size, Δτ | 1/360 |
Computational domain, Ω | (0,300) |
Current stock price, S0 | 100 |
Strike prices, [K1,K2,K3,K4,K5,K6] | [95,95,90,90,85,85] |
Knock-in barrier, D | 50 |
Coupon rates, [c1,c2,c3,c4,c5,c6] | [0.02,0.04,0.06,0.08,0.1,0.12] |
Dummy rate, d | 0.12 |
Face value, F | 100 |
Risk-free interest, r | 0.03 |
Volatility, σ | 0.2 |
The value of ELS exhibits sharp changes near the knock-in barrier and strike price due to rapid variations in the payoff. Therefore, to accurately capture these changes and minimize numerical errors, we calculate its value using a finer grid in these regions, as shown in Figure 6. This approach improves numerical accuracy by mitigating approximation errors, which can become significant when using a coarse grid in areas with large changes (Lyu et al., 2021). In particular, the rapid variations in the payoff near the knock-in barrier and strike price require a finer grid for more precise approximations of derivatives, ensuring a more accurate evaluation of the option price.
The variation in the value of the option due to small fluctuations in the underlying asset price is indicated by the Δ, helping issuers decide on the initial coupon rate and the step-down mechanisms based on their risk tolerance and market outlook. Moreover, Δ plays an important role in measuring the sensitivity of ELS to price, thereby helping to construct a portfolio that reduces risk (Kim et al., 2021). It can be calculated using finite difference approximations as follows.
Δni=(∂u∂S)ni=−hiuni−1hi−1(hi−1+hi)+(hi−hi−1)unihi−1hi+hi−1uni+1hi(hi−1+hi). |
Figure 7 shows the Δ surface as time and the strike price change. The Δ value increases and declines rapidly near the knock-in barrier because of Δ properties. Normally, rises as the price of the underlying asset falls and declines as the price of the underlying asset increases. The majority of ELS products currently available are step-down structures, leading to positive Δ values. Hence, it is essential to consider this and implement hedging strategies (Kanamura (2018); Larguinho et al. (2022)), such as selling the relevant underlying assets.
The variation of Δ is measured by Γ, which reflects the convexity of the payoff of the product. It represents the rate at which Δ changes in response to variations. Essentially, Γ is a partial differentiation of underlying assets for Δ, which is a second-order partial derivative of underlying asset with respect to options. The computation of Γ entails numerical discretization at Si, as outlined below.
Γni=(∂2u∂S2)ni=2(uni−1hi−1(hi−1+hi)−unihi−1hi+uni+1hi(hi−1+hi)). |
Figure 8 shows the profile of the Γ value. As the maturity date approaches, and the underlying asset is near the knock-in barrier, Γ takes on positive values. In contrast, at each early redemption date or when the underlying asset price is around S0, Γ tends to have negative values. When Γ takes on significantly negative or positive values, Δ becomes more sensitive to changes in the price of the underlying asset. For this reason, it can be inferred from Figure 7 that the Δ value responds sensitively when the underlying asset is near the knock-in barrier.
The fluctuation in the value of a portfolio as time t increases is measured by Θ, which is algebraically defined as the partial derivative of the option value with respect to time. We use the central difference method to calculate Θ as follows:
Θni=(∂u∂t)ni=un+1i−un−1i2Δt. |
Figure 9 displays the Θ profile and the Θ surface over time, respectively. The Θ value is generally negative because the value of an option tends to decrease over time as it approaches maturity.
The rate at which option prices fluctuate due to the volatility of the underlying assets is denoted by ν. Some option values may be sensitive to variations in volatility, and ν plays an important role in monitoring markets that are particularly vulnerable to volatility. We discretize ν using the central difference method for numerical computation as follows:
νni=(∂u∂σ)ni=uni(σ+Δσ)−uni(σ−Δσ)2Δσ, |
where Δσ=0.01 is the size of increment of volatility we used.
Figure 10 shows the profile of the ν value. The ν value is large and negative when the underlying asset is near the knock-in barrier and as the redemption date approaches maturity. At each redemption point, the ν value first increases and then decreases. A large absolute ν indicates high sensitivity to changes in volatility. When hedging ν, it is crucial to take into account both the time intervals between rebalancing events and the fluctuations in volatility.
The sensitivity of the option price to variations in the risk-free interest rate is quantified by ρ. A high ρ means that the change in interest has a high impact on the option price. Applying the second order central difference method, ρ can be calculated on the computational grid.
ρni=(∂u∂r)ni=uni(r+Δr)−uni(r−Δr)2Δr. | (15) |
where Δr=0.001 is the size of increment for ρ we used.
Figure 11 shows the ρ surface. Regardless of the redemption point, ρ takes a positive value within the knock-in barrier range. As the price of the underlying asset increases, ρ rapidly decreases and becomes a large negative value. As the underlying asset price approaches S0, ρ starts to increase again and eventually tends toward zero. Similar to ν, ρ increases sharply as the redemption point approaches and then decreases rapidly.
Finally, we calculate the values of different structural ELS products. Unlike the payoffs in Equations (1) and (2) discussed above, we consider an ELS product with the following payoff. If the knock-in barrier is touched,
v(S,Ti)={F(1+c6),if S≥K6,0.5F(S+K6)/S0,otherwise, | (16) |
if the knock-in barrier is not touched,
u(S,Ti)={F(1+c6),if S≥K6,F(1+d),if D<S<K6,0.5F(S+K6)/S0,otherwise. | (17) |
The payoff described above differs from the previously discussed payoffs in Equations (1) and (2) in that it provides a constant profit even if the value of the underlying asset S becomes 0 at maturity. See Figure 12.
When S=0 for products whose future profit is 0, it is calculated using the zero Dirichlet boundary at i=1. However, when S=0 for products whose future profit is not 0, the discounted price of the future profit must be considered at i=1. In this case, either the Dirichlet boundary for the discounted value at i=1 or the linear boundary can be applied. The values and Greeks of ELS calculated by applying the linear boundary are consistent with previous results. Figure 13 shows the values and Greeks of ELS products calculated by applying linear boundary to payoff Equations (16) and (17). The parameters used in this computational test are the same as in Table 1.
In this study, we developed a computational method for the accurate computation of the option prices and Greeks for the ELS near early redemption dates. We used the BS equation and FDM to compute the ELS price and its Greeks. The proposed method solved the BS equation at early redemption dates using the Dirichlet boundary condition at strike prices instead of overwriting the option values on greater than or equal to strike prices by the given option prices. As a result, we observed that our computing algorithm is precise and useful for ELS pricing and its Greeks close to each early redemption date and the maturity time. Greeks have been measured to hedge the risk of derivative products, and Greeks hedging can be applied to manage the risk of ELS products. The computation of Greeks of ELS not only provides better hedging strategies for the finance industry but also facilitates risk- reduction decisions. For future research, a higher-order numerical scheme, such as the Crank–Nicolson method (Roul and Goura, 2020), can be applied to improve the accuracy of time discretization. Moreover, machine learning techniques should be explored to enhance both the efficiency and accuracy of numerical methods. For instance, machine learning approaches should be considered for optimizing hedging strategies based on dynamically changing market conditions (Wang and Yan, 2023). Additionally, deep learning-based surrogate models can be developed to approximate option price surfaces efficiently, which reduces computational costs while maintaining precision (Anderson and Ulrych, 2023).
Yunjae Nam: Conceptualization, Data curation, Resources, Formal analysis, Software, Visualization, Writing – original draft, and Writing – review & editing. Changwoo Yoo: Formal analysis, Resources, Data curation, Validation, and Writing – original draft. Hyundong Kim: Validation, Formal analysis, Resources, Investigation, Visualization, and Writing – review & editing. Jaewon Hong: Validation, Investigation, Data curation, and Writing – original draft. Minjoon Bang: Validation, Investigation, Data curation, and Writing – review & editing. Junseok Kim: Conceptualization, Supervision, Project administration, Funding acquisition, Formal analysis, Resources, Methodology, Writing – original draft, and Writing – review & editing.
This article was not created using Artificial Intelligence (AI) tools.
This work was supported by the National Research Foundation (NRF), Korea, under project BK21 FOUR. The authors are grateful to the reviewers for their insightful and constructive feedback.
The authors declare no conflict of interest in this paper.
[1] |
Anderson D, Ulrych U (2023) Accelerated American option pricing with deep neural networks. Quant Financ Econ 7: 207–228. https://doi.org/10.3934/QFE.2023011 doi: 10.3934/QFE.2023011
![]() |
[2] |
Black F, Scholes M (1973) The pricing of options and corporate liabilities. J Polit Econ 81: 637–654. https://doi.org/10.1086/260062 doi: 10.1086/260062
![]() |
[3] | Cui Y, Li L, Zhang G (2024) Pricing and hedging autocallable products by Markov chain approximation. Rev Deriv Res 27: 1–45. |
[4] |
Hwang Y, Kim I, Kwak S, et al. (2023) Unconditionally stable monte carlo simulation for solving the multi-dimensional Allen–Cahn equation. Electron Res Arch 31. https://doi.org/10.3934/era.2023261 doi: 10.3934/era.2023261
![]() |
[5] |
Jo J, Kim Y (2013) Comparison of numerical schemes on multi-dimensional Black–Scholes equations. Bull Korean Math Soc 50: 2035–2051. https://doi.org/10.4134/BKMS.2013.50.6.2035 doi: 10.4134/BKMS.2013.50.6.2035
![]() |
[6] |
Kanamura T (2018) Diversification effect of commodity futures on financial markets. Quant Financ Econ 2: 821–836. https://doi.org/10.3934/QFE.2018.4.821 doi: 10.3934/QFE.2018.4.821
![]() |
[7] |
Kim ST, Kim HG, Kim JH (2021) ELS pricing and hedging in a fractional Brownian motion environment. Chaos Solitons Fractals 142: 110453. https://doi.org/10.1016/j.chaos.2020.110453 doi: 10.1016/j.chaos.2020.110453
![]() |
[8] |
Kwak S, Kang S, Ham S, et al. (2023) An unconditionally stable difference scheme for the two‐dimensional modified Fisher–Kolmogorov–Petrovsky–Piscounov equation. J Math 2023: 5527728. https://doi.org/10.1155/2023/5527728 doi: 10.1155/2023/5527728
![]() |
[9] |
Larguinho M, Dias JC, Braumann CA (2022) Pricing and hedging bond options and sinking-fund bonds under the CIR model. Quant Financ Econ 6: 1–34. https://doi.org/10.3934/QFE.2022001 doi: 10.3934/QFE.2022001
![]() |
[10] |
Lee H, Ha H, Kong B, et al. (2024) Valuing three-asset barrier options and autocallable products via exit probabilities of Brownian bridge. N Am Econ Financ 73: 102174. https://doi.org/10.1016/j.najef.2024.102174 doi: 10.1016/j.najef.2024.102174
![]() |
[11] |
Lee C, Kwak S, Hwang Y, et al. (2023) Accurate and efficient finite difference method for the Black–Scholes model with no far-field boundary conditions. Comput Econ 61: 1207–1224. https://doi.org/10.1007/s10614-022-10242-w doi: 10.1007/s10614-022-10242-w
![]() |
[12] |
Liu T, Li T, Ullah MZ (2024) On five-point equidistant stencils based on Gaussian function with application in numerical multi-dimensional option pricing. Comput Math Appl 176: 35–45. https://doi.org/10.1016/j.camwa.2024.09.003 doi: 10.1016/j.camwa.2024.09.003
![]() |
[13] |
Liu T, Soleymani F, Ullah MZ (2024) Solving multi-dimensional European option pricing problems by integrals of the inverse quadratic radial basis function on non-uniform meshes. Chaos Solitons Fractals 185: 115156. https://doi.org/10.1016/j.chaos.2024.115156 doi: 10.1016/j.chaos.2024.115156
![]() |
[14] |
Lyu J, Park E, Kim S, et al. (2021) Optimal non-uniform finite difference grids for the Black–Scholes equations. Math Comput Simul 182: 690–704. https://doi.org/10.1016/j.matcom.2020.12.002 doi: 10.1016/j.matcom.2020.12.002
![]() |
[15] |
Roul P, Goura VP (2020) A sixth order numerical method and its convergence for generalized Black–Scholes PDE. J Comput Appl Math 377: 112881. https://doi.org/10.1016/j.cam.2020.112881 doi: 10.1016/j.cam.2020.112881
![]() |
[16] |
Tao L, Lai Y, Ji Y, et al. (2023) Asian option pricing under sub-fractional Vasicek model. Quant Finan Econ 7: 403–419. https://doi.org/10.3934/QFE.2023020 doi: 10.3934/QFE.2023020
![]() |
[17] | Thomas Thomas L (1949) Elliptic Problems in Linear Differential Equations Over a Network: Watson Scientific Computing Laboratory. Columbia University: New York, NY, USA |
[18] |
Wang Y, Yan K (2023) Machine learning-based quantitative trading strategies across different time intervals in the American market. Quant Financ Econ 7: 569–594. https://doi.org/10.3934/QFE.2023028 doi: 10.3934/QFE.2023028
![]() |
[19] |
Wu X, Wen S, Shao W, et al. (2023) Numerical Investigation of Fractional Step-Down ELS Option. Fractal Fract 7: 126. https://doi.org/10.3390/fractalfract7020126 doi: 10.3390/fractalfract7020126
![]() |
Parameter | Value |
Maturity time, T | 3 |
Repayment times, [T1,T2,T3,T4,T5,T6] | [0.5,1,1.5,2,2.5,3] |
Time step size, Δτ | 1/360 |
Computational domain, Ω | (0,300) |
Current stock price, S0 | 100 |
Strike prices, [K1,K2,K3,K4,K5,K6] | [95,95,90,90,85,85] |
Knock-in barrier, D | 50 |
Coupon rates, [c1,c2,c3,c4,c5,c6] | [0.02,0.04,0.06,0.08,0.1,0.12] |
Dummy rate, d | 0.12 |
Face value, F | 100 |
Risk-free interest, r | 0.03 |
Volatility, σ | 0.2 |
Parameter | Value |
Maturity time, T | 3 |
Repayment times, [T1,T2,T3,T4,T5,T6] | [0.5,1,1.5,2,2.5,3] |
Time step size, Δτ | 1/360 |
Computational domain, Ω | (0,300) |
Current stock price, S0 | 100 |
Strike prices, [K1,K2,K3,K4,K5,K6] | [95,95,90,90,85,85] |
Knock-in barrier, D | 50 |
Coupon rates, [c1,c2,c3,c4,c5,c6] | [0.02,0.04,0.06,0.08,0.1,0.12] |
Dummy rate, d | 0.12 |
Face value, F | 100 |
Risk-free interest, r | 0.03 |
Volatility, σ | 0.2 |