A very efficient approach for pricing barrier options on an underlying described by the mixed fractional Brownian motion

We deal with the problem of pricing barrier options on an underlying described by the mixed fractional Brownian model. To this aim, we consider the initial-boundary value partial differential problem that yields the option price and we derive an integral representation of it in which the integrand functions must be obtained solving Volterra equations of the first kind. In addition, we develop an ad-hoc numerical procedure to solve the integral equations obtained. Numerical simulations reveal that the proposed method is extremely accurate and fast, and performs significantly better than the finite difference method.