“…As the random driving source of the risk asset price in the BS model is Brownian motion, it cannot capture many of the characteristic features of the risk asset price, including long-range dependence, heavy tails, and periods of constant values, etc. In order to overcome these shortcomings, scholars adopt other stochastic processes as the random source of the option pricing model, such as fractional Brownian motion [2-7], mixed fractional Brownian motion [8][9][10][11][12], time-changed Brownian motion [13][14][15][16], skew Brownian motion [17][18][19][20][21], etc.…”