The Continuous Galerkin Method for Fractional Order Viscoelasticity