2. Integral doble

Sea R = [a,b] × [c,d], y sea f : 2 una función definida y acotada sobre R. Supongamos que MR = {R11,R12,...Rnm} es un conjunto de nm rectángulos que conforman una malla que cubre R (Figura 6.3). El área de cada celda Rij la denotamos con ΔAij. La malla MR es el conjunto de rectángulos Rij = [xi,xi+1] × [yj,yj+1] de área ΔAij = ΔxiΔyj.

PIC

Fig. 6.3: Malla de nm rectángulos

Una suma de Riemann de f sobre R es una expresión de la forma

i=1n j=1mf(x ij,y ij)ΔA ijdonde(xij,y ij) R ij

Si f es continua y positiva sobre R, entonces f(xij,yij)ΔAij aproxima el volumen de cada paralelepípedo Pij de base Rij y altura f(xij,yij); en este caso la suma de Riemann aproxima el volumen del sólido entre la región R (en el plano XY ) y el gráfico de f.

PIC

Fig. 6.4:

Diámetro de la malla. El diámetro de cada celda Rij es la máxima distancia entre todas las distancias entre cualesquiera dos puntos en Rij y se denota ||Rij||. El diámetro de la malla MR es ||MR|| = Supij{||Rij||}. Conforme ||MR|| 0, el área de cada celda tiende a cero, es decir, ΔAij 0 y la cantidad de celdas se hace infinitamente grande: n .

Volumen. Si f es continua y no negativa en la región R, en el pl;ano XY , entonces el siguiente límite existe,

V Q = lim ||M||0 i=1n j=1mf(x ij,y ij)ΔA ij con nm =  Card (M)

y V Q es el volumen de sólido Q, limitado por la región R en el plano XY , y la superficie S : z = f(x,y). El límite se toma sobre todas las posibles mallas rectángulares MR con (xij,yij) cualquier punto de Rij. En particular se puede tomar como (xij,yij) los extremos (xi,yj).

Propósito: Visualizar el volumen de un sólido como el límite de una suma de volúmenes de paralelepípedos

Caso general. Suponiendo que f es continua sobre una región R en el plano XY y que la región R está limitada por una curva suave a trozos (con un número finito de trozos), entonces conforme ||M|| 0, el área de cada celda tiende a cero (ΔAk 0) y la cantidad de celdas se hace infinitamente grande. La unión de los rectángulos Rij contenidos en R se va ajustando a la región R conforme la malla se refina.

El volumen de Q se obtiene como:

V = lim ||M||0 i=1n j=1mf(x ij,y ij)ΔA ijcon(xij,y ij) R ij

donde el límite se toma sobre todas las mallas M cuya norma tiende a cero, y (xij,yij) es cualquier punto de muestra dentro del rectángulo Rij (la función se debe definir como nula, fuera de la región R). La generalización de estas ideas no requiere que f sea no negativa; solo requiere que el límite exista (es decir, que f sea integrable).

PIC

Fig. 6.5:
Definición 21 — Función integrable.

Si las sumas de Riemann de f convergen a un límite único, independiente de la elección de los puntos de muestra (xij,yij) en cada sub-rectángulo, conforme la norma de la malla M 0 para cualquier sucesión de particiones, entonces decimos que f es integrable sobre R. En este caso, la integral doble se define como:

Rf(x,y)dA = lim M0 i=1n j=1mf(x ij,y ij)Δx iΔyj

donde el límite se toma sobre todas las particiones de la región R cuya norma tiende a cero.

Las propiedades de las funciones integrables en dos variables son similares a las propiedades de las funciones integrables en una variable.

Teorema 24 — Funciones integrables.

1
Si f es continua sobre R (cerrado y acotado), entonces f es integrable sobre R.
2
Si f es acotada y el conjunto de sus discontinuidades en R tiene área cero (en particular, si es continua excepto quizás en un conjunto finito de puntos o a lo largo de un número finito de curvas suaves, como segmentos de recta o arcos de círculo), entonces f es integrable sobre R.
3
Si f es integrable, entonces el valor del límite es el mismo para cualquier elección de los puntos muestra (xij,yij). En consecuencia, para calcular la integral (una vez que sabemos que existe), podemos elegir cualquier conjunto de puntos muestra, por ejemplo los extremos de cada subintervalo o los puntos medios, siempre que tomemos particiones cuya norma tienda a cero.

Si la región es un rectángulo R = [a,b] × [p,q] y se toman particiones regulares (subintervalos de igual longitud), entonces:

Δxi = b a n ,Δyj = q p m ,por lo que ΔAij = ΔxiΔyj = (b a)(q p) nm .

Si la región no es rectangular, debemos definir la función como nula fuera de la región y aplicar la misma fórmula (siempre que la frontera de la región tenga área cero).

Teorema 25 — Propiedades de la funciones integrables.

1
Si f y g son integrables sobre regiones R y S que no se traslapan, entonces f es integrables sobre R S y RSf(x,y)dA =Rf(x,y)dA + Sf(x,y)dA

2
Si f y g son integrables sobre R y f(x,y) g(x,y) para todo (x,y) R, entonces Rf(x,y)dA Rg(x,y)dA

3
Si f es integrable sobre R y M f(x,y) m para todo (x,y) R, entonces, si A(R) es el área de R, mA(R) Rf(x,y)dA MA(R)

La hipótesis "...son integrables sobre regiones R y S que no se traslapan.." es una condición suficiente, pero en realidad no hay problema si se superponen en conjuntos de "medida cero" (por ejemplo, una línea o un conjunto finito de puntos).

Ejemplo145Cálculo de integrales triples

[Aproximación de un volumen]

Consideremos el sólido Q limitado por la superficie S1 : z = 8 x2 y2 y el cilindro S2 : x2 + y2 4 en el primer octante, tal y como se muestra en la Figura 6.6. La región de integración R sería el cuarto de círculo de radio 2 en el primer cuadrante.
La función f(x,y) = 8 x2 y2 es integrable en esta región R. Podemos aproximar el volumen del sólido Q usando la aproximación

Rf(x,y)dA i=1n j=1mf(x ij,y ij)ΔA ij

PIC

Fig. 6.6:

La región de integración no es rectangular, pero para nuestra aproximación podemos tomar los puntos (xij,yij en el interior de la región.

Usando Δx = Δy = 0.5, tendríamos una malla regular de 8 rectángulos, cada uno de área 0.52.

PIC

Fig. 6.7:

Calculamos f(xij,yij) = f(xi,yi) para cada xi = i 0.5,yi = i 0.5 en la malla

xi yi f(xi,yi) xi yi f(xi,yi) xi yi f(xi,yi) xi yi f(xi,yi)
0.5 0.5 7.5 0.5 1. 6.75 0.5 1.5 5.5
1. 0.5 6.75 1. 1. 6. 1. 1.5 4.75
1.5 0.5 5.5 1.5 1. 4.75

Entonces

Rf(x,y)dA f(x1,y1)ΔA11 + f(x1,y2)ΔA12 + ... + f(x1,y3)ΔA13

V Q =Rf(x,y)dA 7.5 0.25 + 6.75 0.25 + 5.5 0.25 + 6.75 0.25 + 6 0.25 + 4.75 0.25 + 5.5 0.25 + 4.75 0.25

Por supuesto, esta aproximación no es muy buena. Se requiere una malla más fina para llegar cerca de V Q = 8π 4 21.1327.

Por ejemplo con Δx = Δy = 0.2, se obtiene V Q 16.2912 y con Δx = Δy = 0.05, se obtiene V Q 18.179

PIC

Fig. 6.8: Malla con Δx = Δy = 0.2

PIC

Fig. 6.9: Malla con Δx = Δy = 0.05

PIC

Fig. 6.10: Malla con Δx = Δy = 0.2

PIC

Fig. 6.11: Malla con Δx = Δy = 0.05

Las integrales dobles se puede aproximar con métodos especiales: Método de Trapecio, Método de Simpson, Método de Monte Carlo (Sección 6.), etc.