Generating sample paths and their convergence of the geometric fractional brownian motion

We derive discrete time model of the geometric fractional Brownian motion. It provides numerical pricing scheme of financial derivatives when the market is driven by geometric fractional Brownian motion. With the convergence analysis, we guarantee the convergence of Monte Carlo simulations. The strong convergence rate of our scheme has order H which is Hurst parameter. To obtain our model we need to convert Wick product term of stochastic differential equation into Wick free discrete equation through Malliavin calculus but ours does not include Malliavin derivative term. Finally, we include several numerical experiments for the option pricing.