\begin{equation*}
\iint\limits_{S}\vec{F}\cdot \hat{n}\,ds= \iint\limits_{OABC}\vec{F}\cdot \hat{n}\,ds+\iint\limits_{DEFG}\vec{F}\cdot \hat{n}\,ds
\end{equation*}
\begin{equation*}
\quad +\iint\limits_{OAFG}\vec{F}\cdot \hat{n}\,ds +\iint\limits_{BCDE}\vec{F}\cdot \hat{n}\,ds
\end{equation*}
\begin{equation}
\qquad +\iint\limits_{ABEF}\vec{F}\cdot \hat{n}\,ds+\iint\limits_{OCDG}\vec{F}\cdot \hat{n}\,ds \tag{1.5.1}
\end{equation}
\begin{equation*}
Now,\quad \iint\limits_{OABC}\vec{F}\cdot \hat{n}\,ds= \iint\limits_{OABC}(4xz\hat{i}-y^{2}\hat{j}+yz\hat{k})\cdot (-\hat{k})\,dx\,dy
\end{equation*}
\begin{equation*}
\quad =-\int\limits_{0}^{1}\int\limits_{0}^{1}yz\,dx\,dy=0.\quad [\because z=0 \quad\text{for face OABC}]
\end{equation*}
\begin{equation*}
\iint\limits_{DEFG}\vec{F}\cdot \hat{n}\,ds= \iint\limits_{DEFG}(4xz\hat{i}-y^{2}\hat{j}+yz\hat{k})\cdot (\hat{k})\,dx\,dy
\end{equation*}
\begin{equation*}
=\int\limits_{0}^{1}\int\limits_{0}^{1}yz\,dx\,dy
=\int\limits_{0}^{1}\int\limits_{0}^{1}y.1 \,dx\,dy = \int\limits_{0}^{1}\,dx\int\limits_{0}^{1}y\,dy
\end{equation*}
\begin{equation*}
= \int\limits_{0}^{1}\,dx\left(\frac{y^{2}}{2}\right)_{0}^{1}
=\frac{1}{2}(x)^{1}_{0}=\frac{1}{2}. \hspace{3pt} \text{(as z=1)}
\end{equation*}
\begin{equation*}
\iint\limits_{OAFG}\vec{F}\cdot \hat{n}\,ds
= \iint\limits_{OAFG}(4xz\hat{i}-y^{2}\hat{j}+yz\hat{k})\cdot (-\hat{j})\,dx\,dz
\end{equation*}
\begin{equation*}
=\int\limits_{0}^{1}\int\limits_{0}^{1}y^{2}\,dx\,dz=0. \quad \text{(as y = 0)}
\end{equation*}
\begin{equation*}
\iint\limits_{BCDE}\vec{F}\cdot \hat{n}\,ds
= \iint\limits_{BCDE}(4xz\hat{i}-y^{2}\hat{j}+yz\hat{k})\cdot (\hat{j})\,dx\,dz=\int\limits_{0}^{1}\int\limits_{0}^{1}-y^{2}\,dx\,dz
\end{equation*}
\begin{equation*}
=-\int\limits_{0}^{1}\int\limits_{0}^{1} dxdz
= -\int\limits_{0}^{1}\,dx\int\limits_{0}^{1}\,dz
= -\left(x\right)^{1}_{0}\left(z\right)^{1}_{0}=-1. \hspace{3pt} \text{(as y=1)}
\end{equation*}
\begin{equation*}
\iint\limits_{ABEF}\vec{F}\cdot \hat{n}\,ds= \iint\limits_{ABEF}(4xz\hat{i}-y^{2}\hat{j}+yz\hat{k})\cdot (\hat{i})\,dy\,dz=\int\limits_{0}^{1}\int\limits_{0}^{1}(4xz\,dy\,dz
\end{equation*}
\begin{equation*}
=\int\limits_{0}^{1}\int\limits_{0}^{1} 4.1.z \,dy\,dz = 2. \quad \text{(as x=1)}
\end{equation*}
\begin{equation*}
\iint\limits_{OCDG}\vec{F}\cdot \hat{n}\,ds= \iint\limits_{OCDG}(4xz\hat{i}-y^{2}\hat{j}+yz\hat{k})\cdot (-\hat{i})\,dy\,dz
\end{equation*}
\begin{equation*}
=-\int\limits_{0}^{1}\int\limits_{0}^{1}(4xz\,dy\,dz= 0. \quad \text{(as x=0)}
\end{equation*}
On putting these values in eqn.
(1.5.1), we get -
\begin{equation*}
\iint\limits_{S}\vec{F}\cdot \hat{n}\,ds = 0+\frac{1}{2}+0-1+2+0=\frac{3}{2}. \hspace{3pt} proved.
\end{equation*}