 Solution

The divergence of the position vector field is $\nabla \cdot \mathbf{r}=3$. At height $z$, the horizontal cross-section is a disk of area $\pi (1-z^2)$. By the divergence theorem, the outward flux integral is three times the enclosed volume:ThereforeA direct check separates the two pieces of the surface integral. On the unit sphere, the outward unit normal is $\mathbf{r}$, so $\mathbf{r}\cdot \mathbf{n}=1$ and the curved contribution is $2\pi (z_0+1)$. On the top disk, $\mathbf{n}={\mathbf{e}}_z$, giving $\pi z_0(1-z_0^2)$. Their sum is the same flux.