Let E be the solid unit cube with diagonally opposite corners at the origin and (1, 1, 1), and faces parallel to the coordinate planes. Let S be the surface of E, oriented with the outward-pointing normal. Use a CAS to find $\iint_S \mathbf{F} \cdot d\mathbf{S}$ using the divergence theorem if $\mathbf{F}(x, y, z) = 3xy\mathbf{i} + 5ye^x\mathbf{j} + x \sin(z)\mathbf{k}$. (Round your answer to four decimal places.)