The value of $\iint_S \vec{F} \cdot \vec{N} \, ds$ where $\vec{F} = 2x^2y \hat{i} - y^2 \hat{j} + 4xz^2 \hat{k}$ and $S$ is the closed surface of the region in the first octant bounded by the cylinder $y^2 + z^2 = 9$ and the planes $x = 0, x = 2, y = 0, z = 0$, is: