A fractional Black-Scholes model with stochastic volatility and European option pricing

•European option pricing is studied under a FMLS model with stochastic volatility.•We have successfully solved a fractional partial differential equation system.•The derived pricing formula is truly explicit, involving no Fourier inversion. In this paper, we introduce the stochastic volatility into the FMLS (finite moment log-stable) model to capture the effect of both jumps and stochastic volatility. However, this additional stochastic source adds another degree of complexity in seeking for analytical formula when pricing European options, as the involved FPDE (fractional partial differential equation) system governing option prices is now of three dimensions. Albeit difficult, we have still managed to present an analytical solution expressed in terms of Fourier cosine series, after a two-step solution procedure is developed for the target FPDE system. This solution is different from the most literature as it is truly explicit, involving no Fourier inversion. It is also shown through the numerical experiments that it converges very rapidly and has potential to be applied in practice.