$${\varrho=\frac{1}{2} \varepsilon_0 (\varepsilon_x E_x^2+ \varepsilon_y E_y^2 +\varepsilon_z E_z^2 )}$$