A semi-implicit method for solving the 3-dimensional magnetohydrodynamic equations on long time scales is presented. Standard explicit methods must use time steps which are constrained by a Courant-Friedrichs-Lewy condition due to the fast compressional and shear alfven motion. This semi-implicit method eliminates both of these restrictions so that very large time steps are permitted. The method is simple to implement and the computation time for one time step is essentially the same as for explicit methods. Numerical test results in slab and cylindrical geometry are presented.