The problem of calculations of functionals of the fractional Brownian motion (fBm) often arises in different applications of modern finance and statistics of random processes, but analytical methods for solving this problem are yet to be developed. At the beginning of the talk we shall describe a fast method for simulation of trajectories of fBm which is based on use of parallel computations. Then we illustrate the performance of the method for two long-standing unsolved problems: calculation of moments of integral functionals of fBm and a distribution of maximum fBm.