La operación 𝑥! es una tan atada a los números naturales que cuesta imaginarla expandida a dominios más allá de ellos.
No obstante, ello es posible utilizando la función gamma, definida como

Γ(𝑛)=0𝑥𝑛1𝑒𝑥d𝑥

Deducirla es un proceso tardado, iniciamos con la integral de ln(𝑥) aplicando integración por partes (por conveniencia y claridad, omitiremos las constantes de integración)

ln(𝑥)d𝑥=𝑥ln(𝑥)𝑑𝑥ln(𝑥)d𝑥=𝑥ln(𝑥)1𝑥

Ahora observemos lo que sucede al integrar ln2(𝑥), poniendo atención en el último término

ln2(𝑥)d𝑥=𝑥ln2(𝑥)2ln(𝑥)𝑑𝑥ln2(𝑥)d𝑥=𝑥ln2(𝑥)2(𝑥ln(𝑥)𝑥)ln2(𝑥)d𝑥=𝑥ln2(𝑥)2𝑥ln(𝑥)+21𝑥

Y al integrar ln3(𝑥)

ln3(𝑥)d𝑥=𝑥ln3(𝑥)3ln2(𝑥)𝑑𝑥ln3(𝑥)d𝑥=𝑥ln3(𝑥)3(𝑥ln2(𝑥)2𝑥ln(𝑥)+2𝑥)ln3(𝑥)d𝑥=𝑥ln3(𝑥)3𝑥ln2(𝑥)+6𝑥ln(𝑥)321𝑥

El patrón es claro para el último término, en donde se comienzan a dar vestigios del 𝑥! en el horizonte. Es posible armar una forma general para la integral de ln𝑛(𝑥)

ln𝑛(𝑥)d𝑥=𝑥ln𝑛(𝑥)𝑛ln𝑛1(𝑥)𝑑𝑥=𝑥ln𝑛(𝑥)𝑛ln𝑛1(𝑥)d𝑥=𝑥ln𝑛(𝑥)𝑛(𝑥ln𝑛1(𝑥)(𝑛1)ln𝑛2(𝑥)d𝑥)==𝑥ln𝑛(𝑥)𝑛𝑥ln𝑛1(𝑥)+Σ𝑛𝑘=2(1)𝑛(𝑛!𝑥ln𝑛𝑘(𝑥))

Nótese que el último término, cuando 𝑘=𝑛, está libre de logaritmos, pues el exponente del logaritmo se hace 0.

Ahora investiguemos lo que sucede con la integral definida 10ln𝑛(𝑥)d𝑥. Para ello usaremos la forma general que acabamos de deducir, comencemos evaluando en 𝑥=0

(0)ln𝑛(0)𝑛(0)ln𝑛1(0)+Σ𝑛𝑘=2(1)𝑛(𝑛!(0)ln𝑛𝑘(0))

Rápidamente nos encontramos con el límite

lim𝑥0𝑥ln(𝑥)=lim𝑥0ln(𝑥)1𝑥=lim𝑥01𝑥1𝑥2=lim𝑥011𝑥=lim𝑥0𝑥=0

Para los logaritmos con exponente 𝑛>1 es necesario aplicar la regla de L’Hôpital 𝑛 veces para llegar al límite mostrado anteriormente.

lim𝑥0𝑥ln𝑛(𝑥)=lim𝑥0ln𝑛(𝑥)1𝑥=lim𝑥0𝑛ln𝑛1(𝑥)1𝑥1𝑥2=lim𝑥0𝑛ln𝑛1(𝑥)1𝑥=lim𝑥0𝑛(𝑛1)ln𝑛2(𝑥)1𝑥1𝑥2=lim𝑥0𝑛(𝑛1)ln𝑛2(𝑥)1𝑥==lim𝑥0(1)𝑛𝑛!1𝑥=lim𝑥0(1)𝑛(𝑛!𝑥)=0

De esta manera, sabemos que si sustituimos 𝑥=0 entonces se cumple que términos con 𝑥ln𝑛𝑘(𝑥)=0|0𝑘<𝑛. De esta manera, el único término de nuestra forma general que nos queda es el independiente de cualquier logaritmo, es decir (1)𝑛𝑛!𝑥, pero evidentemente eso es 0, así que todo el límite inferior de nuestra integral es también 0. Ahora evaluemos el límite superior con 𝑥=1 usando nuestra forma generalizada

(1)ln𝑛(1)𝑛(1)ln𝑛1(1)+Σ𝑛𝑘=2(1)𝑛(𝑛!(1)ln𝑛𝑘(1))

De manera similar, todos los términos con un ln𝑛𝑘(𝑥) se anulan, pero en esta ocasión el término independiente no lo hace, sino que se evalúa a (1)𝑛𝑛!(1)=(1)𝑛𝑛!. Puesto de forma más clara:

10ln𝑛(𝑥)d𝑥=𝑥ln𝑛(𝑥)𝑛𝑥ln𝑛1(𝑥)+Σ𝑛𝑘=2(1)𝑛(𝑛!𝑥ln𝑛𝑘(𝑥))|10=(1)𝑛𝑛!0=(1)𝑛𝑛!

Pero todavía está ese signo alternante molestando, por lo que debemos despejarlo

𝑛!=(1)𝑛10ln𝑛(𝑥)d𝑥𝑛!=10(1)𝑛ln𝑛(𝑥)d𝑥𝑛!=10(1ln(𝑥))𝑛d𝑥𝑛!=10ln𝑛(𝑥)d𝑥

De esta manera llegamos a la forma inicial encontrada por Euler en su carta a Goldbach, el 13 de octubre de 1729, aunque la forma en que la expresó era ligeramente distinta

𝑛!=10(log𝑠)𝑛d𝑠

Aplicando un cambio de variable 𝑢=ln𝑥 encontramos que

𝑢=ln(1𝑥)𝑒𝑢=1𝑥𝑒𝑢=𝑥𝑒𝑢d𝑢=d𝑥

También se tiene que cuando 𝑥=0𝑢=ln(0)= y cuando 𝑥=1𝑢=ln(1)=0. Entonces, reescribiendo en términos de 𝑢:

𝑛!=0(𝑢𝑛)(𝑒𝑢)d𝑢𝑛!=0(𝑢𝑛)(𝑒𝑢)d𝑢𝑛!=0𝑢𝑛𝑒𝑢d𝑢

Pero esto no es la función gamma, antes tenemos que desplazarla por razones que honestamente no entiendo. De esta forma llegamos a

Γ(𝑛)=0𝑢𝑛1𝑒𝑢d𝑢

De tal forma que

Γ(𝑛+1)=𝑛!

De cabrones está muy chingón poder extender funciones limitadas a los naturales utilizando propiedades emergentes del cálculo, se siente como forzar la mano de la lógica para conseguir el resultado deseado. O, visto de otro modo, ser guiado a ciegas por las matemáticas hasta una salida extraña, abriendo un mundo de posibilidades, pues es importante considerar que esta nueva definición del factorial permite evaluar su comportamiento en los complejos, lo cual siempre es divertido.