 Solution

Direct differentiation givesOn the unit disk $D$, polar coordinates give ${\int }_D{y}^{2}dA=\pi /4$, ${\int }_{D}xdA=0$ and ${\int }_{D}dA=\pi$. Consequently the volume integral isFor an independent surface integral calculation, the relevant cross product isThe bottom face contributes ${\int }_{D}xydA=0$ and the top contributes ${\int }_D(-h{y}^2-xy)dA=-\pi h/4$. On the cylindrical side, put $x=\cos \theta$, $y=\sin \theta$, $\hat{\mathbf{n}}=(x,y,0)$ and $dS=d\theta dz$. The flux density is $x^3+2xy{z}^2+2y^{2}z-x{y}^2$. All terms except $2y^{2}z$ integrate to zero because they are odd functions of $x$ or $y$, givingThus minus the total outward flux is $-(\pi h^2-\pi h/4)$, in agreement with the volume integral.