div\vec{E}=\frac{1}{\epsilon}\rho