**The Gamma Function and Fractal Factorials!**This fractal image by

**Thomas Oléron Evans** was created by using iterations of the

**Gamma function**, which is a continuous version of the factorial function.

If n is a positive integer, the

**factorial** of n, n!, is defined to be the product of all the integers from 1 up to n; for example, 4!=1x2x3x4=24. It is clear from the definition that (n+1)! is the product of n+1 and n!, but it is not immediately clear what the “right” way is to extend the factorial function to non-integer values.

If t is a complex number with a positive real part, the Gamma function Γ(t) is defined by integrating the function x^{t–1}e^{–x} from x=0 to infinity. It is a straightforward exercise using integration by parts and mathematical induction to prove that if n is a positive integer, then Γ(n) is equal to (n–1)!, the factorial of (n–1). Since Γ(1)=1, this gives a justification (there are many others) that the factorial of zero is 1.

Using a technique called analytic continuation, the Gamma function can then be extended to all complex numbers except negative integers and zero. The resulting function, Γ(t), is infinitely differentiable, except at the nonpositive integers, where it has simple poles; the latter are the same kind of singularity that the function f(x)=1/x has at x=0. A particularly nice property of the Gamma function is that it satisfies Γ(t+1)=tΓ(t), which extends the recursive property n!=n(n–1)! satisfied by factorials. It is therefore natural to define the factorial of a complex number z by z!=Γ(z+1).

At first, it may not seem very likely that iterating the complex factorial could produce interesting fractals. If n is an integer that is at least 3, then taking repeated factorials of n will produce a sequence that tends to infinity very quickly. However, if one starts with certain complex numbers, such as 1–i, repeated applications of the complex factorial behave very differently. It turns out that (1–i)! is approximately 0.653–0.343i, and taking factorials five times, we find that (1–i)!!!!! is approximately 0.991–0.003i. This suggests that iterated factorials of 1–i may produce a sequence that converges to 1.

It turns out that if one takes repeated factorials of almost any complex number, we either obtain a sequence that converges to 1 (as in the case of 1–i) or a sequence that diverges to infinity (as in the case of 3). However, it is not possible to take factorials of negative integers, and there are some rare numbers, like z=2, that are solutions of z!=z and do not exhibit either type of behaviour.

By plotting the points that diverge to infinity in one colour, and the points that converge to 1 in a different colour, fractal patterns emerge. The image shown here uses an ad hoc method of colouring points to indicate the rate of convergence or divergence. The points that converge to 1 are coloured from red (fast convergence) to yellow (slow convergence), and the points that diverge to infinity are coloured from green (slow divergence) to blue (fast divergence)

**Relevant links**Thomas Oléron Evans discusses these fractals in detail in a blog post (

http://www.mathistopheles.co.uk/2015/05/14/fractal-factorials/) which contains this image and many others. He (and I) would be interested in knowing if these fractals have been studied before.

The applications of the Gamma function in mathematics are extensive. Wikipedia has much more information about the function here:

http://en.wikipedia.org/wiki/Gamma_functionThis post appears in my

*Mathematics* collection at

https://plus.google.com/collection/8zrhX#mathematics #sciencesunday **Various recent posts by me**Camellia flower:

https://goo.gl/8WNrluHorse chestnut tree:

https://goo.gl/FPCGI3A Curious Property of 82000:

https://goo.gl/1rVg8y