En matemáticas, específicamente en cálculo multivariable, una integral múltiple es un tipo de integral definida de una función de varias variables, por ejemplo, f ( x , y ) {\displaystyle f(x,y)} o f ( x , y , z ) {\displaystyle f(x,y,z)} . Integrales de funciones de dos variables sobre una región en R 2 {\displaystyle \mathbb {R} ^{2}} son llamadas integrales dobles mientras que integrales de funciones de tres variables sobre una región en R 3 {\displaystyle \mathbb {R} ^{3}} son llamadas integrales triples.
Así como la integral de una función positiva de una variable se interpreta como el área entre la gráfica de la función y el eje x {\displaystyle x}
, la integral doble de una función positiva de dos variables representa el volumen de la región entre la superficie definida por la función y el plano que contiene el dominio de la función. Para funciones de más de dos variables, la interpretación geométrica de la integral múltiple corresponde a hipervolúmenes.La integral múltiple de una función de n {\displaystyle n} 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 es abreviado por una letra en el signo de integral de más a la derecha:
∫ ⋯ ∫ D f ( x 1 , … , x n ) d x 1 ⋯ d x n . {\displaystyle \int \cdots \int _{D}f(x_{1},\dots ,x_{n})\;dx_{1}\cdots dx_{n}.} variables: f ( x 1 , x 2 , … , x n ) {\displaystyle f(x_{1},x_{2},\dots ,x_{n})} sobre un dominio D {\displaystyle D} típicamente es representada anidandoEs 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 x n + 1 = f ( x 1 , . . . , x n ) {\displaystyle x_{n+1}=f(x_{1},...,x_{n})} y una región T {\displaystyle T} en el espacio definido por los ejes de las variables independientes de la función f {\displaystyle f} (si T {\displaystyle T} es una región cerrada y acotada y f {\displaystyle f} está definida en ésta). Por ejemplo, si n = 2 {\displaystyle n=2} , el volumen situado entre la superficie definida por x 3 = f ( x 1 , x 2 ) {\displaystyle x_{3}=f(x_{1},x_{2})} y una región T {\displaystyle T} en el plano x 1 x 2 {\displaystyle x_{1}x_{2}} es igual a alguna integral doble, si es que, como se mencionó, f {\displaystyle f} está definida en T {\displaystyle T} .
T {\displaystyle T}
puede dividirse en una partición interior Δ {\displaystyle \Delta } formada por m {\displaystyle m} subregiones rectangulares sin solapamiento que estén completamente contenidas en T {\displaystyle T} . La norma | | Δ | | {\displaystyle ||\Delta ||} de esta partición está dada por la diagonal más larga en las m {\displaystyle m} subregiones.Si se toma un punto ( x 1 i , x 2 i , . . . , x n i ) {\displaystyle (x_{1i},x_{2i},...,x_{ni})}
f ( x 1 i , x 2 i , … , x n i ) Δ A i = f ( x 1 i , x 2 i , … , x n i ) Δ x 1 i Δ x 2 i … Δ x n i {\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}} que esté contenido dentro de la subregión con dimensiones Δ x 1 i Δ x 2 i . . . Δ x n i {\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 x n + 1 = f ( x 1 , . . . , x n ) {\displaystyle x_{n+1}=f(x_{1},...,x_{n})} y la subregión i. Este espacio tendrá una magnitud de:Entonces se puede aproximar la magnitud del espacio entero situado entre el objeto definido por la ecuación x n + 1 = f ( x 1 , . . . , x n ) {\displaystyle x_{n+1}=f(x_{1},...,x_{n})} suma de Riemann de las magnitudes de los m {\displaystyle m} espacios correspondientes a cada una de las subregiones:
∑ i = 1 m f ( x 1 i , x 2 i , … , x n i ) Δ x 1 i Δ x 2 i … Δ x n i {\displaystyle \sum _{i=1}^{m}f(x_{1i},x_{2i},\ldots ,x_{ni})\Delta x_{1i}\Delta x_{2i}\ldots \Delta x_{ni}} y la región T {\displaystyle T} mediante laEsta aproximación mejora a medida que el número m {\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:
lim m → ∞ ∑ i = 1 m f ( x 1 i , x 2 i , … , x n i ) Δ x 1 i Δ x 2 i … Δ x n i = lim ‖ Δ ‖ → 0 ∑ i = 1 m f ( x 1 i , x 2 i , … , x n i ) Δ x 1 i Δ x 2 i … Δ x n i {\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 este último límite es que el límite es igual L si y solo si para todo ε > 0 {\displaystyle \varepsilon >0}
| L − ∑ i = 1 m f ( x 1 i , x 2 i , … , x n i ) Δ x 1 i Δ x 2 i … Δ x n i | < ε {\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 } existe un δ > 0 {\displaystyle \delta >0} tal quepara toda partición Δ {\displaystyle \Delta }
de la región T {\displaystyle T} (que satisfaga | | Δ | | < δ {\displaystyle ||\Delta ||<\delta } ), y para todas las elecciones posibles de ( x 1 i , x 2 i , . . . , x n i ) {\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 f {\displaystyle f}
∫ ⋯ ∫ T f ( x 1 , … , x n ) d x 1 ⋯ d x n = lim ‖ Δ ‖ → 0 ∑ i = 1 m f ( x 1 i , x 2 i , … , x n i ) Δ x 1 i Δ x 2 i … Δ x n i {\displaystyle \int \cdots \int _{T}f(x_{1},\dots ,x_{n})\;dx_{1}\cdots 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}} está definida en una región cerrada y acotada T {\displaystyle T} del definido por los ejes de las variables independientes de f, la integral de f {\displaystyle f} sobre T {\displaystyle T} está dada por:siempre que el límite exista. Si el límite existe se dice que f {\displaystyle f}
es integrable con respecto a T.Es común que
∫ ⋯ ∫ T f ( x 1 , … , x n ) d x 1 ⋯ d x n {\displaystyle \int \cdots \int _{T}f(x_{1},\dots ,x_{n})\;dx_{1}\cdots dx_{n}}se denote por
∫ T f ( x ) d n x {\displaystyle \int _{T}f(\mathbf {x} )\;d^{n}\mathbf {x} }Si f {\displaystyle f} funciones continuas en una región cerrada y acotada D ⊂ R n {\displaystyle D\subset \mathbb {R} ^{n}} y c ∈ R {\displaystyle c\in \mathbb {R} } entonces la integral múltiple satisface algunas propiedades
y g {\displaystyle g} sonLa integral es un operador lineal pues satisface
∫ ⋯ ∫ D c f ( x 1 , … , x n ) d x 1 ⋯ d x n = c ∫ ⋯ ∫ D f ( x 1 , … , x n ) d x 1 ⋯ d x n {\displaystyle \int \cdots \int _{D}cf(x_{1},\dots ,x_{n})\;dx_{1}\cdots dx_{n}=c\int \cdots \int _{D}f(x_{1},\dots ,x_{n})\;dx_{1}\cdots dx_{n}} ∫ ⋯ ∫ D d x 1 ⋯ d x n = ∫ ⋯ ∫ D f ( x 1 , … , x n ) d x 1 ⋯ d x n ± ∫ ⋯ ∫ D g ( x 1 … , x n ) d x 1 ⋯ d x n {\displaystyle \int \cdots \int _{D}dx_{1}\cdots dx_{n}=\int \cdots \int _{D}f(x_{1},\dots ,x_{n})dx_{1}\cdots dx_{n}\pm \int \cdots \int _{D}g(x_{1}\dots ,x_{n})dx_{1}\cdots dx_{n}}1. Si f ( x 1 , … , x n ) ≥ 0 {\displaystyle f(x_{1},\dots ,x_{n})\geq 0}
∫ ⋯ ∫ D f ( x 1 , … , x n ) d x 1 ⋯ d x n ≥ 0 {\displaystyle \int \cdots \int _{D}f(x_{1},\dots ,x_{n})\;dx_{1}\cdots dx_{n}\geq 0} entonces2. Si f ( x 1 , … , x n ) ≥ g ( x 1 , … , x n ) {\displaystyle f(x_{1},\dots ,x_{n})\geq g(x_{1},\dots ,x_{n})}
∫ ⋯ ∫ D f ( x 1 , … , x n ) d x 1 ⋯ d x n ≥ ∫ ⋯ ∫ D g ( x 1 , … , x n ) d x 1 ⋯ d x n {\displaystyle \int \cdots \int _{D}f(x_{1},\dots ,x_{n})\;dx_{1}\cdots dx_{n}\geq \int \cdots \int _{D}g(x_{1},\dots ,x_{n})\;dx_{1}\cdots dx_{n}} entonces:3. Si D {\displaystyle D}
∫ ⋯ ∫ D f ( x 1 , … , x n ) d x 1 ⋯ d x n = ∫ ⋯ ∫ D 1 f ( x 1 , … , x n ) d x 1 ⋯ d x n + ∫ ⋯ ∫ D 2 f ( x 1 , … , x n ) d x 1 ⋯ d x n {\displaystyle \int \cdots \int _{D}f(x_{1},\dots ,x_{n})\;dx_{1}\cdots dx_{n}=\int \cdots \int _{D_{1}}f(x_{1},\dots ,x_{n})dx_{1}\cdots dx_{n}+\int \cdots \int _{D_{2}}f(x_{1},\dots ,x_{n})dx_{1}\cdots dx_{n}} es la unión entre las regiones D 1 {\displaystyle D_{1}} y D 2 {\displaystyle D_{2}} que no solapan entre sí entonces:Las integrales múltiples están estrechamente relacionadas con las integrales iteradas, las cuales son necesarias para resolver las integrales múltiples. La diferencia entre integrales múltiples e iteradas consiste en que una se refiere al concepto matemático de integral (aplicado a varias variables) y otra al procedimiento por el cual se resuelve la integral múltiple. Si l}
∫ a b ∫ c d f ( x , y ) d y d x {\displaystyle \int _{a}^{b}\int _{c}^{d}f(x,y)\,dy\,dx}se refiere a una integral iterada, la parte externa
∫ a b ⋯ d x {\displaystyle \int _{a}^{b}\cdots \,dx}es la integral con respecto a x de la función de x:
g ( x ) = ∫ c d f ( x , y ) d y . {\displaystyle g(x)=\int _{c}^{d}f(x,y)\,dy.}Una integral doble, en cambio está definida con respecto a un área en el plano xy. La integral doble existe si y solo si las dos integrales iteradas existen y son iguales. Es decir, si la integral doble existe, entonces es igual a la integral iterada, sin importar si el orden de integración es d y d x {\displaystyle dydx}
∫ a b ∫ c d f ( x , y ) d y d x ≠ ∫ c d ∫ a b f ( x , y ) d x d y . {\displaystyle \int _{a}^{b}\int _{c}^{d}f(x,y)\,dy\,dx\neq \int _{c}^{d}\int _{a}^{b}f(x,y)\,dx\,dy.} o d x d y {\displaystyle dxdy} , y por lo general uno la calcula calculando una sola de estas. Sin embargo, a veces las dos integrales iteradas existen sin ser iguales y en este caso no existe la integral doble, ya que se tiene:De una manera más formal, el Teorema de Fubini afirma que
∫ A × B | f ( x , y ) | d ( x , y ) < ∞ , {\displaystyle \int _{A\times B}|f(x,y)|\,d(x,y)<\infty ,}Esto es, si la integral es absolutamente convergente, entonces la integral doble es igual a la integral iterada.
∫ A × B f ( x , y ) d ( x , y ) = ∫ A ( ∫ B f ( x , y ) d y ) d x = ∫ B ( ∫ A f ( x , y ) d x ) d y . {\displaystyle \int _{A\times B}f(x,y)\,d(x,y)=\int _{A}\left(\int _{B}f(x,y)\,dy\right)\,dx=\int _{B}\left(\int _{A}f(x,y)\,dx\right)\,dy.}Esto ocurre, cuando f {\displaystyle f} función acotada y tanto A como B son regiones acotadas también. Esto se entiende fácilmente pensando que si la función o la región del dominio no están acotadas, la integral múltiple no puede existir.
es unaLa notación
∫ × f ( x , y ) d x d y {\displaystyle \int _{\times }f(x,y)\,dx\,dy}se puede usar si se desea ser enfático al referirse a una integral doble y no a una iterada.
Cuando el integrando es una función constante c {\displaystyle c}
, la integral es igual al producto entre c {\displaystyle c} y la integral sobre la región de integración.En particular si c = 1 {\displaystyle c=1} subconjunto en R 2 {\displaystyle \mathbb {R} ^{2}} entonces la integral doble calcula el área de la región, se denota por A ( D ) {\displaystyle A(D)} y está dada por
A ( D ) = ∬ D d A {\displaystyle A(D)=\iint _{D}dA} y la región de integración D {\displaystyle D} es unmientras que si se la región de integración D {\displaystyle D}
V ( D ) = ∭ D d A {\displaystyle V(D)=\iiint _{D}dA} es un subconjunto en R 3 {\displaystyle \mathbb {R} ^{3}} entonces obtenemos el volumen de la región, se denota por V ( D ) {\displaystyle V(D)} y está dada por EjemploConsidere f ( x , y ) = c {\displaystyle f(x,y)=c\,\!}
D = { ( x , y ) ∈ R 2 : 2 ≤ x ≤ 4 ; 3 ≤ y ≤ 6 } {\displaystyle D=\{(x,y)\in \mathbb {R} ^{2}:2\leq x\leq 4\;;\;3\leq y\leq 6\}} yintegrando f {\displaystyle f}
∬ D c d A = ∫ 3 6 ∫ 2 4 c d x d y = c ∫ 3 6 ∫ 2 4 d x d y = c A ( D ) = c ( 3 × 2 ) = 6 c {\displaystyle {\begin{aligned}\iint _{D}c\;dA&=\int _{3}^{6}\int _{2}^{4}c\;dxdy\\&=c\int _{3}^{6}\int _{2}^{4}dxdy\\&=c\;A(D)\\&=c\;(3\times 2)=6c\end{aligned}}} sobre D {\displaystyle D}En el caso de un dominio en el que exista simetría al menos respecto de uno de los ejes, y donde la función para integrar contenga al menos una función impar con respecto a esa variable, la integral se vuelve nula (ya que la suma de cantidades iguales con signo opuesto es cero).
EjemploConsidere f ( x , y ) = 2 sen ( x ) − 3 y 3 + 5 {\displaystyle f(x,y)=2\operatorname {sen}(x)-3y^{3}+5}
T = { ( x , y ) ∈ R 2 : x 2 + y 2 ≤ 1 } {\displaystyle T=\left\{(x,y)\in \mathbb {R} ^{2}:x^{2}+y^{2}\leq 1\right\}} yesto es T {\displaystyle T}
consiste en los puntos dentro de una circunferencia centrada en el origen de radio 1.Usando la propiedad de linealidad de las integrales, la integral puede ser escrita como
∬ T ( 2 sen ( x ) − 3 y 3 + 5 ) d x d y = ∬ T 2 sen ( x ) d x d y − ∬ T 3 y 3 d x d y + ∬ T 5 d x d y {\displaystyle \iint _{T}(2\operatorname {sen}(x)-3y^{3}+5)\;dxdy=\iint _{T}2\operatorname {sen}(x)\;dxdy-\iint _{T}3y^{3}\;dxdy+\iint _{T}5\;dxdy}Dado que 2 sen ( x ) {\displaystyle 2\operatorname {sen}(x)}
∬ T ( 2 sen ( x ) − 3 y 3 + 5 ) d x d y = ∬ T 2 sen ( x ) d x d y − ∬ T 3 y 3 d x d y + ∬ T 5 d x d y = ∬ T 5 d x d y = 5 A ( T ) = 5 π {\displaystyle {\begin{aligned}\iint _{T}(2\operatorname {sen}(x)-3y^{3}+5)\;dxdy&=\iint _{T}2\operatorname {sen}(x)\;dxdy-\iint _{T}3y^{3}\;dxdy+\iint _{T}5\;dxdy\\&=\iint _{T}5\;dxdy\\&=5A(T)=5\pi \end{aligned}}} y 3 y 3 {\displaystyle 3y^{3}} son funciones impares y existe simetría tanto con respecto al eje x {\displaystyle x} como con respecto al eje y {\displaystyle y} entonces las primeras dos integrales valen, por lo que la integral original se simplifica a sólo la tercera.A menudo, es útil para reducir la complejidad de la integral cambiar una variable por otra que resulte más cómoda, sin embargo esto exige el cambio de la región de integración, además de añadir un factor de corrección al diferencial conocido como determinante jacobiano (en valor absoluto o módulo). El cambio de una variable por otra es en un sentido geométrico, una transformación desde un espacio hasta otro, y es esta transformación la que exige estos ajustes.
Si se utiliza una transformación que siga la relación:
f ( y 1 , … , y n ) → f ( x 1 ( y 1 , y 2 , … , y n ) , … , x n ( y 1 , y 2 , … , y n ) ) {\displaystyle f(y_{1},\ldots ,y_{n})\rightarrow f(x_{1}(y_{1},y_{2},\ldots ,y_{n}),\ldots ,x_{n}(y_{1},y_{2},\ldots ,y_{n}))}Entonces se puede utilizar el jacobiano de la transformación para simplificar la integral
J = D ( y 1 , … , y n ) D ( x 1 , … , x n ) = | ∂ y 1 ∂ x 1 ⋯ ∂ y 1 ∂ x n ⋮ ⋱ ⋮ ∂ y n ∂ x 1 ⋯ ∂ y n ∂ x n | {\displaystyle J={\frac {D(y_{1},\ldots ,y_{n})}{D(x_{1},\ldots ,x_{n})}}={\begin{vmatrix}\displaystyle {\frac {\partial y_{1}}{\partial x_{1}}}&\cdots &\displaystyle {\frac {\partial y_{1}}{\partial x_{n}}}\\\vdots &\ddots &\vdots \\\displaystyle {\frac {\partial y_{n}}{\partial x_{1}}}&\cdots &\displaystyle {\frac {\partial y_{n}}{\partial x_{n}}}\end{vmatrix}}}Integrando la función transformada en el dominio de integración correspondiente a las variables x, y multiplicando por el valor absoluto del determinante jacobiano y por la serie de diferenciales, se obtiene una integral múltiple que es igual a la integral original, si es que esta existe.
∫ ⋯ ∫ D f ( y 1 , … , y n ) d y 1 ⋯ d y n = ∫ ⋯ ∫ T f ( x 1 , … , x n ) | J | d x 1 ⋯ d x n {\displaystyle \int \cdots \int _{D}f(y_{1},\dots ,y_{n})\;dy_{1}\cdots dy_{n}=\int \cdots \int _{T}f(x_{1},\dots ,x_{n})|J|dx_{1}\cdots dx_{n}}A continuación se dan algunos ejemplos de estas transformaciones.
Coordenadas PolaresEn R 2 {\displaystyle \mathbb {R} ^{2}} coordenadas polares, lo que significa que para cada punto genérico P ( x , y ) {\displaystyle P(x,y)} en coordenadas cartesianas cambia a su respectivo punto en las coordenadas polares para simplificar.
, si una región de integración tiene una simetría circular y la función tiene algunas características particulares entonces uno puede aplicar la transformación de coordenadas rectangulares aLa relación para llevar a cabo esta transformación es la siguiente:
f ( x , y ) → f ( r cos θ , r sen θ ) {\displaystyle f(x,y)\to f(r\cos \theta ,r\operatorname {sen} \theta )}esto es
x = r cos θ y = r sen θ {\displaystyle {\begin{aligned}x&=r\cos \theta \\y&=r\operatorname {sen} \theta \end{aligned}}}donde r ≥ 0 {\displaystyle r\geq 0}
y 0 ≤ θ ≤ 2 π {\displaystyle 0\leq \theta \leq 2\pi } .De lo anterior se tiene que
x 2 + y 2 = ( r cos θ ) 2 + ( r sen θ ) 2 = r 2 cos 2 θ + r 2 sen 2 θ = r 2 ( cos 2 θ + sen 2 θ ) = r 2 {\displaystyle {\begin{aligned}x^{2}+y^{2}&=(r\cos \theta )^{2}+(r\operatorname {sen} \theta )^{2}\\&=r^{2}\cos ^{2}\theta +r^{2}\operatorname {sen} ^{2}\theta \\&=r^{2}\left(\cos ^{2}\theta +\operatorname {sen} ^{2}\theta \right)\\&=r^{2}\end{aligned}}}El determinante jacobiano de la transformación es:
∂ ( x , y ) ∂ ( r , θ ) = | cos θ − r sen θ sen θ r cos θ | = r {\displaystyle {\frac {\partial (x,y)}{\partial (r,\theta )}}={\begin{vmatrix}\cos \theta &-r\operatorname {sen} \theta \\\operatorname {sen} \theta &r\cos \theta \end{vmatrix}}=r}el cual fue obtenido insertando las derivadas parciales de x = r cos θ {\displaystyle x=r\cos \theta }
y y = r sen θ {\displaystyle y=r\operatorname {sen} \theta } en la primera columna con respecto a r {\displaystyle r} y en la segunda con respecto a θ {\displaystyle \theta } , por lo que los diferenciales d x d y {\displaystyle dxdy} se transforman en r d r d θ {\displaystyle rdrd\theta }Una vez transformada la función y la región de integración, es posible definir una fórmula para el cambio de variables en coordenadas polares:
∬ D f ( x , y ) d A = ∬ T f ( r cos θ , r sen θ ) r d A ′ {\displaystyle \iint _{D}f(x,y)\;dA=\iint _{T}f(r\cos \theta ,r\operatorname {sen} \theta )r\;dA'}donde d A ′ = d r d θ {\displaystyle dA'=drd\theta }
∬ D f ( x , y ) d A = ∬ T f ( ρ cos θ , ρ sen θ ) ρ d A ′ {\displaystyle \iint _{D}f(x,y)\;dA=\iint _{T}f(\rho \cos \theta ,\rho \operatorname {sen} \theta )\rho \;dA'} o d A ′ = d θ d r {\displaystyle dA'=d\theta dr} . En ocasiones en lugar de utilizar r {\displaystyle r} , es común utilizar ρ {\displaystyle \rho } en cuyo caso obtendríamosPor ejemplo, si la función es f ( x , y ) = x + y {\displaystyle f(x,y)=x+y}
f ( r , θ ) = r cos θ + r sen θ = r ( cos θ + sen θ ) {\displaystyle f(r,\theta )=r\cos \theta +r\operatorname {sen} \theta =r(\cos \theta +\operatorname {sen} \theta )} entonces aplicando la transformación se obtiene EjemploConsidere la región
D = { ( x , y ) ∈ R 2 : x 2 + y 2 ≤ 9 , x 2 + y 2 ≥ 4 , y ≥ 0 } {\displaystyle D=\left\{(x,y)\in \mathbb {R} ^{2}:x^{2}+y^{2}\leq 9,x^{2}+y^{2}\geq 4,y\geq 0\right\}}aplicando la transformación obtenemos la siguiente región
T = { ( r , θ ) : 2 ≤ r ≤ 3 , 0 ≤ θ ≤ 2 π } {\displaystyle T=\left\{(r,\theta ):2\leq r\leq 3,0\leq \theta \leq 2\pi \right\}}Si f ( x , y ) = x {\displaystyle f(x,y)=x}
∬ D x d A = ∬ T r cos θ r d A ′ = ∫ 0 2 π ∫ 2 3 r 2 cos θ d r d θ = ∫ 0 2 π cos θ d θ ∫ 2 3 r 2 d r = 0 {\displaystyle {\begin{aligned}\iint _{D}x\;dA&=\iint _{T}r\cos \theta \;r\;dA'\\&=\int _{0}^{2\pi }\int _{2}^{3}r^{2}\cos \theta \;drd\theta \\&=\int _{0}^{2\pi }\cos \theta \;d\theta \int _{2}^{3}r^{2}dr\\&=0\end{aligned}}} entonces Coordenadas CilíndricasEn R 3 {\displaystyle \mathbb {R} ^{3}} coordenadas cilíndricas. La función se transforma mediante la siguiente relación:
f ( x , y , z ) → f ( r cos θ , r sen θ , z ) {\displaystyle f(x,y,z)\to f(r\cos \theta ,r\operatorname {sen} \theta ,z)} , la integración sobre regiones con base circular puede ser hecha transformando aesto es
x = r cos θ y = r sen θ z = z {\displaystyle {\begin{aligned}x&=r\cos \theta \\y&=r\operatorname {sen} \theta \\z&=z\end{aligned}}}con r ≥ 0 {\displaystyle r\geq 0}
, 0 ≤ θ ≤ 2 π {\displaystyle 0\leq \theta \leq 2\pi } y z ≥ 0 {\displaystyle z\geq 0} .De lo anterior se tiene que
x 2 + y 2 = ( r cos θ ) 2 + ( r sen θ ) 2 = r 2 cos 2 θ + r 2 sen 2 θ = r 2 ( cos 2 θ + sen 2 θ ) = r 2 {\displaystyle {\begin{aligned}x^{2}+y^{2}&=(r\cos \theta )^{2}+(r\operatorname {sen} \theta )^{2}\\&=r^{2}\cos ^{2}\theta +r^{2}\operatorname {sen} ^{2}\theta \\&=r^{2}\left(\cos ^{2}\theta +\operatorname {sen} ^{2}\theta \right)\\&=r^{2}\end{aligned}}}El determinante jacobiano de la transformación es:
∂ ( x , y , z ) ∂ ( r , θ , z ) = | cos θ − r sen θ 0 sen θ r cos θ 0 0 0 1 | = r {\displaystyle {\frac {\partial (x,y,z)}{\partial (r,\theta ,z)}}={\begin{vmatrix}\cos \theta &-r\operatorname {sen} \theta &0\\\operatorname {sen} \theta &r\cos \theta &0\\0&0&1\end{vmatrix}}=r}Una vez transformada la función y la región de integración, es posible definir una fórmula para el cambio de variables en coordenadas cilíndricas:
∭ D f ( x , y , z ) d V = ∭ T f ( r cos θ , r sen θ , z ) r d V ′ {\displaystyle \iiint _{D}f(x,y,z)\;dV=\iiint _{T}f(r\cos \theta ,r\operatorname {sen} \theta ,z)r\;dV'}donde d V ′ {\displaystyle dV'}
∭ D f ( x , y , z ) d V = ∭ T f ( ρ cos θ , ρ sen θ , z ) r d V ′ {\displaystyle \iiint _{D}f(x,y,z)\;dV=\iiint _{T}f(\rho \cos \theta ,\rho \operatorname {sen} \theta ,z)r\;dV'} representa los posibles órdenes de integración. En ocasiones suele utilizarse ρ {\displaystyle \rho } en lugar de r {\displaystyle r} en cuyo caso tendríamos EjemploConsidere la región
D = { ( x , y , z ) ∈ R 3 : x 2 + y 2 ≤ 9 , − 5 ≤ z ≤ 5 } {\displaystyle D=\left\{(x,y,z)\in \mathbb {R} ^{3}:x^{2}+y^{2}\leq 9,-5\leq z\leq 5\right\}}utilizando la transformación obtenemos
T = { ( r , θ , z ) : 0 ≤ r ≤ 3 , 0 ≤ θ ≤ 2 π , − 5 ≤ z ≤ 5 } {\displaystyle T=\left\{(r,\theta ,z):0\leq r\leq 3,0\leq \theta \leq 2\pi ,-5\leq z\leq 5\right\}}si f ( x , y , z ) = x 2 + y 2 + z {\displaystyle f(x,y,z)=x^{2}+y^{2}+z}
∭ D ( x 2 + y 2 + z ) d V = ∭ T ( r 2 + z ) r d V ′ = ∫ − 5 5 ∫ 0 2 π ∫ 0 3 ( r 2 + z ) r d r d θ d z = ∫ 0 2 π d θ ∫ − 5 5 ∫ 0 3 ( r 3 + r z ) d r d z = 405 π {\displaystyle {\begin{aligned}\iiint _{D}\left(x^{2}+y^{2}+z\right)dV&=\iiint _{T}\left(r^{2}+z\right)r\;dV'\\&=\int _{-5}^{5}\int _{0}^{2\pi }\int _{0}^{3}\left(r^{2}+z\right)r\;dr\;d\theta \;dz\\&=\int _{0}^{2\pi }d\theta \int _{-5}^{5}\int _{0}^{3}\left(r^{3}+rz\right)dr\;dz\\&=405\pi \end{aligned}}} entonces Coordenadas EsféricasCuando existe simetría esférica en un dominio en R 3 {\displaystyle \mathbb {R} ^{3}} coordenadas esféricas para simplificar una integral triple. La función es transformada por la relación:
f ( x , y , z ) → f ( ρ cos θ sen ϕ , ρ sen θ sen ϕ , ρ cos ϕ ) {\displaystyle f(x,y,z)\to f(\rho \cos \theta \operatorname {sen} \phi ,\rho \operatorname {sen} \theta \operatorname {sen} \phi ,\rho \cos \phi )} , es posible utilizar una transformación haciaesto es
x = ρ cos θ sen ϕ y = ρ sen θ sen ϕ z = ρ cos ϕ {\displaystyle {\begin{aligned}x&=\rho \cos \theta \operatorname {sen} \phi \\y&=\rho \operatorname {sen} \theta \operatorname {sen} \phi \\z&=\rho \cos \phi \end{aligned}}}con ρ ≥ 0 {\displaystyle \rho \geq 0}
, 0 ≤ θ ≤ 2 π {\displaystyle 0\leq \theta \leq 2\pi } y 0 ≤ ϕ ≤ π {\displaystyle 0\leq \phi \leq \pi } .De lo anterior se tiene que
x 2 + y 2 + z 2 = ( ρ cos θ sen ϕ ) 2 + ( ρ sen θ sen ϕ ) 2 + ( ρ cos ϕ ) 2 = ρ 2 cos 2 θ sen 2 ϕ + ρ 2 sen 2 θ sen 2 ϕ + ρ 2 cos 2 ϕ = ρ 2 ( cos 2 θ sen 2 ϕ + sen 2 θ sen 2 ϕ + cos 2 ϕ ) = ρ 2 ( sen 2 ϕ ( cos 2 θ + sen 2 θ ) + cos 2 ϕ ) = ρ 2 ( sen 2 ϕ + cos 2 ϕ ) = ρ 2 {\displaystyle {\begin{aligned}x^{2}+y^{2}+z^{2}&=(\rho \cos \theta \operatorname {sen} \phi )^{2}+(\rho \operatorname {sen} \theta \operatorname {sen} \phi )^{2}+(\rho \cos \phi )^{2}\\&=\rho ^{2}\cos ^{2}\theta \operatorname {sen} ^{2}\phi +\rho ^{2}\operatorname {sen} ^{2}\theta \operatorname {sen} ^{2}\phi +\rho ^{2}\cos ^{2}\phi \\&=\rho ^{2}\left(\cos ^{2}\theta \operatorname {sen} ^{2}\phi +\operatorname {sen} ^{2}\theta \operatorname {sen} ^{2}\phi +\cos ^{2}\phi \right)\\&=\rho ^{2}\left(\operatorname {sen} ^{2}\phi (\cos ^{2}\theta +\operatorname {sen} ^{2}\theta )+\cos ^{2}\phi \right)\\&=\rho ^{2}\left(\operatorname {sen} ^{2}\phi +\cos ^{2}\phi \right)\\&=\rho ^{2}\end{aligned}}}El determinante jacobiano de la transformación es el siguiente:
∂ ( x , y , z ) ∂ ( ρ , θ , ϕ ) = | ∂ x ∂ ρ ∂ x ∂ θ ∂ x ∂ ϕ ∂ y ∂ ρ ∂ y ∂ θ ∂ y ∂ ϕ ∂ z ∂ ρ ∂ z ∂ θ ∂ z ∂ ϕ | = | cos θ sen ϕ − ρ sen θ sen ϕ ρ cos θ cos ϕ sen θ sen ϕ ρ cos θ sen ϕ ρ sen θ cos ϕ cos ϕ 0 − ρ sen ϕ | = ρ 2 sen ϕ {\displaystyle {\begin{aligned}{\frac {\partial (x,y,z)}{\partial (\rho ,\theta ,\phi )}}&={\begin{vmatrix}{\frac {\partial x}{\partial \rho }}&{\frac {\partial x}{\partial \theta }}&{\frac {\partial x}{\partial \phi }}\\{\frac {\partial y}{\partial \rho }}&{\frac {\partial y}{\partial \theta }}&{\frac {\partial y}{\partial \phi }}\\{\frac {\partial z}{\partial \rho }}&{\frac {\partial z}{\partial \theta }}&{\frac {\partial z}{\partial \phi }}\end{vmatrix}}\\&={\begin{vmatrix}\cos \theta \operatorname {sen} \phi &-\rho \operatorname {sen} \theta \operatorname {sen} \phi &\rho \cos \theta \cos \phi \\\operatorname {sen} \theta \operatorname {sen} \phi &\rho \cos \theta \operatorname {sen} \phi &\rho \operatorname {sen} \theta \cos \phi \\\cos \phi &0&-\rho \operatorname {sen} \phi \end{vmatrix}}\\&=\rho ^{2}\operatorname {sen} \phi \end{aligned}}}Una vez transformada la función y la región de integración, es posible definir una fórmula para el cambio de variables en coordenadas esféricas:
∭ D f ( x , y , z ) d V = ∭ T f ( ρ cos θ sen ϕ , ρ sen θ sen ϕ , ρ cos ϕ ) ρ 2 sen ϕ d V ′ {\displaystyle \iiint _{D}f(x,y,z)\;dV=\iiint _{T}f(\rho \cos \theta \operatorname {sen} \phi ,\rho \operatorname {sen} \theta \operatorname {sen} \phi ,\rho \cos \phi )\rho ^{2}\operatorname {sen} \phi \;dV'} EjemploConsidere la región
D = { ( x , y , z ) ∈ R 3 : x 2 + y 2 + z 2 ≤ 1 } {\displaystyle D=\left\{(x,y,z)\in \mathbb {R} ^{3}:x^{2}+y^{2}+z^{2}\leq 1\right\}}utilizando la transformación obtenemos
T = { ( ρ , θ , φ ) : 0 ≤ ρ ≤ 1 , 0 ≤ θ ≤ 2 π , 0 ≤ φ ≤ π } {\displaystyle T=\left\{(\rho ,\theta ,\varphi ):0\leq \rho \leq 1,0\leq \theta \leq 2\pi ,0\leq \varphi \leq \pi \right\}}si f ( x , y , z ) = exp ( x 2 + y 2 + z 2 ) 3 {\displaystyle f(x,y,z)=\exp {\sqrt {(x^{2}+y^{2}+z^{2})^{3}}}}
∭ D exp ( x 2 + y 2 + z 2 ) 3 d V = ∭ T e ρ 3 ρ 2 sen φ d V ′ = ∫ 0 1 ∫ 0 2 π ∫ 0 π ρ 2 e ρ 3 sen φ d φ d θ d ρ = ∫ 0 2 π d θ ∫ 0 1 ρ 2 e ρ 3 d ρ ∫ 0 π sen φ d φ = 2 π ( e − 1 3 ) ( 2 ) = 4 π ( e − 1 ) 3 {\displaystyle {\begin{aligned}\iiint _{D}\exp {\sqrt {\left(x^{2}+y^{2}+z^{2}\right)^{3}}}\;dV&=\iiint _{T}e^{\rho ^{3}}\rho ^{2}\operatorname {sen} \varphi \;dV'\\&=\int _{0}^{1}\int _{0}^{2\pi }\int _{0}^{\pi }\rho ^{2}e^{\rho ^{3}}\operatorname {sen} \varphi \;d\varphi \;d\theta \;d\rho \\&=\int _{0}^{2\pi }d\theta \int _{0}^{1}\rho ^{2}e^{\rho ^{3}}d\rho \int _{0}^{\pi }\operatorname {sen} \varphi \;d\varphi \\&=2\pi \left({\frac {e-1}{3}}\right)(2)\\&={\frac {4\pi (e-1)}{3}}\end{aligned}}} entoncesEl volumen de un cilindro con altura h {\displaystyle h}
D = { ( x , y , z ) ∈ R 3 : x 2 + y 2 ≤ R 2 , 0 ≤ z ≤ h } {\displaystyle D=\left\{(x,y,z)\in \mathbb {R} ^{3}:x^{2}+y^{2}\leq R^{2},0\leq z\leq h\right\}} y base circular de radio R {\displaystyle R} puede ser calculado con una integral triple utilizando coordenadas cilíndricas, si la región D {\displaystyle D} dada porrepresenta el cilindro entonces utilizando la transformación adecuada obtenemos la región T {\displaystyle T}
T = { ( r , θ , z ) : 0 ≤ θ ≤ 2 π , 0 ≤ r ≤ R , 0 ≤ z ≤ h } {\displaystyle T=\left\{(r,\theta ,z):0\leq \theta \leq 2\pi ,0\leq r\leq R,0\leq z\leq h\right\}} dada porpor lo que el volumen del cilindro puede calcularse como
V ( D ) = ∭ D d V = ∭ T r d V ′ = ∫ 0 2 π ∫ 0 R ∫ 0 h r d z d r d θ = ∫ 0 2 π d θ ∫ 0 R r d r ∫ 0 h d z = 2 π ( R 2 2 ) h = π R 2 h {\displaystyle {\begin{aligned}V(D)&=\iiint _{D}dV\\&=\iiint _{T}r\;dV'\\&=\int _{0}^{2\pi }\int _{0}^{R}\int _{0}^{h}r\;dz\;dr\;d\theta \\&=\int _{0}^{2\pi }d\theta \int _{0}^{R}r\;dr\int _{0}^{h}dz\\&=2\pi \left({\frac {R^{2}}{2}}\right)h\\&=\pi R^{2}h\end{aligned}}} EsferaPara demostrar que el volumen de una esfera de radio r {\displaystyle r}
G = { ( x , y , z ) ∈ R 3 : x 2 + y 2 + z 2 ≤ r 2 } {\displaystyle G=\left\{(x,y,z)\in \mathbb {R} ^{3}:x^{2}+y^{2}+z^{2}\leq r^{2}\right\}} es 4 π r 3 3 {\textstyle {\frac {4\pi r^{3}}{3}}} consideremos la región G ⊂ R 3 {\displaystyle G\subset \mathbb {R} ^{3}} dada poresto es, G {\displaystyle G}
S = { ( ρ , θ , φ ) : 0 ≤ ρ ≤ r , 0 ≤ θ ≤ 2 π , 0 ≤ φ ≤ π } {\displaystyle S=\left\{(\rho ,\theta ,\varphi ):0\leq \rho \leq r,0\leq \theta \leq 2\pi ,0\leq \varphi \leq \pi \right\}} contiene todos los puntos en R 3 {\displaystyle \mathbb {R} ^{3}} dentro de una esfera centrada en el origen y con radio r {\displaystyle r} . Si utilizamos coordenadas esféricas entonces obtendremos la región S {\displaystyle S} dada porpor lo que el volumen de una esfera puede calcularse como
V ( G ) = ∭ G d V = ∭ S ρ 2 sen φ d V ′ = ∫ 0 2 π ∫ 0 π ∫ 0 r ρ 2 sen φ d ρ d φ d θ = ∫ 0 2 π d θ ∫ 0 π sen φ d φ ∫ 0 r ρ 2 d ρ = 2 π ( − cos φ | 0 π ) ( ρ 3 3 | 0 r ) = 2 π ( 2 ) ( r 3 3 ) = 4 π r 3 3 {\displaystyle {\begin{aligned}V(G)&=\iiint _{G}dV\\&=\iiint _{S}\rho ^{2}\operatorname {sen} \varphi \;dV'\\&=\int _{0}^{2\pi }\int _{0}^{\pi }\int _{0}^{r}\rho ^{2}\operatorname {sen} \varphi \;d\rho \;d\varphi \;d\theta \\&=\int _{0}^{2\pi }d\theta \int _{0}^{\pi }\operatorname {sen} \varphi \;d\varphi \int _{0}^{r}\rho ^{2}d\rho \\&=2\pi \left(-\cos \varphi {\bigg |}_{0}^{\pi }\right)\left({\frac {\rho ^{3}}{3}}{\bigg |}_{0}^{r}\right)\\&=2\pi (2)\left({\frac {r^{3}}{3}}\right)\\&={\frac {4\pi r^{3}}{3}}\end{aligned}}}