Stokes’ theorem says that for a open surface $S$ with boundary $c$ positively oriented,
$$\iint_{S} \text{curl} \vec F \cdot d\vec S=\oint_{c} F \cdot d\vec r$$
In our case the boundary of our surface is $c=\{(x,y,z) \in \mathbb{R}^3 : x^2+y^2=1, z=0 \}$.
We may parametrize our surface with $\vec r(\theta)=\langle \cos \theta, \sin \theta,0 \rangle$ and $\theta \in [0,2\pi]$ so that $d \vec r=\langle -\sin \theta, \cos \theta,0 \rangle d\theta$.
On the other hand $\vec F=\langle z^2,x,y^2 \rangle$ and on $c$ we have $\vec F=\langle 0,\cos \theta, \sin^2 \theta \rangle$ as $z=0,x=\cos \theta$, and $y=\sin \theta$.
$$\oint_{c} \vec F \cdot dr$$
$$=\int_{0}^{2\pi} \langle 0,\cos \theta, \sin^2 \theta \rangle \cdot \langle -\sin \theta, \cos \theta, 0 \rangle d\theta$$
$$=\int_{0}^{2\pi} \cos^2 \theta d\theta$$
$$=\pi$$
On the other hand we can directly compute,
$$\iint_{S} \text{curl} F \cdot n dS=\iint_{D} \text{curl} F \cdot n |\vec r_{\theta} \times \vec r_{\phi}| dA$$
Parametrize the surface $S=W$ with,
$$\vec r(\theta,\phi)=\langle \sin \phi \cos \theta, \sin \phi \sin \theta, \cos \phi \rangle$$
And $D=\{(\theta,\phi) : (\theta,\phi) \in [0,2\pi] \times [0,\frac{\pi}{2}] \}$.
On this surface $S$ we have,
$$\text{curl} \vec F=\langle 2\sin \phi \sin \theta, 2 \cos \phi,1 \rangle$$
Also, $|\vec r_{\theta} \times \vec r_{\phi}|=\sin \phi$. (It is is standard practice to show that when a sphere of radius $\rho$ is considered this quantity turns out to be $\rho^2 \sin \phi$).
A normal to a sphere centered at the origin is $\langle x,y,z \rangle$, hence the outward pointing unit normal to a sphere of radius $\rho$ is,
$$\vec n=\frac{\langle x,y,z \rangle}{\rho}$$
Since $\rho=\sqrt{x^2+y^2+z^2}$. In our case $\rho=1$, so the outward unit normal is,
$$\vec n=\langle x,y,z \rangle$$
$$=\langle \sin \phi \cos \theta, \sin \phi \sin \theta, \cos \phi \rangle$$
Putting everything together,
$$\text{curl} F \cdot \vec n |\vec r_{\theta} \times \vec r_{\phi}|=2\sin^3 \phi \sin \theta \cos \theta+2 \sin^2 \phi \cos \phi \sin \theta+\sin \phi \cos \phi$$
And the surface integral we want to compute is equal to,
$$\int_{0}^{2\pi} \int_{0}^{\frac{\pi}{2}} \left(2\sin^3 \phi \sin \theta \cos \theta+2 \sin^2 \phi \cos \phi \sin \theta+\sin \phi \cos \phi \right) d\phi d\theta$$
We can safely ignore the term with $\sin \theta$ being the only part with $\theta$ because $\int_{c}^{d} \int_{a}^{b} f(x)g(y) dx dy=\left(\int_{c}^{d} g(y) dy \right) \left( \int_{a}^{b} f(x) dx \right)$ and $\int_{0}^{2\pi} \sin \theta d\theta=0$. Similarly we can ignore the term with $\sin \theta \cos \theta$ being the part with $\theta$ for pretty much the same reason. We end with,
$$\int_{0}^{2\pi} \int_{0}^{\frac{\pi}{2}} \sin \phi \cos \phi d\phi d\theta$$
$$=2\pi \int_{0}^{\frac{\pi}{2}} \sin \phi \cos \phi d\phi$$
$$=\pi \int_{0}^{\frac{\pi}{2}} \sin (2 \phi) d\phi$$
$$=\pi$$
