Solution (source code)

= Solution

The specified unit <normal vector> is $\mathbf e_z$, while $F_z=0$ everywhere. Thus the normal component of the field vanishes pointwise, giving <flux>
$$
\boxed{\int_S\mathbf F\cdot d\mathbf S=0}.
$$