Calcular el área bajo la campana de Gauss
En la cita Ouch! del pasado mes de mayo dani comentó aquí que la demostración de que el área entre la curva descrita por la campana de Gauss y el eje X es 1 era hermosa. En este artículo os dejo la demostración que conozco de este hecho. Si conocéis otra no dudéis en comentarla.
Introducción

La función definida de la siguiente forma:
se denomina función gaussiana y su gráfica tiene forma de campana. Tomando ciertos valores de y
obtenemos que esta función es la función de densidad de una variable aleatoria normal:
Si
es una variable aleatoria que sigue una distribución normal de media
y desviación típica
(suele escribirse también de la forma
), entonces
tiene como función de densidad a:
Para comprobar que es una función de densidad debemos comprobar estas dos condiciones:
1.-
2.-
La primera condición es evidente, al ser (por definición de desviación típica) y por ser la exponencial siempre positiva. La comprobación de la segunda condición consiste simplemente en el cálculo de esa integral impropia…¿Simplemente?
Quienes no conozcan esta integral pueden intentar calcularla encontrando una primitiva y evaluando después en y en
. Tiempo perdido. No podemos encontrar dicha primitiva. Pero, por suerte, sí podemos calcular esa integral, aunque no con cálculo integral de una variable sino con cálculo integral de dos variables.
Explicación previa
Para empezar es evidente que es simétrica respecto de la recta
(sustituid
por
y
y veréis que se obtiene el mismo resultado). Por simplicidad para los cálculos tomaremos
y
, es decir, la función de densidad de la conocida
. Esto es, calcularemos la integral siguiente:
Con ello, por tanto, obtenemos una función par (simétrica respecto del eje Y), por lo que bastará con calcular la integral entre e
y después multiplicar por dos el resultado.
Hemos dicho que vamos a utilizar cálculo integral de dos variables. Concretamente calcularemos una integral doble (en dos variables) y deduciremos de ese cálculo cuánto vale la integral buscada. Para el cálculo de dicha integral doble serán esenciales las coordenadas polares.
Cálculo de la integral
La integral doble que vamos a calcular es la siguiente:
Para ello realizamos un cambio a coordenadas polares:
Al realizar este cambio los intervalos en los que toman valores y
son los siguientes:
Por otra parte se tiene que en polares , por lo que el exponente de la exponencial quedará
. Y por realizar un cambio de variable en una integral doble debemos multiplicar por el valor absoluto del jacobiano del cambio de variable, que en este caso es
.
Ya estamos preparados para plantear la integral:
Al quedarnos una función independiente de obtenemos:
Calculamos dicha integral:
Obtenemos entonces la primera parte de nuestro cálculo:
Por otro lado:
Sacamos la función independiente de (la segunda exponencial) de la integral respecto de
:
Como la integral que queda dentro es independiente de puede sacarse completamente de la integral respecto de
, quedando:
Obtenemos entonces un producto de dos integrales iguales simplemente renombrando la de la segunda integral a
(esto no supone ningún problema, es simplemente un cambio de nombre de la variable). Tenemos entonces la integral al cuadrado. Tomando el principio y el final de este paso llegamos a:
Uniendo los resultados obtenidos de estos dos pasos llegamos a:
Tomando raíces cuadradas se tiene lo siguiente:
Recordemos que la función a integrar era par, por lo que el valor de la integral completa será el resultado de multiplicar lo obtenido por dos, es decir:
Ahora tomamos la definición inicial de la función y llegamos al resultado buscado:
Nota:
La imagen que ilustra este post fue la propuesta de logo de nuestro lector Piponauta para el Concurso: Logo para Gaussianos que organicé el pasado año 2008.






Trackback | 11 Jun, 2009
Bitacoras.com
Cristobal | 11 de Junio de 2009 | 9:36
Publiqué sobre la resolución de la integral gaussiana por el método que explicas y me sorprendió saber que mucha gente no lo conocía el método.
Dichosas oposiciones
Propongo un ejercicio que me propusieron en los comentarios, resolverla cuando la x esta elevada al cubo. Aún no he tenido tiempo de calcularla
trollll | 11 de Junio de 2009 | 11:07
Lo único que no entiendo es que, el ángulo theta del cambio a polares, se mueva entre 0 y Pi/2… No sé por qué se mueve sólo hasta Pi/2
Entendería que se moviera hasta Pi, por estar calculando la mitad del plano, pero no entiendo el que se mueva sólo hasta Pi/2
Un saludo y gracias!
Dani | 11 de Junio de 2009 | 11:25
jejeje sí, esa es la demostración de la que hablaba. Es simplemente genial meterte en el cálculo de dos variables para resolver un problema de una variable. Como no, la prueba es de Gauss. Me la hicieron en clase este año (Cálculo II) y me quedé con una sonrisa de oreja o oreja jajaja. Cristobal, así a bote pronto se me ocurre reciclar la prueba de la campana pero para 3 variables:
Si ahora acemos el cambio a coordenadas esféricas tenemos 
Dani | 11 de Junio de 2009 | 11:33
entonces
Algo así será. (^DiAmOnD^, qué comando usas para la evaluación de las primitivas?? yo uso \big \bracevert (
) pero no queda muy bien la verdad :S)
Dani | 11 de Junio de 2009 | 11:37
Con lo cual te quedaría
Por supuesto luego cambiando constantes, la fracción del exponente de e blah blah blah te quedaría distinto, y por supuesto unitario si lo ajustas.
Dani | 11 de Junio de 2009 | 11:47
trollll, te digo: la integral en cartesianas tiene limites de integracion
Eso quiere decir que estamos integrando en la región del plano
, es decir en el primer cuadrante sólo. Por lo tanto para expresar todo el primer cuadrante en coordenadas polares necesitamos (obviamente que
) y que el ángulo polar desde el eje x positivo esté entre
y
.
Dani | 11 de Junio de 2009 | 12:01
arghhhh! no! está mal mi cálculo. se me olvidó dividir por todos los senos y cosenos que multiplican a r^3 en el exponente. buf… eso va a ser más complicado :S…
trollll | 11 de Junio de 2009 | 12:14
Correcto Dani, no me di cuenta de lo de los límites de integración, creí que sólo integrábamos la primera mitad del plano (entonces los límites serían 0<x< inf ; – inf<y< inf)
Todo arreglado!
M | 11 de Junio de 2009 | 13:05
Cuidado: la integral
es divergente (integrado no acotado). Por otro lado
.
Dani | 11 de Junio de 2009 | 13:28
(se me olvidó el “-” en el exponente
)
^DiAmOnD^ | 11 de Junio de 2009 | 14:43
Dani, he usado esto:
\Bigg ]_0^{\infty}
De toas formas, si pinchas en una imagen generada con
con el botón derecho del ratón y miras las propiedades puedes ver el código utilizado.
M | 11 de Junio de 2009 | 16:37
Tras desempolvar apuntes de la carrera, comento una prueba que sólo hace uso de cálculo integral de una variable. Aunque evita el cálculo de integrales dobles, usa la regla de Leibniz para derivar integrales paramétricas.
Sea
, (
>0). Derivando la integral,
, y haciendo el cambio
, vemos que
Así que
. Ya que
, vemos que
. Finalmente, ya que
, obtenemos el valor
.
hernan | 11 de Junio de 2009 | 17:42
Linda la última, no la conocía.
Pero por lo que veo no requiere saber la regla de Leibnitz completa, porque la variable de derivación está sólo en el límite de integración y no en el integrando (con lo cual, con el teorema fundamentental del cálculo más la regla de la cadena ya está).
El caso de la derivada “full”, en mi experiencia, suele dejar perplejos a mis estudiantes (ingeniería, últimos años de la carrera) que no lo han visto (y tras cursar Análisis I creen que saben derivar cualquier cosa), y que no saben qué hacer cuando la variable aparece en el límite de integración Y en el integrando. Por supuesto, “haciendo derivada por definición” sale.
vengoroso | 12 de Junio de 2009 | 12:30
Huy, me acuerdo de esta integral… nos la pusieron como ejercicio en estadística de primero :-/ Recuerdo que un compañero la resolvió usando las funciones Gamma y Beta de Euler, pero hace muchos años de aquello y no recuerdo los detalles.
Otra opción es pasarse a variable compleja y usar el teorema de los residuos. El haber perdido tanto tiempo con estos problemas es lo que te hace apreciar de verdad la variable compleja
Gaby | 12 de Junio de 2009 | 16:00
Pienso que mis compañeros de la carrera (Matemáticas) deberían visitar páginas como Gaussianos más a menudo, parece que en mi clase estas inquietudes brillan por su ausencia.
Américo Tavares | 13 de Junio de 2009 | 0:46
Sou um admirador vosso.
Um integral muito mais simples que se pode calcular usando coordenadas polares é:

a partir do qual podemos chegar ao resultado clássico
Angus Taylor, Advanced Calculus, p.663
Jose | 13 de Junio de 2009 | 1:11
” Sea f:[a,b]–> R tal que todo subintervalo cerrado cumple el Teorema del valor intermedio(no necesariamente continua). ¿Se cumple que la funcion es continua?, caso contrario dar un contraejemplo.”
Agradecería mucho que me echaran una mano en este asunto. Saludos
Dani | 13 de Junio de 2009 | 2:02
define
Se comprueba f\’acilmente que definida por ejemplo en el
cumple el teorema del valor intermedio en todo subintervalo cerrado, sin embargo no es continua en
.
Américo Tavares | 13 de Junio de 2009 | 10:19
Em Problema do mês #1 coloquei o seguinte enunciado:
Um leitor já demonstrou que 5 é um majorante. Não sei provar sem computador que
Alguém saberá?
fede | 13 de Junio de 2009 | 11:48
Américo, creo que el primer comentario de PB a tu post ya demuestra que en la descomposición de
en factores primos, cada componente
, y en consecuencia
.
Américo Tavares | 13 de Junio de 2009 | 12:14
fede,
, com
.
De acordo!
Enganei-me, visto que foi demonstrado o máximo possível. A questão que falta demonstrar (sem o computador) é que
Espero agora não ter errado!
Obrigado.
Américo Tavares | 13 de Junio de 2009 | 12:21
… foi demonstrado que 5 é um majorante …
Jose | 13 de Junio de 2009 | 16:35
Hola Dani, gracias por responder. Aunque no me queda del todo claro. La funcion es sen1/x cuando x es distinto de 0. ¿Y en 0?. ¿Y como pruebo que en una subcaja que contiene a 0, la funciòn efectivamente cumple TVI?.
Dani | 13 de Junio de 2009 | 17:10
da igual como definas
di por ejemplo
. Probemos que todo subintervalo cerrado del
comprueba el TVI. sea
Si el intervalo no contiene al
la función es continua en ese intervalo y entonces el TVI se cumple trivialmente. Por lo tanto podemos suponer
. Tenemos que probar que la función toma todos los valores entre
y
. Si
sea
suficientemente grande como para que
Entonces también tenemos
pero entre
y
toma todos los valores entre
y
, pues el seno tiene periodo
. En particular es seguro que toma todos los valores entre
y
. Si
y
el razonamiento es idéntico pero por “la izquierda del cero”. Esto demuestra que
(con digamos
) cumple el TVI en todo subintervalo cerrado del
, y sin embargo no es continua en
. Espero que esto resuelva tus dudas.
Dani | 13 de Junio de 2009 | 19:29
quería decir:
y
el razonamiento es idéntico pero por…
*si
Américo Tavares | 14 de Junio de 2009 | 10:25
fede,
Para ficar mais claro, repito, ligeiramente editado:
«De acordo!
Enganei-me, visto que foi demonstrado que 5 é um majorante.
A questão que falta demonstrar (sem o computador) é que para
,
.»
Desculpem a repetição.
Cumprimentos
fede | 14 de Junio de 2009 | 11:56
Américo, si
es el exponente de p en la descomposición de n en factores primos y k no divide a
,
.
Porque, usando la fórmula de Legendre ,
Por tanto
divide a
, si k no divide a
.
Américo Tavares | 14 de Junio de 2009 | 11:59
fede
Obrigado pela demonstração.
Américo
Américo Tavares | 14 de Junio de 2009 | 12:08
Posso publicar no meu blog esta demonstração, indicando a fonte?
Américo
fede | 14 de Junio de 2009 | 12:14
Sí, por supuesto.
fede | 14 de Junio de 2009 | 12:32
Me doy cuenta ahora que en mi comentario anterior la condición para que todos los sumandos sean 1 es errónea.
La condición correcta es que p no divida a k.
Jose | 16 de Junio de 2009 | 0:53
Gracias Dani ya lo contemplo. Aunque me gustaría saber como fue que llegaste a esa función o si ya lo habías visto antes
Américo Tavares | 16 de Junio de 2009 | 8:29
fede,
Com a sua correcção fica:
«Si
es el exponente de
en la descomposición de
en factores primos y
no divide a
,
.
Porque, usando la fórmula de Legendre,

Por tanto
divide a
, si
no divide a
.»
Obrigado novamente!
Dani | 16 de Junio de 2009 | 10:30
pues no sé, la función es un clásico ejemplo de discontinuidad patológica, y me pareció que cumplía los requisitos del enunciado… es como si lees algún problema de función discontinua en todos o casi todos los puntos: en seguida se te pasa por la cabeza alguna función que tenga que ver con la de dirichlet! de nada!
fede | 16 de Junio de 2009 | 13:01
Américo, debe decir “… y p no divide a k,” y no “… y k no divide a p,”.
Y la última frase debe ser “si 13 no divide a k”.
En general se tiene que
, porque la expresión
vale 0 cuando
divide a k, y 1 cuando
no divide a k. Entonces
Américo Tavares | 17 de Junio de 2009 | 11:58
Informo ter hoje publicado no meu blog(ue) um link para
13 de Junio de 2009 | 10:19
«Em 13.06.09 divulguei este problema nos Gaussianos.
On July 13, 2009 I posted this Problem in Gaussianos.»
Trackback | 19 Jun, 2009
problemas | teoremas
Trackback | 22 Jun, 2009
Usando a fórmula de Legendre « problemas | teoremas
Trackback | 20 Jul, 2009
Problema do mês, Problem of the month. Resolução do Problema #1, Solution of Problem #1 « problemas | teoremas
Cristobal | 30 de Agosto de 2009 | 10:57
@Dani-> Con x^3+y^3+z^3 no es viable el cambio a esféricas, además en tu desarrollo falta la derivada del exponente de la exponencial.
@M-> Si no me he equivocado al ver tu razonamiento con el método de Leibniz a mi el término:
f(b(t),t)*b’(t) = exp(-2t^2)*1
Y a ti te sale que es cero.
Sigo con ello
Cristobal | 30 de Agosto de 2009 | 11:01
@M-> Ah ya, haces tender t a infinito y da cero, lo que pasa es que haces el límite un poco pronto
juancho | 17 de Octubre de 2009 | 23:36
Saludos a todos… entonces si tengo la integral desde – infinito hasta + infinito de exp(-ax^2)dx… como pruebo que el resultado es (pi/a)^(1/2), me ayudan por favor????
juancho | 17 de Octubre de 2009 | 23:37
Es urgente!!!!!!!!!!!! De antemano gracias
Dani | 18 de Octubre de 2009 | 0:20
simplemente copia la demostración que está puesta arriba cambiando las constantes…
Trackback | 2 Nov, 2009
Funciones sin primitiva elemental | Gaussianos
alexander | 5 de Noviembre de 2009 | 3:00
por favor necesito ayuda que alguien explique un ejemplo de como aplico la formula en la estadistica de un numero con la campanade gauss ej tengo 100 numeros de dos cifras y en los proximos 100 sorteos futuros saldra alguno o lagunos de los dos numeros de dos cifras como se hace, le agradeceria en el alma gracias