Consider the vector field $\mathbb{v} = \operatorname{curl}\mathbb{u}$, where $\mathbb{u} = (xy, xz^2, x^2y)$.
Find the flux of $\mathbb{v}$ downwards (negative z-component) through the paraboloid $z = 1 - x^2 - y^2$, where $x \geq 0$ and $z \geq 0$.
So, $\mathbb{v} = (x^2-2xz, -2xy, z^2-x)$. Tried then apply the divergence theorem, but found out that it probably wouldn't work as the surface is open.
Kind of lost on where to go from here. I think the next step would be to just compute $$\iint_{S} \mathbb{v}\cdot\mathbb{n}\;dS,$$ but I have to parametrize the surface first? No idea how to do that, or if I'm going in the right direction.

$\mathbb{v} = (x^2-2xz, -xy, z^2-x)$.
Then, when $z=0\implies(x^2, -2xy, -x)$.
\begin{align} \iint_S(x^2, -2xy, -x)\cdot(0, 0, -1)\;dS &= \iint_Sx\;dS \\ &= \int_0^1\int_{-\pi/2}^{\pi/2} (r\cos{\theta})r\;d\theta dr \\ &= \int_0^1 2r^2 \;d\theta dr \\ &= 2/3. \end{align}
The integral for the other wall also becomes zero like the divergence of the $\operatorname{curl}\mathbb{u}$. This leaves us with the answer $2/3$.
Check out @giobrach's answer to see some more details to the solution.