A Numerical Scheme for Dynamic Systems Containing Fractional Derivatives