The functions $f(x)$ and $g(x)$ have Laplace transforms $F(p)$ and $G(p)$ respectively, and $f(x)=g(x)=0$ for $x \leqslant 0$. The convolution $h(x)$ of $f(x)$ and $g(x)$ is defined by

$$h(x)=\int_{0}^{x} f(y) g(x-y) d y \quad \text { for } \quad x>0 \quad \text { and } \quad h(x)=0 \quad \text { for } \quad x \leqslant 0$$

Express the Laplace transform $H(p)$ of $h(x)$ in terms of $F(p)$ and $G(p)$.

Now suppose that $f(x)=x^{\alpha}$ and $g(x)=x^{\beta}$ for $x>0$, where $\alpha, \beta>-1$. Find expressions for $F(p)$ and $G(p)$ by using a standard integral formula for the Gamma function. Find an expression for $h(x)$ by using a standard integral formula for the Beta function. Hence deduce that

$$\frac{\Gamma(z) \Gamma(w)}{\Gamma(z+w)}=\mathrm{B}(z, w)$$

for all $\operatorname{Re}(z)>0, \operatorname{Re}(w)>0$.