Stochastic differential equations (SDEs) are very useful in modeling many problems in biology, economic data, turbulence, and medicine. Fractional Brownian motion (fBm) and variable‐order fractional Brownian motion (vofBm) are suitable alternatives to standard Brownian motion (sBm) for describing and modeling many phenomena, since the increments of these processes are dependent of the past and for
these increments have the property of long‐term dependence. Classical mathematical techniques such as Ito's calculus do not work for stochastic computations on fBm and vofBm due to they are not semi‐Martingale for
. Therefore, solving these equations is much more difficult than solving SDEs with sBm. On the other hand, these equations do not have an analytical solution, so we have to use numerical methods to find their solution. In this paper, a computational approach based on hybrid of block‐pulse and parabolic functions (HBPFs) has been introduced for simulating vofBm and solving a modern class of SDEs. The mechanism of this approach is based on stochastic and fractional integration operational matrices, which transform the intended problem to a nonlinear system of algebraic equations. Thus, the complexity of solving the mentioned problem is reduced significantly. Also, convergence analysis of the expressed method has been theoretically examined. Finally, the accuracy and efficiency of the proposed algorithm have been experimentally investigated through some test problems and comparison of obtained results with results of previous papers. High accurate numerical results are obtained by using a small number of basic functions. Therefore, this method deals with smaller matrices and vectors, which is one of the most important advantage of our suggested method. Also, presenting an applicable procedure to construct vofBm is another innovation of this work. To gain this aim, at first, discretized sBm is generated via fundamental features of this process, and afterward, block‐pulse functions (BPFs) and HBPFs are utilized for simulating discretized vofBm. Finally, spline interpolation method has been employed to provide a continuous path of vofBm.