- Find the flux of the vector field $\vec{G}=\operatorname{curl} \vec{F}$, where
$$
\vec{F}(x, y, z)=y^{3} \vec{i}+x^{3} \vec{j}+z^{3} \vec{k}
$$
through the upward oriented part of paraboloid
$$
z=1-x^{2}-y^{2}
$$
lying over the plane $z=0$.
What i evaluated is

Is that the right way to evaluate
Gauss's Law aka the Divergence Theorem states that the volume integral of the divergence of a vector field equals the total flux of that vector field through the entire closed boundary surface of the volume.
You have identified a volume between the plane $z=0$ and the function $z = 1-x^2 - y^2$. That volume has two pieces to the boundary: the flat circle within the plane $z=0$, and the piece of the paraboloid above it. So the outward (downward) flux through the circle plus the flux outward (upward) through the paraboloid equals the volume integral of the divergence, which you have already found to be $0$.
Your goal is the upward flux through the paraboloid, so find the downward flux through the circle and reverse the sign to get your answer.
For the flat circle, the unit normal vector is $-\hat z$, so $\overrightarrow F \cdot \hat n = <y^3,x^3, z^3>\cdot<0,0,-1> = -z^3$. Since you are in the plane $z=0$, $\overrightarrow F \cdot \hat n = 0$ and the flux through the circle is $0$, meaning that the flux through the paraboloid is also $0$.