I was messing around this morning and I discovered the following, which seemed cute enough to share. In this post, I’ll make what strikes me as a very reasonable attempt to define for not an integer. Will I get the function? Wait and see!

We have . So, by basic complex analysis, , where the integral is taken along a loop around the origin. This formula is also morally right for a negative integer: wants to be $\infty$ when (because , so should be infinity, and likewise for the other negative integers). So wants to be zero for and, sure enough, this integral has no poles and vanishes in that case.

We can’t use this formula for not an integer, because has a branch cut and the path of integration would have to cross it. But we can fix that by taking the branch cut of to be along the negative real axis, and drawing our loop out to stretch very far in the negative real direction. Then will be very small at the point where the integration path crosses the real axis, so the branch cut will contribute very little. In the limit, we can define

where is a path that comes in from the negative real direction below the real axis, circles around the origin, and returns to infinity in the negative real direction above the axis. This integral will converge for all complex

So, how does this do as a definition of ? Well, it obeys the right recursion. A quick integration by parts gives , so .

Let’s take our path and shrink it towards the negative real axis. As we approach from above (for a positive real), approaches . As we approach from below, approaches . The difference between the two is . So one might think that our integral was equal to .

If you are more careful, you’ll see that this argument only works for ; otherwise, the pole at the origin is too wild to permit the limiting process. So we get that our previous definition is equivalent to

for .

This is where a person who has seen the function defined before will say “well, you’re on the right track, but that sure looks funky.” Writing for the standard complex extension of the factorial function^{1}, we have . So I’ve got the right integral, but it’s being evaluated at the wrong place, and there is this strange extra factor of floating around.

But it all works out! We have the functional equation of the function:

So the integral I have above really is the standard extension, but gotten at from the other side.

One wants to turn this into a proof of the functional equation, but as yet I don’t see how…

^{1}For historical reasons, . So I’m writing for the function .

Small correction.

“so (-1)! should be zero, and likewise for the other negative integers” should read “so (-1)! should be infinity…”

Fixed, thanks.

Thank you David. I am curious how you came to working this out. What were you “messing around”, or what was your frame of mind, or mood? (I don’t really know what to ask but any answer will be welcome.)

In any case it was very interesting. I wonder whether it is possible rather to _not_ arrive at the Gamma function, given any “reasonably convergent” formula that matches n! on integers.

I was basically messing around. I noticed that the integrand in looked a lot like the integrand defining the function, but the contours didn’t match up, so I was trying to get them to match up.

I was also rocking my daughter to sleep for the first 20 minutes or so, so I was trying to do this in my head and there was a lot of going through the same computations over and over to catch missing ‘s and ‘s. When doing this sort of routine, hand-occupying child care, my main ways of occupying my mind are doing mental computations and reciting Dylan lyrics to myself.

Thanks. :) This should deter me from asking next time, for fear of catching you doing something embarassing.

Thanks for the interesting post. for that first e^z expression should it read as sum z^n/n! (there’s an a _n)

Fixed, thanks!