ε-Strong Simulation of Fractional Brownian Motion and Related Stochastic Differential Equations