Research article

Accurate computation of Greeks for equity-linked security (ELS) near early redemption dates

  • Received: 14 January 2025 Revised: 17 March 2025 Accepted: 21 April 2025 Published: 24 April 2025
  • JEL Codes: C63, G13, G12

  • 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

    Related Papers:

    [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 K1K2K3K4K5K6, and the coupon rates set as c1c2c3c4c5c6 for each early redemption date T1T2T3T4T5T6.

    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 SK6,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 SK6,F(1+d),if D<S<K6,FS/S0,otherwise. (2)
    Figure 1.  Function of payoff at early redemption dates.

    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)222v(S,t)S2+rSv(S,t)Srv(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 τ=Tt. Then, Equation (3) becomes

    v(S,τ)τ=(σS)222v(S,τ)S2+rSv(S,τ)Srv(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+1ivniΔτ=(σSi)22(2vn+1i1hi1(hi1+hi)2vn+1ihi1hi+2vn+1i+1hi(hi1+hi))+rSi(hivn+1i1hi1(hi1+hi)+(hihi1)vn+1ihi1hi+hi1vn+1i+1hi(hi1+hi))rvn+1i, (5)

    where vniv(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=2SNxSNx1 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.

    Figure 2.  Schematic of discrete domain.

    Equation (5) can be rewritten as follows:

    αivn+1i1+βivn+1i+γivn+1i+1=fi,i=2,,Nx, (6)

    where

    αi=(σSi)2+rSihihi1(hi1+hi),βi=1Δτ+(σSi)2rSi(hihi1)hi1hi+r, (7)
    γi=(σSi)2rSihi1hi(hi1+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+1Nxvn+1Nx1 and get

    (αNxγNx)vn+1Nx1+(βNx+2γNx)vn+1Nx=fNx. (9)

    Then, the solution vector vn+12:Nx=[vn+12vn+13vn+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γ200000α3β3γ300000α4β4γ40000000αNx1βNx1γNx100000αNxγNxβNx+2γNx)

    and

    f2:Nx=(f2f3f4fNx1fNx).

    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:Kj1=[vn+12vn+13vn+1Kj1]T can be determined by solving the tridiagonal system

    Avn+12:Kj1=f2:Kj1, (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:Kj1=(β2γ200000α3β3γ300000α4β4γ40000000αKj2βKj2γKj200000αKj1βKj1)

    and

    f2:Kj1=(f2f3f4fKj1fKjγKj1vn+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)
    αKj1un+1Kj2+βKj1un+1Kj1=gKj1γKj1un+1Kj. (13)
    Figure 3.  Option values v at before redemption date and redemption date.
    Figure 4.  Option values u without knock and v with knock.

    Then, the solution vector un+1D+1:Kj1=[un+1D+1un+1D+2un+1Kj1]T can be obtained by resolving the tridiagonal system

    Bun+1D+1:Kj1=gD+1:Kj1, (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:Kj1=(βD+1γD+100000αD+2βD+2γD+200000αD+3βD+3γD+30000000αKj2βKj2γKj200000αKj1βKj1) and fD+1:Kj1=(fD+1αD+1un+1DfD+2fD+2fKj2fKj1γKj1un+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.

    Figure 5.  Option values u at before redemption date and redemption date.

    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.

    Table 1.  Parameters used in Section 3.
    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

     | Show Table
    DownLoad: CSV

    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.

    Figure 6.  Non-uniform mesh grid for computational experiments.

    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=(uS)ni=hiuni1hi1(hi1+hi)+(hihi1)unihi1hi+hi1uni+1hi(hi1+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.

    Figure 7.  Δ surface of ELS at each early redemption dates and until maturity under specific conditions and its corresponding conditioned price.

    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=(2uS2)ni=2(uni1hi1(hi1+hi)unihi1hi+uni+1hi(hi1+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.

    Figure 8.  Γ surface of ELS at each early redemption dates and until maturity under specific conditions and its corresponding conditioned price.

    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=(ut)ni=un+1iun1i2Δ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.

    Figure 9.  Θ surface of ELS at each early redemption dates and until maturity under specific conditions and its corresponding conditioned price.

    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.

    Figure 10.  ν surface of ELS at each early redemption dates and until maturity under specific conditions and its corresponding conditioned price.

    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=(ur)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.

    Figure 11.  ρ surface of ELS at each early redemption dates and until maturity under specific conditions and its corresponding conditioned price.

    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 SK6,0.5F(S+K6)/S0,otherwise, (16)

    if the knock-in barrier is not touched,

    u(S,Ti)={F(1+c6),if SK6,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.

    Figure 12.  Payoff function for the early redemption date across different scenarios.

    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.

    Figure 13.  Numerical price and Greeks for payoffs in Equations (1) and (2).

    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
  • Reader Comments
  • © 2025 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(405) PDF downloads(69) Cited by(0)

Figures and Tables

Figures(13)  /  Tables(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog