Hace unas semanas, publiqué en la cuenta de Instagram del blog un par de integrales cuyos resultados era el producto y el cociente de dos de los números irracionales más conocidos de las matemáticas: pi (trascendente) y phi, el número áureo (algebraico). Hoy vamos a calcularlas, y confirmar así que dichos resultados son correctos, con la ayuda de dos funciones integrales muy interesantes.
Antes de nada, os dejo a continuación la publicación de Instagram que comentaba al comienzo:
(En el momento de escribir esta entrada, acaba de superar los 1500 likes.)
Y, antes de resolverlas, vamos a comentar las dos funciones integrales que vamos a usar en las correspondientes resoluciones. Como seguro que más de uno ya había imaginado, son las denominadas funciones Gamma () y Beta (
) de Euler:
Sobre la función Gamma ya hemos hablado en Gaussianos, y la función Beta está íntimamente relacionada con ella mediante la siguiente propiedad:
Usaremos algunas otras propiedades y expresiones alternativas en la resolución de las dos integrales, comentando cada una de ellas en el momento en el que las necesitemos.
Por cierto, ambas integrales pueden resolverse mediante integración compleja. No son especialmente difíciles mediante ese método, pero se requiere conocer el teorema de los residuos, con toda la teoría de que dicho teorema lleva detrás. Por ello, me he decantado por utilizar las funciones Gamma y Beta.
Si conoces alguna otra forma de resolver estas integrales, háblanos de ella en los comentarios.
Vamos con la primera integral. Se asegura que
En este caso, vamos a comenzar usando una expresión alternativa para la función Beta. Se tiene que:
(A partir de la definición inicial de la Beta, podemos llegar a esta expresión con el cambio de variable . Os dejo los detalles a vosotros, y si algo no os sale podéis preguntarlo en los comentarios.)
Ésta ya se parece un poco más a la nuestra, al menos los límites de integración ya coinciden. Vamos a hacer el cambio de variable en la integral que queremos calcular:
Sustituimos lo que hemos calculado, simplificamos y recolocamos todo:
Como se puede ver fácilmente, la integral que nos ha quedado corresponde con la expresión alternativa de la Beta que hemos presentado antes. Para calcular los parámetros y
, igualamos y resolvemos el sistema:
Tenemos entonces que nuestra integral inicial es, exactamente
Usamos ahora la propiedad que relaciona la Beta y la Gamma:
Y ahora necesitamos un par de propiedades de la Gamma, que ya vimos en el post sobre ella que enlacé en los primeros párrafos. Son éstas:
Tomando , se tiene que en el numerador de nuestra integral podemos aplicar la segunda propiedad de la Gamma, por lo que nos queda lo siguiente:
Después de todo el esfuerzo realizado, sólo nos queda calcular . El cálculo no es inmediato, pero hay varias formas de hacerlo. Por ejemplo, se puede usar que
, como calculamos en esta entrada, y después usar la fórmula para calcular el seno del ángulo mitad (bueno, y alguna que otra manipulación algebraica). Vamos con ello:
Multiplicamos numerador y denominador por , y en el numerador usamos que
en dos ocasiones (recordemos que
es una de las soluciones de la ecuación
):
Por tanto, tenemos que
Con esto, llegamos al resultado buscado para nuestra integral:
Una de las integrales ya ha caído, vayamos ahora a por la otra. Veremos que:
La estrategia va a ser la misma que en el caso anterior: convertir nuestra integral en una expresión que involucre a una Beta (con el mismo cambio de variable, ), pasar a la Gamma con la propiedad que las relaciona, usar las propiedades de la Gamma y calcular la razón trigonométrica que nos quede. Como los cálculos son muy parecidos, me voy a saltar algunos pasos:
¿Cómo calcular ? Pues, en este caso, es bien sencillo. Al ser
un ángulo positivo menor que
, tenemos que
, que es
, cuyo valor es
. Por tanto, ya tenemos el resultado de la segunda integral:
¿Te ha parecido interesante esta entrada? ¿Tienes alguna duda de alguno de los pasos seguidos? ¿Conoces alguna manera más eficiente/rápida/sencilla de hacer alguno de los cálculos realizados? ¿Has conseguido relacionar las expresiones de la Beta o, por el contrario, te has atrancado en algún punto? Por cualquiera de estas razones, o por alguna otra que se te ocurra, te animo a que comentes a continuación. Muchas gracias.
¿Te ha gustado la entrada? Puedes invitarme a un café, Gauss te lo agradecerá 😉
A quien pueda interesar, un programita sencillísimo, escrito en Pari gp (Universidad de Burdeos), cuya App es gratuita y está en Play Store, para encontrar en segundis todas kas ternas pitagóricas hasta donde se precise: for(r=2,30,for(s=1,r-1,if(gcd(r,s)==1&((r%2==0&s%2==1)||(r%2==1&s%2==0)),a=2*r*s;n=r^2-s^2;m=r^2+s^2;print([r,s,a,n,m])))) [2, 1, 4, 3, 5] [3, 2, 12, 5, 13] [4, 1, 8, 15, 17] [4, 3, 24, 7, 25] [5, 2, 20, 21, 29] [5, 4, 40, 9, 41] [6, 1, 12, 35, 37] [6, 5, 60, 11, 61] [7, 2, 28, 45, 53] [7, 4, 56, 33, 65] [7, 6, 84, 13, 85] [8, 1, 16, 63, 65] [8, 3, 48, 55,… Lee más »
Triángulo rectángulo de lados racionales y área 65, con las fracciones con mayor número de dígitos hasta ahora halladas: [212953601/9673620, 1257570600/212953601, 46952601370867201/2060032213705620, 65] He encontrado 10 soluciones para una área de s = 65, con un móvil barato de 150 €, de hace 2 años. Pienso que hay infinitos triángulos racionales para algunas áreas dadas, aunque son escasos. Un infinito escaso. Habrá alguno, con fracciones de miles de dígitos en numerador y denominador. Y es posible que para un número también infinito de valores enteros del área, no exista ningún triángulo racional. Pero no sé cómo demostrarlo. Quizás los valores… Lee más »
Este programa me ha llevado un montón de esfuerzo y de tiempo, sin paga ni recompensa ni «like» alguno:
t=5;for(r=2,5050,for(s=1,r-1,if(gcd(r,s)==1&((r%2==0&s%2==1)||(r%2==1&s%2==0)),l=2*r*s;n=r^2-s^2;m=r^2+s^2;i=l*m/t;j=n*m/t;if(issquare(i)||issquare(j),print([r,s,l,n,m]);if(floor(i)==ceil(i)&sqrtint(i)==sqrt(i),b=sqrtint(i);a=n;h=(l^2+m^2)/(a*b);c=2*t*b/a,b=sqrtint(j);a=l;h=(n^2+m^2)/(a*b);c=2*t*b/a);print([a/b,c,h,t])))));gettime()
[2, 1, 4, 3, 5]
[3/2, 20/3, 41/6, 5]
[40, 9, 720, 1519, 1681]
[1519/492, 4920/1519, 3344161/747348, 5]
[5041, 242, 2439844, 25353117, 25470245]
[25353117/3525434, 35254340/25353117, 654686219104361/89380740677778, 5]
%3 = 78684 milisegundos
Existe otro algoritmo distinto, pero que aún no he desarrollado para Pari gp, que obtiene todos los triángulos rectángulos de lados racionales, para las áreas enteras de la forma 2(2k + 1), es decir áreas pares pero no múltiplos de 4. No sé si lo desarrollaré porque soy demasiado lento y ya me estoy cansando. Por eso es posible que hasta un max(r,s) dado, s = 34 tenga más soluciones que s = 65, los dos campeones por ahora. Una es dos veces un primo de la forma 4k+1 y el otro el producto de 2 primos de esta forma.… Lee más »
Corrección ligera : Ha de ser ;
….max(r,s) dado, el área t = 34 ( y no s = 34). Porque r y s son los generadores de las ternas pitagóricas.