We propose a numerical procedure for nonlinear electromechanical problems in smart ceramics. Two essential features of the present procedure are the use of a Landau's phenomenological potential energy as nonlinear constitutive law and the finite element discretization of an alternative governing equation with vector potential. For solving the resulting semi-discrete nonlinear system, the Newmark-β method is employed. Numerical examples in two dimensions are shown to validate the present procedure.