We introduce a stabilized conforming mixed finite element method for a macroscopic model in micromagnetics. We show well-posedness of the discrete problem for higher order elements in two and three dimensions, develop a full a priori analysis for lowest order elements, and discuss the extension of the method to higher order elements. We introduce a residual-based a posteriori error estimator and present an adaptive strategy. Numerical examples illustrate the performance of the method.