  We derive the equations of nonlinear electroelastostatics using three
different variational formulations involving the deformation function and an
independent field variable representing the electric character - considering
either one of the electric field $\mathbb{E}$, electric displacement
$\mathbb{D}$, or electric polarization $\mathbb{P}$. The first variation of the
energy functional results in the set of Euler-Lagrange partial differential
equations which are the equilibrium equations, boundary conditions, { and
certain constitutive equations} for the electroelastic system. The partial
differential equations for obtaining the bifurcation point have been also found
using the second variation based bilinear functional. We show that the
well-known Maxwell stress in vacuum is a natural outcome of the derivation of
equations from the variational principles and does not depend on the
formulation used. As a result of careful analysis it is found that there are
certain terms in the bifurcation equation which appear difficult to obtain by
an ordinary perturbation based analysis of the Euler-Lagrange equation. From a
practical viewpoint, the formulations based on $\mathbb{E}$ and $\mathbb{D}$
result in simpler equations and are anticipated to be more suitable for
analysing problems of stability as well as post-buckling behaviour.
