Solution
= Solution
The <divergence theorem> is $\iint_{\partial V}F\cdot n\,dS=\iiint_V\nabla\cdot F\,dV$. Here $\nabla\cdot F=3z+2yz$. The $y$ term integrates to zero by symmetry, while the ellipse has area $\pi/\sqrt{ab}$. Hence the flux is
$$
\boxed{\frac\pi{\sqrt{ab}}\int_0^3 3z\,dz=\frac{27\pi}{2\sqrt{ab}}.}
$$