Evaluate $1$-loop contribution to the $4$-point Green's function I am trying to evaluate the following integral
\begin{equation}
I = \int \frac{d^d p_\text{E}}{(2 \pi)^d} \frac{1}{(p_\text{E}^2+m^2)((q_\text{E}-p_\text{E})^2 + m^2)} \tag{1}
\end{equation}
where $p_\text{E}$ is the momentum in Euclidean space (i.e. after having performed a Wick rotation). Using Feynman's parametrization trick, and some effort, equation $(1)$ can be written as:
\begin{equation}
I = \int\limits^1_0 \mathrm{d} x \int \frac{d^d p'_\text{E}}{(2 \pi)^d} \frac{1}{\left[p_\text{E}'^2 + q_\text{E}^2 x(1-x) + m^2 \right]^2} \tag{2}
\end{equation}
where $p'_\text{E} = p_\text{E} - q(1-x)$. I know that equation $(2)$ is correct as I have verified by checking 2 different sources. Now, according to the two aforementioned sources, equation $(2)$  can be written evaluated as:
\begin{equation}
I = \frac{\pi^{d/2}}{(2 \pi)^d} \Gamma(2-d/2) \int\limits^1_0 \mathrm{d} x \; (q_\text{E}^2 x(1-x) + m^2)^{d/2-2} \tag{3}
\end{equation}
where $\Gamma(2-d/2)$ is the gamma function. However, I don't know where this comes from.
So far, the closest I have come is by writing equation $(2)$ as:
\begin{equation}
I = \int\limits^1_0 \mathrm{d} x \int \frac{d^d p'_\text{E}}{(2 \pi)^d} \int\limits^\infty_0 \mathrm{d}u \int\limits^\infty_0 \mathrm{d}u \; e^{-u(p_\text{E}'^2 + q_\text{E}^2 x(1-x) + m^2 )}  \tag{4}
\end{equation}
We can perform the integration over $p'_\text{E}$ by using Gaussian integration:
\begin{equation}
\int \mathrm{d}^d p_{\text{E}} \;  e^{-up'^2_{\text{E}}} = \left(\frac{\pi}{u}\right)^{d/2}
\end{equation}
Thus equation $(4)$ becomes:
\begin{equation}
I = \frac{\pi^{d/2}}{(2 \pi)^d} \int\limits^1_0 \mathrm{d} x \int\limits^\infty_0 \mathrm{d}u \int\limits^\infty_0 \mathrm{d}u \; u^{-d/2}e^{-u(q_\text{E}^2 x(1-x) + m^2 )}  \tag{5}
\end{equation}
Now, letting:
\begin{equation}
w = \left(q_\text{E}^2 x(1-x) + m^2 \right) u \Rightarrow \mathrm{d} u = \left(q_\text{E}^2 x(1-x) + m^2 \right)^{-1} \; \mathrm{d} w
\end{equation}
and substituting the above into equation $(5)$ yields:
\begin{equation}
\begin{aligned}
I & = \frac{\pi^{d/2}}{(2 \pi)^d} \int\limits^1_0 \mathrm{d} x \int\limits^\infty_0 \mathrm{d}w \int\limits^\infty_0 \mathrm{d}w \; w^{-d/2}e^{-w} \left(q_\text{E}^2 x(1-x) + m^2 \right)^{d/2-2} \\&
= \frac{\pi^{d/2}}{(2 \pi)^d} \int\limits^1_0 \mathrm{d} x \int\limits^\infty_0 \mathrm{d}w \; \Gamma(1-d/2) \left(q_\text{E}^2 x(1-x) + m^2 \right)^{d/2-2} \tag{6}
\end{aligned}
\end{equation}
At this point I'm stuck. I have no idea if the steps taken to get to equation $(6)$ are correct, but equation $(6)$ does kind of look like equation $(3)$ so I hope they are. Could anybody help me and tell me how I can derive equation $(3)$?
 A: You can prove a general formula
$$
  \int\frac{d^dp_{\mathrm{E}}}{(2\pi)^d}
  \frac{(p_{\mathrm{E}}^2)^m}{(p_{\mathrm{E}}^2+\Delta)^n}=
  \frac{1}{(4\pi)^{d/2}}\frac{\Gamma(m+d/2)\Gamma(n-m-d/2)}{\Gamma(d/2)\Gamma(n)}
  \left(\frac{1}{\Delta}\right)^{n-m-d/2},\quad n>m+d/2
$$
by using Gaussian integral and Euler integral of the first kind.
A: Although the answer given by soliton is sufficient, I've found a way to explicitly evaluate this integral (in case anybody might be interested). Let us start from equation $(2)$ in the original message:
\begin{equation}
I = \int\limits^1_0 \mathrm{d} x \int \frac{d^d p'_\text{E}}{(2 \pi)^d} \frac{1}{\left[p_\text{E}'^2 + q_\text{E}^2 x(1-x) + m^2 \right]^2} \tag{2}
\end{equation}
This can be written as:
\begin{equation}
I = \int\limits_0^1 \mathrm{d}x \int \frac{\mathrm{d}^d p_{\text{E}}'}{(2\pi)^d} \int\limits^\infty_0 \mathrm{d} u_1 \; e^{- u_1 \left( p_{\text{E}}'^2 + q_\text{E}^2 x(1-x) + m^2 \right)}  \int\limits^\infty_0 \mathrm{d} u_2 \; e^{- u_2 \left( p_{\text{E}}'^2 +q_\text{E}^2 x(1-x) + m^2 \right)} 
\end{equation}
The integration over $p_{\text{E}}'$ becomes an ordinary Gaussian integral:
\begin{equation}
\int \mathrm{d}^d p_{\text{E}}' \; e^{-p_{\text{E}}'^2(u_1+u_2)} = \left( \frac{\pi}{u_1+u_2} \right)^{d/2}
\end{equation}
and so:
\begin{equation}
I = \frac{1}{(4 \pi)^{d/2}} \int\limits_0^1 \mathrm{d}x \int\limits^\infty_0 \mathrm{d} u_1 \int\limits^\infty_0 \mathrm{d} u_2 \; \left(\frac{1}{u_1+u_2}\right)^{d/2} e^{- \left( q_\text{E}^2 x(1-x) + m^2 \right)\left(u_1+u_2 \right)} \tag{3}
\end{equation}
We can make the following substitution:
\begin{equation}
\begin{array}{cc}
u_1 = sx \; ,& u_2 = s(1-x)
\end{array}
\end{equation}
where $0 < s < \infty$ (and of course $0 < x < 1$), such that:
\begin{equation}
u_1+u_2=s \; \; \; \text{and} \; \; \; J = \left|\frac{\partial(u_1,u_2)}{\partial(x,s)} \right| = s
\end{equation}
Thus, equation $(3)$ becomes:
\begin{equation}
I = \frac{1}{(4 \pi)^{d/2}} \int\limits_0^1 \mathrm{d}x \int\limits^\infty_0 \mathrm{d} s \; s^{1-d/2} e^{- \left( q_\text{E}^2 x(1-x) + m^2 \right)s}
\end{equation}
Now, we make the substitution:
\begin{equation}
w = \left( q_\text{E}^2 x(1-x) + m^2 \right)s \Rightarrow \mathrm{d} s = \left( q_\text{E}^2 x(1-x) + m^2 \right)^{-1} \mathrm{d} w
\end{equation}
such that:
\begin{equation}
\begin{aligned}
I & = \frac{1}{(4 \pi)^{d/2}} \int\limits_0^1 \mathrm{d}x \int\limits^\infty_0 \mathrm{d} w \; w^{1-d/2} e^{- w} \left( q_\text{E}^2 x(1-x) + m^2 \right)^{d/2-2} \\& 
= \frac{1}{(4 \pi)^{d/2}} \Gamma(2-d/2) \int\limits_0^1 \mathrm{d}x \;  \left( q_\text{E}^2 x(1-x) + m^2 \right)^{d/2-2}
\end{aligned}
\end{equation}
