What Murphrid describes is usually known as "Gauß's Theorem" and refers to a mathematical theorem in vector calculus in threedimensional euclidean space. It says
\int_{V} \mathrm{d}^3 \vec{x} \vec{\nabla} \cdot \vec{A}(\vec{x})=\int_{\partial V} \mathrm{d}^2 \vec{F} \cdot \vec{A}(\vec{x}).
Here \vec{A} is a sufficiently well-behaved vector field, V a volume with sufficiently "nice" shape and boundary in threedimensional Euclidean space, and \partial V its boundary. The orientation of the surface-normal vectors \mathrm{d}^2 \vec{F} is such that they point out of the volume V.
Gauß's Law is one fundamental equation in electromagnetism and thus part of Maxwell's equations. It says (in Heaviside-Lorentz Units)
\vec{\nabla} \cdot \vec{E}=\rho.
Here \vec{E} are the electric components of the electromagnetic field (mostly just called "the electric field") and \rho the density of electric charges.
Using Gauß's theorem, you get the integral form of the same law,
\int_{\partial_V} \mathrm{d}^2 \vec{F} \cdot \vec{E} = \int_V \mathrm{d}^3 \vec{x} \rho(\vec{x})=Q_{V}.
Again the orientation of the surface vectors is as explained above. It tells you that the electric flux, i.e., the surface integral over the electric field on the left-hand side of the equation always equals the total charge contained in the volume, bounded by the surface.
Sometimes, if a problem is sufficiently symmetric, you can use the integral form to determine the electric field. The most simple case is the field of a point charge Q, which we put for simplicity in the origin of the coordinate system. Due to spherical symmetry, we expect that the electric field is radial and its magnitude only depends on the distance from the charge, i.e., we make the ansatz
\vec{E}=E(r) \vec{e}_r,
where \vec{e}_r=\vec{x}/|\vec{x}| is the radial unit vector.
Now take a sphere S_R of radius R around the origin and use spherical coordinates for the surface integral. The surface normal vectors in standard spherical coordinates (r,\vartheta,\varphi) are
\mathrm{d}^2 \vec{F}=\mathrm{d} \vartheta \; \mathrm{d}\varphi \; \sin \vartheta \vec{e}_r.
From this you get for Gauß's Law
4 \pi R^2 E_r(R)=Q
and thus immediately the Coulomb-field solution
E_r(R)=\frac{Q}{4 \pi R^2}.
It is pretty easy to show that indeed
\vec{\nabla} \cdot \vec{E}=0 \quad \text{for} \quad \vec{x} \neq 0.
That the factor is correct, we have just proven with the integral form of Gauß's Law. At this step you necessarily need the integral form, because a point charge is a singularity, because the corresponding charge density is given by a Dirac \delta distribution,
\rho(\vec{x})=Q \delta^{(3)}(\vec{x}).