 Solution

The divergence theorem states that a continuously differentiable vector field on a bounded region with piecewise smooth boundary satisfieswith the unit normal pointing out of $V$. Apply it separately to $\mathbf{u}=f{\mathbf{e}}_i$, where ${\mathbf{e}}_i$ is each constant Cartesian basis vector. Since $\nabla \cdot (f{\mathbf{e}}_i)={\partial }_{i}f$ and $(f{\mathbf{e}}_i)\cdot \mathbf{n}=f{n}_i$, the three component identities combine to give

For the specified scalar field, $\nabla f=(z,0,x)$. The $x$ term integrates to zero by reflection symmetry. At height $z$, the region has a disk cross-section of area $\pi (a^2-cz)$, and $0\le z\le a^2/c$. Therefore the volume integral isTo check the surface integral independently, split the boundary into the base disk and the curved graph $z=g(x,y)=(a^2-x^2-y^2)/c$. On the base $f=0$, so its contribution vanishes. The region lies below the curved graph, whose outward vector surface element of a graph isOver the disk $x^2+y^2\le a^2$, the $y$ and $z$ components of $fd\mathbf{S}$ are odd in $x$ and integrate to zero. Its $x$ component gives, in polar coordinates,Thus the two explicitly evaluated integrals agree, with the outward orientation included.