Hybrid Taylor and block-pulse functions operational matrix algorithm and its application to obtain the approximate solution of stochastic evolution equation driven by fractional Brownian motion.