Numerical Methods for Two-Fluid Plasma Equations