{"title":"A new ADI framework for oil option pricing with cubic B-spline and Gauss-Hermite integration","authors":"Fares Alazemi, Abdulaziz Alsenafi","doi":"10.1016/j.cam.2026.117439","DOIUrl":null,"url":null,"abstract":"<div><div>This paper investigates the pricing of European options within the commodity market, specifically focusing on the dynamics of oil spot prices. We employ a long memory stochastic volatility model to effectively capture and predict the behavior of oil prices. Initially, we utilize statistical tests such as the Augmented Dickey-Fuller (ADF) test and other relevant techniques to extract essential characteristics from oil spot price data, thereby establishing a foundation for our analysis. To model these features, we propose a two-factor financial model grounded in a mixed fractional Brownian motion process. This approach is crucial, as it recognizes that the increments of the process exhibit dependence and retain long memory characteristics for values of <em>H</em> within the range <span><math><mrow><mo>(</mo><mfrac><mn>3</mn><mn>4</mn></mfrac><mo>,</mo><mn>1</mn><mo>)</mo></mrow></math></span>. However, this complexity poses challenges in deriving an explicit formula for forward contract pricing in commodity markets due to the inherent non-normal behavior of the underlying asset. To address this issue, we implement an advanced Monte Carlo method to approximate the prices of the contracts. This approach enables us to simulate a large number of potential price paths for oil, thereby facilitating more accurate pricing under the observed stochastic volatility conditions. Additionally, we incorporate stochastic calculus techniques alongside a delta hedging strategy to derive the partial differential equation (PDE) governing the option price. To solve this equation, we employ the Alternating Direction Implicit (ADI) method, known for its efficiency and stability in managing high-dimensional problems. Unlike explicit methods, which can suffer from stability issues and require very small time steps, the ADI method allows for larger time steps while maintaining stability, making it computationally efficient. This is particularly important in option pricing, where the financial markets can exhibit rapid changes. By leveraging these methodologies, the project aims to enhance the understanding of option pricing dynamics in the commodity market, providing valuable insights for traders and financial analysts involved in oil market operations.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"484 ","pages":"Article 117439"},"PeriodicalIF":2.4000,"publicationDate":"2026-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational and Applied Mathematics","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0377042726001044","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"2026/2/11 0:00:00","PubModel":"Epub","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 0
Abstract
This paper investigates the pricing of European options within the commodity market, specifically focusing on the dynamics of oil spot prices. We employ a long memory stochastic volatility model to effectively capture and predict the behavior of oil prices. Initially, we utilize statistical tests such as the Augmented Dickey-Fuller (ADF) test and other relevant techniques to extract essential characteristics from oil spot price data, thereby establishing a foundation for our analysis. To model these features, we propose a two-factor financial model grounded in a mixed fractional Brownian motion process. This approach is crucial, as it recognizes that the increments of the process exhibit dependence and retain long memory characteristics for values of H within the range . However, this complexity poses challenges in deriving an explicit formula for forward contract pricing in commodity markets due to the inherent non-normal behavior of the underlying asset. To address this issue, we implement an advanced Monte Carlo method to approximate the prices of the contracts. This approach enables us to simulate a large number of potential price paths for oil, thereby facilitating more accurate pricing under the observed stochastic volatility conditions. Additionally, we incorporate stochastic calculus techniques alongside a delta hedging strategy to derive the partial differential equation (PDE) governing the option price. To solve this equation, we employ the Alternating Direction Implicit (ADI) method, known for its efficiency and stability in managing high-dimensional problems. Unlike explicit methods, which can suffer from stability issues and require very small time steps, the ADI method allows for larger time steps while maintaining stability, making it computationally efficient. This is particularly important in option pricing, where the financial markets can exhibit rapid changes. By leveraging these methodologies, the project aims to enhance the understanding of option pricing dynamics in the commodity market, providing valuable insights for traders and financial analysts involved in oil market operations.
期刊介绍:
The Journal of Computational and Applied Mathematics publishes original papers of high scientific value in all areas of computational and applied mathematics. The main interest of the Journal is in papers that describe and analyze new computational techniques for solving scientific or engineering problems. Also the improved analysis, including the effectiveness and applicability, of existing methods and algorithms is of importance. The computational efficiency (e.g. the convergence, stability, accuracy, ...) should be proved and illustrated by nontrivial numerical examples. Papers describing only variants of existing methods, without adding significant new computational properties are not of interest.
The audience consists of: applied mathematicians, numerical analysts, computational scientists and engineers.