I know that with the case
$$\int_0^{\infty}xe^{-ax}dx$$
You just substitute $u=ax$ and use the gamma function after so that it evaluates to $\frac{1}{a^2}$. However this clearly involves switching the bounds. So in this case, if $a$ were a positive imaginary number we would have our upper bound be $u=i\infty$. So I suppose my question is to show why
$$\int_0^{\infty}xe^{-ax}dx$$
Obeys the equality to $\frac{1}{a^2}$ for complex $a$.
$EDIT$
I know how to do Integration by parts and u substitution, that’s not the problem. The problem is what complex analysis backs up the fact that the integral obeys the identity
$$\int_0^{\infty}xe^{-ax}dx = \frac{1}{a^2}$$
For some complex $a$.