Use the Divergence Theorem to calculate the surface integral ∫S​F⋅dS; that is, calculate the flux of F across S.
F(x,y,z)=3xy²i+xe²j+z³k
S is the surface of the solid bounded by the cylinder.