Calculate the flux of the vector field $\vec{F}(x,y,z) = \cos(x^2+y^2) \vec{k}$ through the disk $x^2+y^2 \leq 16$ oriented upward in the plane $z = 3$.

Flux =