I've asked this question in every math class where the teacher has introduced the Gamma function, and never gotten a satisfactory answer. Not only does it seem more natural to extend the factorial directly, but the integral definition $\Gamma(z) = \int_0^\infty  t^{z-1} e^{-t}\,dt$, makes more sense as $\Pi(z) = \int_0^\infty  t^{z} e^{-t}\,dt$. Indeed Wikipedia says that this function was [introduced by Gauss][1], but doesn't explain why it was supplanted by the Gamma function. As that section of the Wikipedia article demonstrates, it also makes its functional equations simpler: we get $\Pi(z) \; \Pi(-z) = \frac{\pi z}{\sin( \pi z)} = \frac{1}{\operatorname{sinc}(z)}$ instead of $\Gamma(1-z) \; \Gamma(z) = \frac{\pi}{\sin{(\pi z)}}$, the duplication formula is simpler, and the Riemann zeta functional equation reduces from $\zeta(s) = 2^s\pi^{s-1}\ \sin\left(\frac{\pi s}{2}\right)\ \Gamma(1-s)\ \zeta(1-s)$ to $\zeta(s) = 2^s\pi^{s-1}\ \sin\left(\frac{\pi s}{2}\right)\ \Pi(-s)\ \zeta(1-s)$.

I suspect that it's just a historical coincidence, in the same way $\pi$ is defined as circumfrence/diameter instead of the much more natural circumfrence/radius. Does anyone have an actual reason why it's better to use $\Gamma(z)$ instead of $\Pi(z)$?


  [1]: http://en.wikipedia.org/wiki/Gamma_function#Pi_function