A finite-element A-V formulation is developed to solve EMQS (electromagnetic quasistatic), aka Darwin, problems in the time and frequency domains by making use of a two-domain continuity gauging strategy. The formulation utilizes the magnetic vector and the electric scalar potentials, and an additional scalar in the conducting regions to make the discretized Galerkin system symmetric. The resulting FE (finite element) matrix is LF (low-frequency) stable in the frequency domain. As a result, the time-step size is not limited by the characteristics of the system matrix in the backward Euler time-stepping scheme. The system matrixes are also regular, which allows choosing a direct solver to speed up the temporal solution significantly. The formulation is verified in both the frequency and the time domains by solving benchmark problems, and the results are compared to full-wave and eddy-current solutions.