This chapter discusses 2 X 2 symmetric variational formulations and associated finite element methods for the EMI equations. We demonstrate that the presented methods converge at expected rates, and compare the approaches in terms of approximation of the transmembrane potential. Overall, the choice of which formulation to employ for solving EMI models becomes largely a matter of desired accuracy and available computational resources.