Compute the flux of the vector field $\vec F=x\vec i+y\vec j+z\vec k$ through the surface $S$, which is a closed cylinder of radius 1, centered on the $x$-axis, with $-3\le x\le 3$, and oriented outward.

flux =