$\newcommand{\angles}[1]{\left\langle\, #1 \,\right\rangle}
\newcommand{\braces}[1]{\left\lbrace\, #1 \,\right\rbrace}
\newcommand{\bracks}[1]{\left\lbrack\, #1 \,\right\rbrack}
\newcommand{\ceil}[1]{\,\left\lceil\, #1 \,\right\rceil\,}
\newcommand{\dd}{{\rm d}}
\newcommand{\ds}[1]{\displaystyle{#1}}
\newcommand{\expo}[1]{\,{\rm e}^{#1}\,}
\newcommand{\fermi}{\,{\rm f}}
\newcommand{\floor}[1]{\,\left\lfloor #1 \right\rfloor\,}
\newcommand{\half}{{1 \over 2}}
\newcommand{\ic}{{\rm i}}
\newcommand{\iff}{\Longleftrightarrow}
\newcommand{\imp}{\Longrightarrow}
\newcommand{\pars}[1]{\left(\, #1 \,\right)}
\newcommand{\partiald}[3][]{\frac{\partial^{#1} #2}{\partial #3^{#1}}}
\newcommand{\pp}{{\cal P}}
\newcommand{\root}[2][]{\,\sqrt[#1]{\vphantom{\large A}\,#2\,}\,}
\newcommand{\sech}{\,{\rm sech}}
\newcommand{\sgn}{\,{\rm sgn}}
\newcommand{\totald}[3][]{\frac{{\rm d}^{#1} #2}{{\rm d} #3^{#1}}}
\newcommand{\verts}[1]{\left\vert\, #1 \,\right\vert}$
\begin{align}
&\color{#66f}{\large\int_{0}^{a}\int_{0}^{x}\int_{0}^{y}\int_{0}^{z}
\fermi\pars{u}\,\dd u\,\dd z \,\dd y\,\dd x}
\\[3mm]&=\int_{0}^{a}\int_{0}^{a}\int_{0}^{a}\int_{0}^{a}
\fermi\pars{u}\Theta\pars{z - u}\Theta\pars{y - z}\Theta\pars{x - y}
\,\dd u\,\dd z\,\dd y\,\dd x
\\[5mm]&=\int_{0}^{a}\int_{0}^{a}\int_{0}^{a}
\fermi\pars{u}\Theta\pars{z - u}\Theta\pars{y - z}\
\overbrace{\bracks{\int_{0}^{a}\Theta\pars{x - y}\,\dd x}}
^{\ds{=\ \color{#c00000}{a - y}}}\ \,\dd y\,\dd z\,\dd u
\\[3mm]&=-\int_{0}^{a}\int_{0}^{a}
\fermi\pars{u}\Theta\pars{z - u}\bracks{\int_{0}^{a}\pars{y - a}
\Theta\pars{y - z}\,\dd y}\,\dd z\,\dd u
\\[3mm]&=-\int_{0}^{a}\int_{0}^{a}
\fermi\pars{u}\Theta\pars{z - u}\bracks{\int_{z}^{a}\pars{y - a}
\,\dd y}\,\dd z\,\dd u
\\[3mm]&=\int_{0}^{a}\int_{0}^{a}
\fermi\pars{u}\Theta\pars{z - u}{\pars{z - a}^{2} \over 2}\,\dd z\,\dd u
=\int_{0}^{a}\fermi\pars{u}\int_{u}^{a}{\pars{z - a}^{2} \over 2}\,\dd z\,\dd u
\\[3mm]&=-\int_{0}^{a}\fermi\pars{u}{\pars{u - a}^{3} \over 3!}\,\dd u
=\color{#66f}{\large\int_{0}^{a}{\pars{a - u}^{3} \over 3!}\fermi\pars{u}\,\dd u}
\end{align}
$\ds{\Theta\pars{\xi}}$ is the Heaviside Step Function.