Problema de Basilea

Problema de Basilea

Problema de Basilea

El Problema de Basilea es un famoso problema de teoría de números, planteado por primera vez por Pietro Mengoli, y resuelto por Leonhard Euler en 1735. Puesto que el problema había resistido los ataques de los matemáticos más importantes de la época, la solución llevó a Euler rápidamente a la fama cuando tenía veintiocho años. Euler generalizó el problema considerablemente, y sus ideas fueron tomadas años después por Bernhard Riemann en su artículo de 1859 On the Number of Primes Less Than a Given Magnitude, en donde definió su función zeta y demostró sus propiedades básicas. El problema debe su nombre a la ciudad de residencia de Euler (Basilea), ciudad donde vivía también la familia Bernoulli, que atacó el problema sin éxito.

El problema de Basilea consiste en encontrar la suma exacta de los inversos de los cuadrados de los enteros positivos, esto es, la suma exacta de la serie infinita:


\sum_{n=1}^\infin \frac{1}{n^2} =
\lim_{n \to \infty}\left(\frac{1}{1^2} + \frac{1}{2^2} + \cdots + \frac{1}{n^2}\right) = \frac{\pi ^2}{6}

Numéricamente, se puede obtener que la serie es aproximadamente igual a 1,644934. Sin embargo, el problema de Basilea busca la suma exacta de la serie, de forma cerrada, así como una demostración de que dicha suma es correcta. Euler encontró que la suma exacta era π2/6 y anunció su descubrimiento en 1735. Sus argumentos estaban basados en manipulaciones que no estaban aun justificadas, y no fue hasta 1741 cuando pudo dar una demostración verdaderamente rigurosa.

Aunque es poco conocido, esta suma puede escribirse, de manera integral, como función de dos variables. A menudo, se coloca como ejercicio para los estudiantes de matemáticas:
\int_0^1\int_0^1 {1 \over\ 1-xy}\,dy\,dx = \frac{\pi ^2}{6}

Contenido

Euler ataca el problema

El método inicial de Euler para la obtención del valor π2/6 es original e ingenioso. En esencia, lo que hizo fue extender resultados aplicables a polinomios finitos, considerándolos también válidos para series infinitas. Claro está que el razonamiento de Euler requiere justificación, pero aun sin ella, simplemente obteniendo el valor correcto, pudo verificarlo numéricamente frente a sumas parciales de la serie. La concordancia observada le dio suficiente confianza como para anunciar su resultado a la comunidad matemática.

Para seguir el razonamiento de Euler, hay que recordar el desarrollo en serie de Taylor de la función seno:

\sin(x) = x - \frac{x^3}{3!} + \frac{x^5}{5!} - \frac{x^7}{7!} + \cdots

Dividiendo por x, tenemos

 \frac{\sin(x)}{x} = 1 - \frac{x^2}{3!} + \frac{x^4}{5!} - \frac{x^6}{7!} + \cdots

Ahora bien, las raíces (ceros) de sin(x)/x se encuentran precisamente en x = n\cdot\pi, donde n = ±1, ±2, ±3... Asumamos que podemos expresar esta serie infinita como producto de factores lineales dados por las raíces, de la misma forma que hacemos con los polinomios finitos:


\frac{\sin(x)}{x} =
\left(1 - \frac{x}{\pi}\right)\left(1 + \frac{x}{\pi}\right)\left(1 - \frac{x}{2\pi}\right)\left(1 + \frac{x}{2\pi}\right)\left(1 - \frac{x}{3\pi}\right)\left(1 + \frac{x}{3\pi}\right) \cdots

= \left(1 - \frac{x^2}{\pi^2}\right)\left(1 - \frac{x^2}{4\pi^2}\right)\left(1 - \frac{x^2}{9\pi^2}\right) \cdots

Si hacemos este producto y agrupamos todos los términos en x2, vemos que el coeficiente de x2 para la función sin(x)/x es


-\left(\frac{1}{\pi^2} + \frac{1}{4\pi^2} + \frac{1}{9\pi^2} + \cdots \right) =
-\frac{1}{\pi^2}\sum_{n=1}^{\infty}\frac{1}{n^2}.

Pero del desarrollo en serie de Taylor original de sin(x)/x, tenemos que el coeficiente de x2 es −1/(3!) = −1/6. Estos dos coeficientes deben ser iguales (por el teorema de unicidad del desarrollo en serie); por tanto,


-\frac{1}{6} =
-\frac{1}{\pi^2}\sum_{n=1}^{\infty}\frac{1}{n^2}.

Multiplicando ambos miembros de la igualdad por −π2 obtenemos la suma de los inversos de los cuadrados de los enteros positivos.

La función zeta de Riemann

La función zeta de Riemann ζ(s) es una de las funciones más importantes en matemáticas, debido a su relación con la distribución de los números primos. La función está definida para todo número complejo s cuya parte real sea mayor que la unidad (Re(s) > 1) por la siguiente fórmula:


\zeta(s) =
\sum_{n=1}^\infin \frac{1}{n^s}.

Haciendo s = 2, vemos que ζ(2) es igual a la suma de los inversos de los cuadrados de los enteros positivos.


\zeta(2) =
\sum_{n=1}^\infin \frac{1}{n^2} =
\frac{1}{1^2} + \frac{1}{2^2} + \frac{1}{3^2} + \frac{1}{4^2} + \cdots \approx 1.644934.

¿Cómo sabemos que la suma converge? Podemos demostrarlo con la siguiente desigualdad:


\zeta(2) =
\sum_{n=1}^\infin \frac{1}{n^2} < 1 + \sum_{n=2}^\infin \frac{1}{(n-1)n} = 1 + \sum_{n=2}^\infin \left( \frac{1}{n-1} - \frac{1}{n} \right) = 2.

Esto marca el límite superior ζ(2) < 2, y como la suma infinita tiene sólo términos positivos, debe converger. Puede demostrarse que ζ(s) se expresa en términos de los números de Bernoulli siempre que s sea un número par positivo.

Una demostración rigurosa

El siguiente razonamiento demuestra la identidad ζ(2) = π2/6, donde ζ(s) es la función zeta de Riemann. Es con mucho la demostración más sencilla hasta ahora disponible; mientras otras demostraciones emplean resultados de matemática avanzada, como análisis de Fourier, análisis complejo, y cálculo multivariable, la siguiente no requiere siquiera cálculo de una variable (aunque es necesaria la obtención de un límite al final).

Historia de la demostración

El origen de la demostración es incierto. Apareció en la revista Eureka en 1982, atribuida a John Scholes, pero Scholes dice haber aprendido la demostración de Peter Swinnerton-Dyer, y en cualquier caso mantiene que ésta era "de dominio público en la Universidad de Cambridge a finales de los 60".

Lo que se necesita saber

Para comprender la demostración, es necesario conocer los siguientes resultados:

(\cos x + i\sin x)^n = \cos (nx) + i \sin (nx). \,\!
Demostración: se puede demostrar a partir de la fórmula de Euler.
(x+y)^n=\sum_{k=0}^n{n \choose k}x^ky^{n-k}
donde los coeficientes binomiales son
{n \choose k}=\frac{n!}{k!(n-k)!} (ver factorial)
Demostración: se necesitan inducción matemática y algunas propiedades de los coeficientes binomiales.
  • La función cot2 x es inyectiva en el intervalo (0, π/2).
    • Demostración: supongamos cot2 x = cot2 y para algún x y algún y en el intervalo (0, π/2). Usando la definición de cotangente cot x = (cos x)/(sin x) y la identidad trigonométrica cos2 x = 1 − sin2 x, vemos que (sin2 x)(1 − sin2 y) = (sin2 y)(1 − sin2 x). Restando (sin2 x)(sin2 y) a cada miembro, tenemos sin2 x = sin2 y. Dado que la función seno es no negativa en el intervalo (0, π/2), esto significa que sin x = sin y, pero es geométricamente evidente (observando el círculo unidad, por ejemplo) que la función seno es inyectiva en el intervalo (0, π/2), por lo que x = y, como queríamos demostrar.
  • Si p(t) = amtm + am − 1tm − 1 +... + a1t + a0, donde am ≠ 0, entonces la suma de las raíces de p (contando la múltiples) es −am − 1/am.
    • Demostración: Si am = 1, entonces p(t) = ∏(ts), donde s abarca todas las raíces de p. Desarrollando este producto, vemos que el coeficiente de tm − 1 es menos la suma de todas las raíces. Si am ≠ 1, podemos dividir cada término por este valor, obteniendo un nuevo polinomio con las mismas raíces, cuyo coeficiente principal es ahora 1; aplicando el razonamiento precedente vemos que la suma de las raíces del polinomio original p(t) = suma de las raíces del nuevo polinomio = −am − 1/am.
  • La identidad trigonométrica csc2 x = 1 + cot2 x.
    • Demostración: se sigue de la identidad fundamental 1 = sin2 x + cos2 x tras dividir ambos miembros por sin2 x.
  • Para todo número real x siendo 0 < x < π/2, se cumplen las desigualdades cot2 x < 1/x2 < csc2 x.
    • Demostración: fijémonos primero en que 0 < sin x < x < tan x. Esto puede verse en el siguiente esquema:

Circle-trig6.svg

Para ver que 0 < sin x < x, hay que observar que en el dibujo, sin θ es la longitud del segmento AC, y θ la longitud del arco AD.
Para ver que x < tan x, obsérvese que el área del triángulo OAE es tan(θ)/2, el área del sector OAD es θ/2, y que el sector está contenido en el triángulo.
Ahora, tomemos los inversos y elevemos al cuadrado en las desigualdades, recordando que éstas deben cambiar el sentido. Obtenemos así las desigualdades buscadas.


Otra demostración, menos geométrica, puede conseguirse con: sea f(x) = sin x - x, vemos que f(0) = 0 y f'(x) = cos x - 1 < 0 (en casi todo punto), lo que implica que f es estrictamente decreciente, deduciendo que f(x) < 0 y, por tanto, que sin x < x.
Para ver la otra desigualdad, definimos g(x) = tan x - x y aplicamos el mismo razonamiento.


  • Sean a, b, y c números reales, con a y c no nulos; entonces el límite de la función (am + b)/(am + c) cuando m tiende a infinito es 1.
    • Demostración: dividamos cada término por m, obteniendo (a + b/m)/(a + c/m). Si dividimos un número fijo por un número cada vez mayor, el cociente se aproxima a cero; por tanto, tanto el numerador como el denominador tienden a a, por lo que su cociente tenderá a 1.
  • El criterio del sandwich, que establece que si una función se encuentra acotada entre otras dos funciones, y estas dos funciones tienden ambas a un mismo límite, entonces la función acotada por ellas tiende igualmente a dicho límite.
    • Demostración: ver el artículo squeeze theorem (en inglés) para una discusión y demostración detalladas.

La demostración

La idea principal en la que descansa la demostración es la de limitar las sumas parciales

\sum_{k=1}^m \frac{1}{k^2} = \frac{1}{1^2} + \frac{1}{2^2} + \cdots + \frac{1}{m^2}

entre dos expresiones, las cuales tenderán ambas a π2/6 cuando m tienda a infinito. Las dos expresiones se obtienen de identidades que incluyen las funciones cotangente y cosecante. Estas identidades son a su vez obtenidas de la fórmula de De Moivre, y eso es lo que nos disponemos a hacer ahora.

Sea x un número real tal que 0 < x < π/2, y sea n un entero positivo. Entonces de la fórmula de De Moivre y de la definición de cotangente, tenemos

\frac{\cos (nx) + i \sin (nx)}{(\sin x)^n} = \frac{(\cos x + i\sin x)^n}{(\sin x)^n} = \left(\frac{\cos x + i \sin x}{\sin x}\right)^n = (\cot x + i)^n.

Del binomio de Newton, tenemos

(\cot x + i)^n = {n \choose 0} \cot^n x + {n \choose 1} (\cot^{n-1} x)i + \cdots + {n \choose {n-1}} (\cot x)i^{n-1} + {n \choose n} i^n
= \left[ {n \choose 0} \cot^n x - {n \choose 2} \cot^{n-2} x \pm \cdots \right] \; + \; i\left[ {n \choose 1} \cot^{n-1} x - {n \choose 3} \cot^{n-3} x \mp \cdots \right].

Combinando ambas ecuaciones e igualando las partes imaginarias, obtenemos la identidad

\frac{\sin (nx)}{(\sin x)^n} = \left[ {n \choose 1} \cot^{n-1} x - {n \choose 3} \cot^{n-3} x \mp \cdots \right].

Tomamos esta identidad y establecemos que n = 2m + 1, siendo m un entero positivo, y x = r π/(2m + 1), siendo r = 1, 2,..., m. Entonces nx = r π, por lo que sin(nx) = 0 y, de esta forma,

0 = {{2m+1} \choose 1} \cot^{2m} x - {{2m+1} \choose 3} \cot^{2m-2} x \mp \cdots + (-1)^m{{2m+1} \choose {2m+1}}.

Esta ecuación es válida para todos los valores de x = r π/(2m + 1), tales que r = 1, 2,..., m. Estos valores de x son números distintos estrictamente entre 0 y π/2. Dado que la función cot2(x) es inyectiva en el intervalo (0, π/2), los números cot2(x) = cot2(r π/(2m + 1)) son distintos para cada valor de r = 1, 2,..., m. Pero por la ecuación anterior, cada uno de estos m números distintos es una raíz del polinomio de grado m

p(t) := {{2m+1} \choose 1}t^m - {{2m+1} \choose 3}t^{m-1} \mp + (-1)^m{{2m+1} \choose {2m+1}}.

Esto significa que los números x = cot2(r π/(2m + 1)), para r = 1, 2,..., m son precisamente las raíces del polinomio p(t). Pero podemos calcular directamente la suma de las raíces examinando los coeficientes, y la comparación muestra que

\cot ^2 \left(\frac{\pi}{2m+1}\right) + \cot ^2 \left(\frac{2 \pi}{2m+1}\right) + \cdots + \cot ^2 \left(\frac{m \pi}{2m+1}\right)
= {{2m+1} \choose 3} / {{2m+1} \choose 1} = \frac{(2m)(2m-1)}{6}.

Sustituyendo la identidad csc2 x = cot2 x + 1, tenemos

\csc ^2 \left(\frac{\pi}{2m+1}\right) + \csc ^2 \left(\frac{2 \pi}{2m+1}\right) + \cdots + \csc ^2 \left(\frac{m \pi}{2m+1}\right)
= {{2m+1} \choose 3} / {{2m+1} \choose 1} + m = \frac{(2m)(2m+2)}{6}.

Consideremos ahora la desigualdad cot2 x < 1/x2 < csc2 x. Si sumamos todas estas desigualdades para cada uno de los valores x = r π/(2m + 1), y empleamos las dos identidades anteriores, obtenemos

\frac{(2m)(2m-1)}{6} < \left( \frac{2m+1}{\pi} \right) ^2 + \left( \frac{2m+1}{2 \pi} \right) ^2 + \cdots + \left( \frac{2m+1}{m \pi} \right) ^2 < \frac{(2m)(2m+2)}{6}.

Multiplicando todo por (π/(2m + 1))2, pasamos a

\frac{\pi ^2}{6}\left(\frac{2m}{2m+1}\right)\left(\frac{2m-1}{2m+1}\right) < \frac{1}{1^2} + \frac{1}{2^2} + \cdots + \frac{1}{m^2} < \frac{\pi ^2}{6}\left(\frac{2m}{2m+1}\right)\left(\frac{2m+2}{2m+1}\right).

Cuando m tiende a infinito, las expresiones a la izquierda y a la derecha tienden ambas a π2/6, luego por el criterio del sandwich,

\zeta(2) =
\sum_{k=1}^\infin \frac{1}{k^2} =
\lim_{m \to \infty}\left(\frac{1}{1^2} + \frac{1}{2^2} + \cdots + \frac{1}{m^2}\right) = \frac{\pi ^2}{6}

con lo que queda completa la demostración. C.Q.D.

Aplicación al cálculo de π

En 1881, Ernesto Cesaro demostró que la probabilidad de que dos enteros sean primos relativos es 6 / π2, que es el inverso de ζ(2). Mediante la demostración de más arriba, el teorema de Cesaro permite calcular un valor para π a partir de una muestra grande de enteros aleatorios, determinando la proporción de éstos que son primos relativos.

Referencias

  • Number Theory: An Approach Through History, Andre Weil, Springer, ISBN 0-8176-3141-0
  • Euler: The Master of Us All, William Dunham, MAA, ISBN 0-88385-328-0
  • Prime Obsession: Bernhard Riemann and the Greatest Unsolved Problem in Mathematics, John Derbyshire, Joseph Henry Press, ISBN 0-309-08549-7
  • Proofs From the Book, Martin Aigner, Gunter Ziegler, Springer, ISBN 3-540-67865-4
  • Riemann's Zeta Function, Harold M. Edwards, Dover, ISBN 0-486-41740-9

Enlaces externos

Obtenido de "Problema de Basilea"

Wikimedia foundation. 2010.

Mira otros diccionarios:

  • Basilea (desambiguación) — Saltar a navegación, búsqueda El término Basilea puede referirse a: Basilea ciudad suiza, capital del semi cantón de Basilea Campaña Basilea II una normativa internacional para el control del riesgo bancario Basilea cantón suizo, dividido en:… …   Wikipedia Español

  • Basilea, Concilio de — ► RELIGIÓN Concilio convocado en 1431 por el papa Martín V y trasladado por Eugenio IV a Ferrara (1437) y Florencia. Fue disuelto en 1442. Se propuso afrontar la problemática presentada por los husitas. * * * (1431–49). Concilio de la Iglesia… …   Enciclopedia Universal

  • Museos de Basilea — Indicador de museos en Basilea. Por museos de Basilea se entiende un grupo de museos de la ciudad de Basilea y su área de influencia. Cubren un amplio espectro de exposiciones, con el punto fuerte puesto en las artes plásticas. Albergan una gran… …   Wikipedia Español

  • Leonhard Euler — Retrato de Leonhard Euler, pintado por Johann Georg Bruck …   Wikipedia Español

  • Número π — π (pi) es la relación entre la longitud de una circunferencia y su diámetro, en geometría euclidiana. Es un número irracional y una de las constantes matemáticas más importantes. Se emplea frecuentemente en matemáticas, física e ingeniería. El… …   Wikipedia Español

  • 1 − 2 + 3 − 4 + · · · — Los primeros miles de términos y sumas parciales de 1 − 2 + 3 − 4 + · · ·. En matemáticas, la expresión 1 − 2 + 3 − 4 + · · · es una serie infin …   Wikipedia Español

  • Número primo — Un número primo es un número natural mayor que 1, que tiene únicamente dos divisores distintos: él mismo y el 1. Se contraponen así a los números compuestos, que son aquellos que tienen algún divisor natural aparte de sí mismos y del 1. El número …   Wikipedia Español

  • Constante zeta — En matemática, una constante zeta es el resultado de la función zeta de Riemann cuando esta se aplica sobre un número entero. Contenido 1 La función zeta de Riemann en 0 y en 1 2 Enteros positivos 2.1 Enteros positivos pares …   Wikipedia Español

  • Serie convergente — En matemáticas, una serie se llama serie convergente si la sucesión de sumas parciales tiene un límite en el espacio considerado. De otro modo se llama serie divergente. Contenido 1 Definición formal 2 Ejemplos 3 Convergencia absoluta …   Wikipedia Español

  • Carl Gustav Jung — «Jung» redirige aquí. Para otras acepciones, véase Jung (desambiguación). Carl Gustav Jung …   Wikipedia Español

Compartir el artículo y extractos

Link directo
Do a right-click on the link above
and select “Copy Link”