Considering the change of cooordinates \begin{align*}
\left\{\begin{array}{l}
x=a\,r\sin\vartheta\cos\varphi\\
y=b\,r\sin\vartheta\sin\varphi\\
z=c\,r\cos\vartheta
\end{array}\right\}\,,\quad r\in[0,+\infty)\,,\; \vartheta\in[0,\pi], \; \varphi\in[0,2\pi]\,,
\end{align*} with Jacobian
\begin{align*}
\biggl|{\frac{\partial(x,y,z)}{\partial(r,\vartheta,\varphi)}}\biggr|&=\left|{\begin{array}{ccc}
a\sin\vartheta\cos\varphi & a\,r\cos\vartheta\cos\varphi & -a\,r\sin\vartheta\sin\varphi\\
b\sin\vartheta\sin\varphi & b\,r\cos\vartheta\sin\varphi & b\,r\sin\vartheta\cos\varphi\\
c\cos\vartheta & -c\,r\sin\vartheta & 0
\end{array}
}\right|\\
&=abc\,r^2\sin\vartheta\,,
\end{align*}
- the solid ellipsoid ${\rm{E}}=\Big\{{(x,y,z)\in{\mathbb{R}}^3\;\big|\;\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}\leqslant1\,}\Big\}$ can be represented as \[{\rm{E}}=\Big\{{\big(a\,r\sin\vartheta\cos\varphi,\, b\,r\sin\vartheta\sin\varphi,\,\,c\,r\cos\vartheta\big)\in{\mathbb{R}}^3\;\big|\;r\in[0,1]\,,\; \vartheta\in[0,\pi], \; \varphi\in[0,2\pi]\,}\Big\}.\]
So \begin{align*}
\mathop{\iiint}\limits_{\rm{E}}{xyz\,d(x,y,z)}&=\mathop{\iiint}\limits_{\rm{E}}{abc\,r^3\sin^2\vartheta\cos\vartheta\cos\varphi\sin\varphi\,\big|abc\,r^2\sin\vartheta\big|\,d(r,\vartheta,\varphi)}\\
&=(abc)^2\int_{0}^{1}\int_{0}^{\pi}\int_{0}^{2\pi}{r^5\sin^3\vartheta\cos\vartheta\cos\varphi\sin\varphi\,d\varphi\,d\vartheta\,dr}\\
&=\frac{(abc)^2}{2}\int_{0}^{1}r^5\int_{0}^{\pi}\sin^3\vartheta\cos\vartheta\cancelto{0}{\bigg(\int_{0}^{2\pi}{\sin(2\varphi)\,d\varphi}\bigg)}\,d\vartheta\,dr\\
&=0\,.\end{align*}
- Similarly the part ${\rm{E}}_{+}=\Big\{{(x,y,z)\in{\mathbb{R}}^3\;\big|\;\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}\leqslant1\,,\; x\geqslant0\,,\; y\geqslant0\,,\; z\geqslant0\,}\Big\}$ of ${\rm{E}}$ can be represented as \[{\rm{E}}_{+}=\Big\{{\big(a\,r\sin\vartheta\cos\varphi,\, b\,r\sin\vartheta\sin\varphi,\,\,c\,r\cos\vartheta\big)\in{\mathbb{R}}^3\;\big|\;r\in[0,1]\,,\; \vartheta\in\big[0,\tfrac{\pi}{2}\big], \; \varphi\in\big[0,\tfrac{\pi}{2}\big]\,}\Big\}.\] So
\begin{align*}
\mathop{\iiint}\limits_{{\rm{E}}_{+}}{xyz\,d(x,y,z)}&=\mathop{\iiint}\limits_{{\rm{E}}_+}{abc\,r^3\sin^2\vartheta\cos\vartheta\cos\varphi\sin\varphi\,\big|abc\,r^2\sin\vartheta\big|\,d(r,\vartheta,\varphi)}\\
&=(abc)^2\int_{0}^{1}\int_{0}^{\frac{\pi}{2}}\int_{0}^{\frac{\pi}{2}}{r^5\sin^3\vartheta\cos\vartheta\cos\varphi\sin\varphi\,d\varphi\,d\vartheta\,dr}\\
&=\frac{(abc)^2}{2}\int_{0}^{1}r^5\int_{0}^{\frac{\pi}{2}}\sin^3\vartheta\cos\vartheta\cancelto{1}{\bigg(\int_{0}^{\frac{\pi}{2}}{\sin(2\varphi)\,d\varphi}\bigg)}\,d\vartheta\,dr\\
&=\frac{(abc)^2}{2}\int_{0}^{1}r^5\cancelto{\frac{1}{4}}{\bigg(\int_{0}^{\frac{\pi}{2}}\sin^3\vartheta\cos\vartheta\,d\vartheta\bigg)}\,dr\\
&=\frac{(abc)^2}{8}\int_{0}^{1}r^5\,dr\\
&=\frac{(abc)^2}{8}\,\frac{1}{6}\\
&=\frac{(abc)^2}{48}\,.\end{align*}
Additional question: For $n\in\mathbb{N}$ evaluate
\[\displaystyle\mathop{\iiint}\limits_{\rm{E}}\big(a^2b^2-b^2x^2-a^2y^2\big)^{n-\frac{1}{2}}\,d(x,y,z)\,,\] where ${\rm{E}}$ is the above solid ellipsoid.