After the discovery of the fractal structures of financial markets, enormous effort has been dedicated to finding accurate and stable numerical schemes to solve fractional Black-Scholes partial differential equations. This work, therefore, proposes a numerical scheme for pricing double-barrier options, written on an underlying stock whose dynamics are governed by a non-standard fractal stochastic process. The resultant model is time-fractional and is herein referred to as a time-fractional Black-Scholes model. The presence of the time-fractional derivative helps to capture the time-decaying effects of the underlying stock while capturing the globalized change in underlying prices and barriers. In this paper, we present the construction of the proposed scheme, analyse it in terms of its stability and convergence, and present two numerical examples of pricing double knock-in barrier-option problems. The results suggest that the proposed scheme is unconditionally stable and convergent with order $\mathcal{O}(h^2+k^2)$.

Since the early 1970s, the study of Black–Scholes (BS) partial differential equations (PDEs) under the Efficient Market Hypothesis (EMH) has been a subject of active research in financial engineering. It has now become obvious, even to casual observers, that the classical BS models derived under the EMH framework fail to account for a number of realistic price evolutions in real-time market data. An alternative approach to the EMH framework is the Fractal Market Hypothesis (FMH), which proposes better and clearer explanations of market behaviours during unfavourable market conditions. The FMH involves non-local derivatives and integral operators, as well as fractional stochastic processes, which provide better tools for explaining the dynamics of evolving market anomalies, something that classical BS models may fail to explain. In this work, using the FMH, we derive a time-fractional Black–Scholes partial differential equation (tfBS-PDE) and then transform it into a heat equation, which allows for ease of implementing a high-order numerical scheme for solving it. Furthermore, the stability and convergence properties of the numerical scheme are discussed, and overall techniques are applied to pricing European put option problems.

We validate the new Basel liquidity standards as encapsulated by the net stable funding ratio in a quantitative manner. In this regard, we consider the dynamics of inverse net stable funding ratio as a measure to quantify the bank’s prospects for a stable funding over a period of a year. In essence, this justifies how Basel III liquidity standards can be effectively implemented in mitigating liquidity problems. We also discuss various classes of available stable funding and required stable funding. Furthermore, we discuss an optimal control problem for a continuous-time inverse net stable funding ratio. In particular, we make optimal choices for the inverse net stable funding targets in order to formulate its cost. This is normally done by obtaining analytic solution of the value function. Finally, we provide a numerical example for the dynamics of the inverse net stable funding ratio to identify trends in which banks behavior convey forward looking information on long-term market liquidity developments.