Summary: We introduce a finite element approximation of a variational formulation of Bean’s model for the physical configuration of an infinitely long cylindrical superconductor subject to a transverse magnetic field. We prove an error between the exact solution and the approximate solution for the current density and the magnetic field in appropriate norms of order \(h^{1/2}+\Delta t\). Numerical simulations for a variety of applied magnetic fields are also presented.