Integrales Multiples I

De la misma manera en que la integral de una función positiva ${\displaystyle f(x)}$ de una variable definida en un intervalo puede interpretarse cómo el área entre la gráfica de la función y el eje x en ese intervalo, la doble integral de una función positiva ${\displaystyle f(x,y)}$ de dos variables, definida en una región del plano xy, se puede interpretar como el volumen entre la superficie definida por la función y el plano xy en ese intervalo. Al realizar una "integral triple" de una función ${\displaystyle f(x,y,z)}$ definida en una región del espacio xyz, el resultado es un hipervolumen, sin embargo es bueno notar que si ${\displaystyle f(x,y,z)=1}$ el resultado se puede interpretar como el volumen de la región de integración. Para integrales de órdenes superiores, el resultado geométrico corresponde a hipervolúmenes de dimensiones cada vez superiores.
La manera más usual de representar una integral múltiple es anidando signos de integración en el orden inverso al orden de ejecución (el de más a la izquierda es el último en ser calculado), seguido de la función y los diferenciales en orden de ejecución. El dominio de integración se representa sobre cada signo de integral, o a menudo es abreviado por una letra en el signo de integral de más a la derecha:

${\displaystyle \iint \ldots \int _{\mathbf {D} }\;f(x_{1},x_{2},\ldots ,x_{n})\;\mathbf {d} x_{1}\mathbf {d} x_{2}\!\ldots \mathbf {d} x_{n}}$

Es importante destacar que no es posible calcular la función primitiva o antiderivada de una función de más de una variable por lo que las integrales múltiples indefinidas no existen.

Una forma relativamente sencilla de definir las integrales múltiples es mediante su representación geométrica como la magnitud del espacio entre el objeto definido por la ecuación ${\displaystyle x_{n+1}=f(x_{1},...,x_{n})}$ y una región ${\displaystyle T}$ en el espacio definido por los ejes de las variables independientes de la función ${\displaystyle f}$ (si ${\displaystyle T}$ es una región cerrada y acotada y ${\displaystyle f}$ está definida en ésta). Por ejemplo, si ${\displaystyle n=2}$, el volumen situado entre la superficie definida por ${\displaystyle x_{3}=f(x_{1},x_{2})}$ y una región ${\displaystyle T}$ en el plano ${\displaystyle x_{1}x_{2}}$ es igual a alguna integral doble, si es que, como se mencionó, ${\displaystyle f}$ está definida en ${\displaystyle T}$.
${\displaystyle T}$ puede dividirse en una partición interior ${\displaystyle \Delta }$ formada por ${\displaystyle m}$ subregiones rectangulares sin solapamiento que estén completamente contenidas en ${\displaystyle T}$. La norma ${\displaystyle ||\Delta ||}$ de esta partición está dada por la diagonal más larga en las ${\displaystyle m}$ subregiones.
Si se toma un punto ${\displaystyle (x_{1i},x_{2i},...,x_{ni})}$ que esté contenido dentro de la subregión con dimensiones ${\displaystyle \Delta x_{1i}\Delta x_{2i}...\Delta x_{ni}}$ para cada una de las m subregiones de la partición, se puede construir un espacio con una magnitud aproximada a la del espacio entre el objeto definido por ${\displaystyle x_{n+1}=f(x_{1},...,x_{n})}$ y la subregión i. Este espacio tendrá una magnitud de:

${\displaystyle f(x_{1i},x_{2i},\ldots ,x_{ni})\Delta A_{i}=f(x_{1i},x_{2i},\ldots ,x_{ni})\Delta x_{1i}\Delta x_{2i}\ldots \Delta x_{ni}}$

Entonces se puede aproximar la magnitud del espacio entero situado entre el objeto definido por la ecuación ${\displaystyle x_{n+1}=f(x_{1},...,x_{n})}$ y la región${\displaystyle T}$mediante la suma de Riemann de las magnitudes de los ${\displaystyle m}$ espacios correspondientes a cada una de las subregiones:

${\displaystyle \sum _{i=1}^{m}f(x_{1i},x_{2i},\ldots ,x_{ni})\Delta x_{1i}\Delta x_{2i}\ldots \Delta x_{ni}}$

Esta aproximación mejora a medida que el número ${\displaystyle m}$ de subregiones se hace mayor. Esto sugiere que se podría obtener la magnitud exacta tomando el límite. Al aumentar el número de subregiones disminuirá la norma de la partición:

${\displaystyle \lim _{m\to \infty }\sum _{i=1}^{m}f(x_{1i},x_{2i},\ldots ,x_{ni})\Delta x_{1i}\Delta x_{2i}\ldots \Delta x_{ni}=\lim _{\|\Delta \|\to 0}\sum _{i=1}^{m}f(x_{1i},x_{2i},\ldots ,x_{ni})\Delta x_{1i}\Delta x_{2i}\ldots \Delta x_{ni}}$

El significado riguroso de éste último límite es que el límite es igual L si y sólo si para todo ${\displaystyle \varepsilon >0}$ existe un ${\displaystyle \delta >0}$ tal que

${\displaystyle \left|L-\sum _{i=1}^{m}f(x_{1i},x_{2i},\ldots ,x_{ni})\Delta x_{1i}\Delta x_{2i}\ldots \Delta x_{ni}\right\vert <\varepsilon }$

para toda partición ${\displaystyle \Delta }$ de la región${\displaystyle T}$(que satisfaga ${\displaystyle ||\Delta ||<\delta }$), y para todas las elecciones posibles de ${\displaystyle (x_{1i},x_{2i},...,x_{ni})}$ en la iésima subregión. Esto conduce a la definición formal de una integral múltiple:

Si${\displaystyle f}$está definida en una región cerrada y acotada${\displaystyle T}$del definido por los ejes de las variables independientes de f, la integral de${\displaystyle f}$sobre${\displaystyle T}$está dada por:

${\displaystyle \iint \ldots \int _{T}\,f(x_{1},\ldots ,x_{n})dx_{1}\ldots \,dx_{n}=\lim _{\|\Delta \|\to 0}\sum _{i=1}^{m}f(x_{1i},x_{2i},\ldots ,x_{ni})\Delta x_{1i}\Delta x_{2i}\ldots \Delta x_{ni}}$

siempre que el límite exista. Si el límite existe se dice que${\displaystyle f}$es integrable con respecto a T.