[guided]We need to show that the product satisfies the same kind of estimates as a symbol of order $m+m'$. Define
\begin{align*}
p: U \times \mathbb{R}^n \to \mathbb{C}
\end{align*}
by
\begin{align*}
p(x,\xi) = a(x,\xi)b(x,\xi).
\end{align*}
Since $a$ and $b$ are smooth on $U \times \mathbb{R}^n$, the product $p$ is smooth on $U \times \mathbb{R}^n$. To prove $p \in S^{m+m'}(U \times \mathbb{R}^n)$, we fix a compact set $K \subset U$ and arbitrary multi-indices $\alpha,\beta \in \mathbb{N}_0^n$, and we must prove an estimate of the form
\begin{align*}
|\partial_x^\alpha\partial_\xi^\beta p(x,\xi)|
\leq
C_{K,\alpha,\beta}\langle \xi \rangle^{m+m'-|\beta|}.
\end{align*}
The only point to track is how the $\xi$-derivatives are distributed. The multi-index Leibniz formula gives
\begin{align*}
\partial_x^\alpha\partial_\xi^\beta p
=
\sum_{\alpha_1+\alpha_2=\alpha}
\sum_{\beta_1+\beta_2=\beta}
\binom{\alpha}{\alpha_1}\binom{\beta}{\beta_1}
(\partial_x^{\alpha_1}\partial_\xi^{\beta_1}a)
(\partial_x^{\alpha_2}\partial_\xi^{\beta_2}b).
\end{align*}
This formula is useful because the defining estimates for $S^m_{1,0}$ are stable under arbitrary $x$-derivatives but lose one order for each $\xi$-derivative.
For a fixed summand, the hypothesis $a \in S^m(U \times \mathbb{R}^n)$ gives a constant $A_{K,\alpha_1,\beta_1} > 0$ such that
\begin{align*}
|\partial_x^{\alpha_1}\partial_\xi^{\beta_1}a(x,\xi)|
\leq
A_{K,\alpha_1,\beta_1}\langle \xi \rangle^{m-|\beta_1|}.
\end{align*}
Similarly, the hypothesis $b \in S^{m'}(U \times \mathbb{R}^n)$ gives a constant $B_{K,\alpha_2,\beta_2} > 0$ such that
\begin{align*}
|\partial_x^{\alpha_2}\partial_\xi^{\beta_2}b(x,\xi)|
\leq
B_{K,\alpha_2,\beta_2}\langle \xi \rangle^{m'-|\beta_2|}.
\end{align*}
Multiplying the two estimates gives
\begin{align*}
|(\partial_x^{\alpha_1}\partial_\xi^{\beta_1}a)(x,\xi)(\partial_x^{\alpha_2}\partial_\xi^{\beta_2}b)(x,\xi)|
\leq
A_{K,\alpha_1,\beta_1}B_{K,\alpha_2,\beta_2}
\langle \xi \rangle^{m+m'-|\beta_1|-|\beta_2|}.
\end{align*}
Because $\beta_1+\beta_2=\beta$, we have $|\beta_1|+|\beta_2|=|\beta|$. Hence every summand has exactly the order
\begin{align*}
m+m'-|\beta|.
\end{align*}
There are only finitely many decompositions $\alpha_1+\alpha_2=\alpha$ and $\beta_1+\beta_2=\beta$, so we may define
\begin{align*}
C_{K,\alpha,\beta}
=
\sum_{\alpha_1+\alpha_2=\alpha}
\sum_{\beta_1+\beta_2=\beta}
\binom{\alpha}{\alpha_1}\binom{\beta}{\beta_1}
A_{K,\alpha_1,\beta_1}B_{K,\alpha_2,\beta_2}.
\end{align*}
Adding the finitely many bounded summands then gives
\begin{align*}
|\partial_x^\alpha\partial_\xi^\beta p(x,\xi)|
\leq
C_{K,\alpha,\beta}\langle \xi \rangle^{m+m'-|\beta|}
\end{align*}
for all $(x,\xi) \in K \times \mathbb{R}^n$. This is the defining estimate for $p \in S^{m+m'}(U \times \mathbb{R}^n)$.[/guided]