Let $\mathbf F$ be the vector field $\mathbf F(x, y, z) = (xz, yz, -(x^2 + y^2))$, and let the region $\mathscr H$ be delimited by the cone
$x^2 + y^2 = z^2$, and the (multivalued) function $z = \arctan (y / x)$, with $0 \leq z \leq 2\pi$. Verify Stokes' theorem for this scenario by showing
$$
\iint\limits_{\mathscr H} (\nabla \times \mathbf F) \cdot \mathbf n \diff A = \oint\limits_{\partial \mathscr H} \mathbf F \cdot \diff \mathbf s .
$$
Proposed by professor Pablo Álvarez Caudevilla as a challenge problem in the subject Vector Calculus.
We first require the curl of $\mathbf F$, which follows from a cross product with the nabla "vector": $$ \nabla \times \mathbf F = \begin{vmatrix} \ihat & \jhat & \khat\\ \partial_x & \partial_y & \partial_z\\ xz & yz & -(x^2 + y^2) \end{vmatrix} = -3y \ihat + 3x \jhat = (-3y, 3x, 0), $$ and now we parametrize the surface, letting $x = r\cos\theta$, $y = r\sin\theta$ and therefore $z = \theta$, (notice that this eliminates the multiple branches of the inverse tangent, so we can integrate $\theta$ from $0$ to $2\pi$ without any issues.) The limits of integration are $0 \leq \theta \leq 2\pi$, and $0 \leq r \leq \theta$. Under this parametrization, the direction vectors become $$ \mathbf T_r = (\cos \theta, \sin \theta, 0),\qquad \mathbf T_\theta = (-r\sin\theta, r\cos\theta, 1), $$ and consequently the normal vector is $$ \mathbf n = \mathbf T_r \times \mathbf T_\theta = \begin{vmatrix} \ihat & \jhat & \khat\\ \cos\theta & \sin\theta & 0\\ -r\sin\theta & r\cos\theta & 1 \end{vmatrix} = \sin\theta \ihat - \cos \theta \jhat + r\khat = (\sin\theta, -\cos\theta, r); $$ since the $z$-coordinate is $r$, which is always positive, we are integrating in the outer direction. This is what we want, so that the boundary is then integrated along the positive way. Finally, putting everything together we get $$ \begin{align*} \iint\limits_{\mathscr H} (\nabla \times \mathbf F) \cdot \mathbf n \diff A &= \lint_0^{2\pi} \lint_0^r (-3r\sin\theta, 3r\cos\theta, 0) \cdot (\sin\theta, -\cos\theta, r) \diff r \diff\theta\\ &= -3 \lint_0^{2\pi} \lint_0^r r (\cancel{\sin^2\theta + \cos^2\theta}) \diff r \diff \theta = -\frac{3}{2} \lint_0^{2\pi} \theta^2 \diff\theta = \color{green}\boxed{-4\pi^3.} \end{align*} $$