An accurate and stable numerical method for option hedge parameters

Junhyun Cho, Yejin Kim, Sungchul Lee

Research output: Contribution to journalArticlepeer-review

2 Citations (Scopus)


We propose an unconditionally stable numerical algorithm, which uses the Feynman-Kac formula of the Black-Scholes equation to obtain accurate option prices and hedge parameters. We discretize the asset and time using uniform grid points. We approximate the option values by piecewise quadratic polynomials for each time step and integrate them analytically over the log-normal distribution. The piecewise quadratic approximation gives the third-order convergence in the asset direction, and the analytic integration reduces truncation error in the time direction. The estimation errors are propagated backward in time following the convection and diffusion characteristics of the Black-Scholes equation, which assures the unconditional stability of our method. The vectorized code implementation reduces the time complexity. The convergence test shows that our approach outperforms the Crank-Nicolson scheme of the finite difference method in both time and asset directions, and the stability test verifies that our method is stable as the Crank-Nicolson. Furthermore, we show that our algorithm reduces the price errors and hedge parameter errors by more than 50% from the benchmark.

Original languageEnglish
Article number127276
JournalApplied Mathematics and Computation
Publication statusPublished - 2022 Oct 1

Bibliographical note

Publisher Copyright:
© 2022 Elsevier Inc.

All Science Journal Classification (ASJC) codes

  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'An accurate and stable numerical method for option hedge parameters'. Together they form a unique fingerprint.

Cite this