The proof of all the above formulae are analogous to the classical product rule. As an example we give a proof of the first. By definition of the partial derivative we need to compute the limit of the difference quotient
\begin{equation*}
D(t):=\frac{1}{t} \Bigl(u(\vect x+t\vect e_i)\vect f(\vect x+t\vect e_i) -u(\vect x)\vect f(\vect x)\Bigr)
\end{equation*}
as \(t\to 0\text{,}\) where \(\vect e_i\) is the vector
\begin{equation*}
\vect e_i:=(0,\dots,0,1,0,\dots,0)\text{,}
\end{equation*}
where the \(1\) appears in the \(i\)-th coordinate
By adding and subtracting a term we can rewrite \(D(t)\) by
\begin{equation*}
D(t)=\frac{1}{t} \Bigl(u(\vect x+t\vect e_i)-u(\vect x)\Bigr)\vect f(\vect x+t\vect e_i) +u(x)\frac{1}{t} \Bigl(\vect f(\vect x+t\vect e_i)-\vect f(\vect x)\Bigr)\text{.}
\end{equation*}
Using the continuity of the functions and the definition of the partial derivatives we get from the above that
\begin{equation*}
\lim_{t\to 0}D(t)
=\Bigl(\frac{\partial u}{\partial x_i}(\vect x)\Bigr)
\vect f(\vect x)
+u(\vect x)\frac{\partial \vect f}{\partial x_i}(\vect x)
\end{equation*}
as required.