Anisotropic polynomial constitutive equations for electroelastic crystals