Numerics for the fractional Langevin equation driven by the fractional Brownian motion