En este capítulo extendemos el concepto de integral de Riemann visto en el capítulo de integrales a funciones de varias variables, en particular en y donde las integrales múltiples tienen multitud de aplicaciones en Ciencias e Ingenierías.
También veremos que el Teorema Fundamental del Cálculo sigue siendo válido y que el cálculo de integrales con respecto a una variable es el proceso inverso del cálculo de la derivada parcial con respecto a esa misma variable.
Sumas de Riemann
Las sumas de Riemann que se introdujeron en la Sección 8.1 pueden generalizarse fácilmente a funciones de varias variables, pero primero debemos generalizar el concepto de intervalo cerrado en .
Definición 13.1 (Intervalo dimensional) Dados intervalos cerrados de , su producto cartesiano se conoce como intervalo dimensional cerrado o hiperrectángulo.
En , el intervalo define un rectángulo, mientras que en el intervalo define una caja.
Definición 13.2 (Partición de un intervalo dimensional) Dado un intervalo dimensional cerrado , con , y una partición de cada intervalo , el producto cartesiano se conoce como partición dimensional de .
El conjunto de todas las particiones dimensionales de se denota .
Definición 13.3 (Sumas de Riemann) Dada una función acotada en el intervalo dimensional y una partición dimensional de , se define la suma inferior de Riemann de respecto de , y se denota , como
donde .
Del mismo modo, se define la suma superior de Riemann de respecto de , y se denota , como
donde .
Aunque la definición es un poco enrevesada para , si es una función de dos variables acotada en el intervalo , y es una partición de con y , las sumas de Riemann inferior y superior de respecto de son
donde y son el ínfimo y el supremo de .
Gráficamente, si es una función positiva en , la suma inferior se puede interpretar como la suma de los volúmenes de los prismas de base rectangular con lados y y altura .
Ejemplo 13.1 Veamos cómo calcular las sumas de Riemann de la función en el rectángulo para la partición
Esta partición descompone el intervalo bidimensional $ en rectángulos (en realidad cuadrados) como los que se ven la siguiente figura.
Como es creciente tanto en como en en el intervalo , alcanzará el mínimo en el extremo inferior izquierdo de cada subintervalo y el máximo en el extremo superior derecho. Por tanto, la suma inferior de Riemann es
Y la suma superior de Riemann es
Integral de Riemann múltiple
Cuando la función es positiva en el intervalo , la suma inferior de Riemann nos da una aproximación por defecto del volumen encerrado entre la superficie de la función y el plano en la región definida por , mientras que la suma superior de Riemann nos da una aproximación por exceso. Al hacer refinamientos de la partición , cada vez con mayor número de subintervalos, las sumas inferiores crecen y las superiores decrecen, dando aproximaciones cada vez mejores, de manera que en el límite, cuando el número de subintervalos tiende a , obtendremos el volumen real.
Definición 13.4 (Integral inferior y superior de Riemann) Dada una función acotada en el intervalo dimensional y una partición dimensional de , se define la integral inferior de Riemann de en , como
Y se define la integral superior de Riemann de en , como
Definición 13.5 Dada una función acotada en el intervalo dimensional , se dice que es integrable Riemann en si
y a este valor se le llama integral de Riemann o integral definida de en y se denota por
donde .
Proposición 13.1 Si es integrable Riemann en el intervalo dimensional , entonces
donde es una partición de en subintervalos de igual tamaño y es el hipervolumen de los subintervalos de la partición.
En el caso de una función de dos variables la proposición anterior establece que si es integrable Riemann en un intervalo , entonces
donde y , es el área de los rectángulos definidos por cada subintervalo de la partición. Es decir, la integral de Riemann aparece al tomar particiones cada vez más refinadas con subintervalos de igual tamaño y hacer la suma de las áreas de los subintervalos multiplicadas por el valor de la función en el extremo superior derecho del subintervalo.
Cuando se quieren hacer explícitas las variables de las que depende la función , la integral de Riemann se suele escribir
para funciones de dos variables y
para funciones de tres variables.
Propiedades de las Integrales de Riemann múltiples
Las propiedades de la integral de Riemann para funciones de una variable también se pueden generalizar para integrales de Riemann múltiples.
Teorema 13.1 Si son dos funciones integrables Riemann en un intervalo dimensional , entonces
es integrable Riemann en y
Para cualquier , es integrable Riemann en y .
Si , entonces .
Prueba. La demostración es similar a la vista para funciones de una variable y se deja como ejercicio.
Integrales múltiples iteradas
El cálculo de integrales de Riemann tomando límites de sumas de Riemann resulta poco operativo en la práctica, pero afortunadamente, podemos utilizar el teorema fundamental del cálculo para calcular integrales de manera mucho más rápida y sencilla. En esta sección veremos como calcular integrales dobles y triples como integrales iteradas, que se pueden calcular realizando varias integrales simples.
El procedimiento para reducir el cálculo de una integral múltiple al cálculo de varias integrales simples es, en el fondo, el mismo que se utilizó para calcular derivadas parciales, es decir, realizar integrales simples de las funciones parciales que se obtienen al fijar como constantes todas las variables excepto aquella con respecto a la que se integra.
Veamos el procedimiento para una función de dos variables. Supongamos que es integrable Riemann en el intervalo . Si mantenemos constante, podemos calcular la integral definida simple con respecto a en el intervalo , es decir, . Para cada valor de , esta integral tomará un valor distinto, por lo que podemos definir la función
que solo depende de . Si ahora se integra esta función con respecto a en el intervalo , se tiene
Esta integral se conoce como integral doble iterada, primero con respecto a y después con respecto a .
De manera similar, podríamos haber integrado primero con respecto a y luego con respecto a , y tendríamos esta otra integral doble iterada
El siguiente teorema establece que una integral múltiple se puede calcular mediante integrales iteradas.
Teorema 13.2 (Fubini) Si es una función continua en un intervalo , entonces
Prueba. Daremos una demostración intuitiva. Si es positiva, la integral doble es el volumen del sólido encerrado entre la superficie de y en plano en el intervalo . Si fijamos el valor de a un valor constante e intersecamos este sólido con el plano , se obtiene una región plana cuya área se puede calcular mediante la integral
De esta manera, para cada la función
nos dará el área de la región plana que se obtiene de la intersección del sólido con el plano que pasa por y es perpendicular al eje . Por tanto, si integramos estas áreas con respecto a en el intervalo ,
obtendremos el volumen del sólido, tal y como se vió en la Sección 8.8.
El mismo razonamiento se puede usar cambiando el orden de integración.
Ejemplo 13.2 La integral doble de la función en el intervalo es
Podríamos haber llegado a este mismo resultado haciendo la otra integral iterada
Aunque cuando la función es integrable Riemann el orden de integración de las integrales iteradas no importa, a menudo, una integral iterada suele ser más fácil de calcular que la otra.
Ejemplo 13.3 Para calcular la integral de la función en el intervalo , si primero integramos con respecto a y después con respecto a se tiene
Sin embargo, si primero integramos con respecto a y luego con respecto a se tiene una integral mucho más difícil de calcular.
Calculadora de sumas de Riemann
En el caso de que pueda factorizarse como producto de funciones que sólo dependen de o de , entonces el cálculo de la integral doble se simplifica.
Proposición 13.2 Si es integrable Riemann en el intervalo , entonces
Prueba. (1) no depende de .
(2) no depende de .
Ejemplo 13.4 La integral doble de la función en el intervalo es
Integrales sobre regiones no regulares
En muchos problemas reales, como por ejemplo en el cálculo de volúmenes, la región de integración no es un integración no es un intervalo dimensional o hiperrectángulo, sino una región no regular, es decir, con fronteras curvas. En esta sección veremos cómo calcular integrales dobles de funciones de dos variables para algunos tipos de regiones irregulares acotadas del plano .
En general, la integral de Riemann de una función sobre un intervalo dimensional puede extenderse a una región irregular acotada simplemente tomando un intervalo dimensional que incluya la región, es decir y redefiniendo como
.
y calculando la integral de Riemann de sobre como se ha visto en la sección anterior, es decir,
Este procedimiento tiene sentido siempre que sea continua sobre la frontera de la región , ya que el valor de en los puntos que están fuera de la región es , y por tanto, no contribuyen la suma de la integral.
En el caso de funciones de dos variables, existen varios tipos de regiones en el plano donde la integral de Riemann se puede calcular mediante integrales iteradas donde los límites de integración dependen de funciones que delimitan la región de integración.
Si es una región del plano delimitada por las gráficas de dos funciones y , es decir,
la integral de Riemann sobre esta región se puede calcular mediante la integral iterada
Del mismo modo, si es una región del plano delimitada por las gráficas de dos funciones y , es decir,
la integral de Riemann sobre esta región se puede calcular mediante la integral iterada
Ejemplo 13.5 Veamos cómo calcular la integral doble de la función sobre la región encerrada entre las curvas e . La región encerrada entre las curvas se muestra en el siguiente gráfico.
Esta región puede expresarse como , por lo que la integral doble de sobre esta región puede calcularse mediante la integral iterada
En este caso, la región también puede expresarse como , por lo que la integral doble de sobre esta región también puede calcularse mediante la integral iterada
En este tipo de integrales iteradas, hay que integrar siempre primero con respecto a la variable que está acotada por las funciones que definen la región y después con respecto a la variable que está acotada por valores fijos. Cuando interese realizar la integral iterada en el otro orden, porque resulte más sencilla, previamente hay que cambiar la expresión de la región de integración mediante las funciones inversas, siempre y cuando sea posible.
Ejemplo 13.6 Para calcular la integral de la función en la región , como está acotada por una función de , habría que hacer la integral iterada
Sin embargo, la función no tiene primitiva elemental, por lo que necesariamente hay que realizar la integral iterada en el otro orden, pero para ello previamente hay que expresar la región como , con lo que la integral doble sobre esta región es
Integrales dobles en coordenadas polares
Al que igual que se vio en la Sección 8.6.3, para integrar determinadas funciones es preferible hacer un cambio a coordenadas polares. Recordemos que para pasar de coordenadas cartesianas a coordenadas polares se utilizan las fórmulas
y para pasar de coordenadas polares a cartesianas se utiliza la fórmula
En coordenadas polares un intervalo bidimensional define un sector de disco como el que se muestra a continuación.
Si la función es continua sobre un intervalo polar con , la integral doble de sobre esta región se puede calcular mediante la integral iterada
Ejemplo 13.7 La gráfica de la función es una semiesfera de radio uno centrada en el origen de coordenadas. Para calcular su volumen, la región de integración es el círculo de radio uno centrado en el origen de coordenadas, que puede expresarse como . En este caso la integral en coordenadas rectangulares es complicada y merece la pena hacer la integral en coordenadas polares.
Y por tanto, el volumen de la esfera de radio uno centrada en el origen es , que coincide con la fórmula habitual.
Cuando la región de integración es irregular pero puede expresarse de la forma , la integral doble en coordenadas polares puede calcularse mediante la integral iterada
Ejemplo 13.8 Veamos cómo calcular el volumen encerrado entre la gráfica integral del paraboloide y el plano sobre la circunferencia de ecuación . En coordenadas polares, la circunferencia que delimita la región de integración puede escribirse como . Por tanto, podemos calcular la integral de sobre la región mediante la integral iterada
Cálculo de volúmenes
Cuando es positiva en el intervalo de integración , la integral doble da el volumen del sólido encerrado entre la superficie de y el plano en la región delimitada por . Ahora bien, si es negativa en , entonces el valor de la integral doble es negativo, por lo que habrá que cambiarla de signo para poder interpretarla como un volumen.
La cosa se complica cuando puede tomar valores positivos y negativos en el intervalo . En tal caso, para obtener el volumen tendremos que calcular la integral del valor absoluto de la función, o bien descomponer el intervalo de integración en subintervalos donde la función sea positiva o negativa, integrar la función en cada uno de estos subintervalos y luego sumar el valor absoluto de estas integrales.
Ejemplo 13.9 Veamos cómo calcular el volumen de la función en el intervalo . En el intervalo de integración, la función toma valores negativos en el subintervalo y valores positivos en el intervalo , por lo que el volumen encerrado entre la gráfica de la función y el plano viene dado por la suma de las integrales
y
es decir,
Obsérvese que también se puede usar la integral doble para calcular el área de una región plana integrando sobre esta región la función , ya que cuanto la altura , el volumen coincide con el área de la base.
Del mismo modo se puede usar una integral tripe para calcular el volumen de una región tridimensional.
Cálculo del area de una superficie
Otra aplicación menos intuitiva de la las integrales dobles es el cálculo del área de la superficie de una función sobre una región . Si es una función de dos variables , con derivadas parciales continuas sobre un intervalo bidimensional , para aproximar el área de la superficie de sobre podemos tomar una partición del intervalo en subintervalos de igual tamaño, y otra del intervalo , y descomponer el intervalo en subintervalos rectangulares mediante la partición . De esta manera, todos los subintervalos tendrán el mismo área .
Si para cada subintervalo , podemos calcular el plano tangente a la superficie de en el punto , el paralelogramo que define el plano tangente sobre el intervalo nos dará una aproximación del área de la superficie de sobre el subintervalo , y por tanto, la suma de las áreas de estos paralelogramos tangentes será una aproximación de la superficie de sobre el intervalo .
Si para cada paralelogramo consideramos los vectores y con origen en y con la misma dirección y longitud que los lados del paralelogramo, tal y como se aprecia en la gráfica de más abajo, se tiene que el área del paralelogramo es , tal y como se vió en la Proposición 10.5.
Por otro lado, como está contenido en la intersección plano tangente a en , con el plano , se tiene que su pendiente es , por lo que sus coordenadas son . Del mismo modo, como está contenido en la intersección plano tangente, con el plano , se tiene que su pendiente es , por lo que sus coordenadas son .
Así pues, el producto vectorial de y es
y el área del paralelogramo es
De este modo, el área de la superficie de sobre el intervalo se puede aproximar mediante la suma
Tomando particiones cada vez más refinadas, en el límite cuando se tiene que el área de la superficie de sobre el intervalo viene dado por la integral
Ejemplo 13.10 La gráfica de la función es la semiesfera de radio centrada en el origen de coordenadas. El área de la superficie de esta semiesfera puede calcularse mediante la integral
donde .
Como las derivadas parciales de valen
la integral anterior se puede calcular mediante la integral iterada
Esta integral, no obstante, es difícil de calcular, por lo que merece la pena hacer el cambio a coordenadas polares, de manera que la región de integración polar es y la integral se convierte en
Esta última integral es impropia ya que cuando , pero la integral es finita, ya que
y siguiendo con el cálculo de la integral anterior se tiene
por lo que el área de la esfera de radio centrada en el origen de coordenadas es , que coincide con la fórmula habitual.
Aplicaciones físicas
Centro de masas
Una aplicación física inmediata de las integrales dobles o triples es el cálculo de la masa de una sólido plano o del espacio con densidad variable en cada punto. Si es la densidad en cada uno de los puntos de un sólido que ocupa una región del plano , la masa total del sólido puede calcularse mediante la integral
siempre y cuando la función sea integrable en la región .
Ejemplo 13.11 Una placa metálica tiene la forma de región y la densidad en cada punto viene dada por la función en gr/cm. Entonces su masa es
También podemos utilizar las integrales dobles para calcular el momento de cada punto de un sólido del plano con respecto a cada uno de los ejes. Tal y como se vio en la Sección 8.11.3.2, el momento con respecto a cada eje en un punto de un sólido plano viene dado por el producto de la masa en el punto y la distancia del punto al eje, es decir, el momento con respecto al eje es y el momento con respecto al eje es , donde es la masa en el punto . La suma de los momentos en cada uno de los puntos del sólido región nos dará el momento de todo el sólido. Así pues, si es la densidad en cada uno de los puntos de un sólido plano , los momentos de con respecto a los ejes e pueden calcularse, respectivamente, mediante las integrales
A partir de los momentos resulta sencillo calcular el centro de masas o centroide del sólido plano, dividiendo cada momento por la masa del sólido, es decir,
Un sólido plano se comporta como si toda su masa se concentrara en el centroide, de manera que se mantendría en equilibrio horizontal al apoyarlo sobre este punto.
Ejemplo 13.12 El momento con respecto al eje de la placa metálica del ejemplo anterior es
Del mismo modo, el momento con respecto al eje es
Por lo tanto, el centroide de la placa metálica tendrá coordenadas
Aplicaciones estadísticas
Valor medio de una función
En la Sección 8.12.1 vimos cómo utilizar la integral definida para calcular el valor medio de una función de una variable en un intervalo unidimensional. Para funciones de dos variables podemos usar la misma idea para calcular el valor medio de una función sobre una región dividiendo la integral por el área de la región de integración, y para funciones de tres variables dividiendo por el volumen de la región de integración. En general, para calcular el valor medio de una función de varias variables en una región de utilizaremos la integral de la siguiente definición.
Definición 13.6 (Valor medio de una función de varias variables) Dada una función integrable sobre una región , el valor medio de sobre es
Obsérvese que la integral del denominador , para una función de dos variables da el área de la región , mientras que para una función de tres variables da el volumen de la región .
Ejemplo 13.13 En el Ejemplo 13.5 vimos que la integral doble de la función sobre la región era . Como el área de la región es
el valor medio de en es