El teorema de Fubini permite calcular la integral doble de una función continua de dos variables sobre un rectángulo mediante dos integrales sucesivas sobre un intervalo, no importando el orden de integración, lo [br]que se suele llamar como integrales iteradas.[br][br]Aún cuando el dominio no es rectangular, una integral doble se puede calcular mediante integración iterada, aunque en este caso los límites de las integrales interiores no tienen por qué ser constantes. Vamos a ver un ejemplo donde el dominio es triangular. Consideremos un dominio formado por un triángulo, [math]T[/math] , de vértices [math]A=\left(a,c\right)[/math] , [math]B=\left(b,c\right)[/math] y [math]\text{C=(b,d)}[/math]. Los lados del triángulo se escriben en términos de las variables [math]x[/math] e [math]y[/math] (variables de integración) como:[list][*] Segmento [math]AB[/math], corresponde a [math]y=c[/math] (constante) y [math]a\le x\le b[/math].[/*][*]Segmento [math]BC[/math], corresponde a [math]x=b[/math] (constante) y [math]c\le y\le d[/math].[/*][*]Segmento [math]AC[/math], está en la recta [math]\left(d-c\right)\cdot\left(x-a\right)-\left(b-a\right)\left(y-c\right)=0[/math], para valores de [math]x[/math] con [math]a\le x\le b[/math] (si queremos escribir [math]y[/math] en función de [math]x[/math]) o si la queremos escribir [math]x[/math] en términos de [math]y[/math], con valores de [math]y[/math] con [math]c\le y\le d[/math].[/*][/list][list][/list]Para calcular la integral de una función [math]f[/math] sobre el triángulo [math]T[/math], se puede [i]barrer[/i] el triángulo mediante segmentos horizontales de longitud variable. La longitud de cada segmento, depende de la altura [math]y[/math] a la que se encuentre. Por la elección de triángulo que hemos hecho, cada segmento horizontal a altura [math]y_0[/math] se apoya en el segmento [math]AC[/math] (límite inferior) y en segmento [math]BC[/math] (límite superior). En términos de las variables, se tiene [math]y=y_0[/math] (constante) y [math]x[/math] varía entre [math]\psi\left(y_0\right)=a+\frac{b-a}{d-c}\left(y_0-c\right)[/math] y [math]b[/math] (la función [math]\psi\left(y_{ }\right)[/math] se obtiene despejando [math]x[/math] de la recta que contiene al segmento [math]AC[/math]. El triángulo se barre complentamente variando la altura [math]y[/math] de [math]c[/math] a [math]d[/math]. Es decir[br][br][math]\displaystyle\int\int_{T}f(x,y)\,dxdy=\displaystyle\int_c^d\left(\int_{a+\frac{b-a}{d-c}(y-c)}^{b}f(x,y)\,dx\right)dy[/math].[br][br]La integral interior, cuando [math]y=y_0[/math], es el área con signo bajo la curva [math]f\left(x,y_0\right)[/math] al variar [math]x[/math] entre [math]\psi\left(y_0\right)[/math] y [math]b[/math].[br][br]De manera similar, se puede barrer el triángulo [math]T[/math] con segmentos verticales de longitud variable. Estos segmentos se apoyan inferiormente en el lado del triángulo [math]AB[/math] y superiormente en el lado [math]AC[/math]. En términos de las variables de integración en este caso [math]x=x_0[/math] (constante) y la variable [math]y[/math] está entre [math]c[/math] y [math]\varphi\left(x_0\right)=c+\frac{d-c}{b-a}\left(x_0-a\right)[/math], donde la función [math]\varphi\left(x\right)[/math] se encuentra despejando la variable [math]y[/math] en la ecuación de la recta que contiene al segmento [math]AC[/math]. Para barrer el triángulo completamente la [math]x[/math] tiene que variar entre [math]a[/math] y [math]b[/math]. En esta situación la integral doble se escribe[br][br][br][math]\displaystyle\int\int_Tf(x,y)\,dxdy=\displaystyle\int_{a}^{b}\left(\int_{c}^{c+\frac{d-c}{b-a}(x-a)}f(x,y)\,dy\right)dx[/math].[br][br]En este caso, la integral interior es el área con signo bajo la curva [math]f\left(x_{0,}y\right)[/math], al variar [math]y[/math] entre [math]c[/math] y [math]\varphi\left(x_0\right)[/math].[br][br]Obsérvese que tanto en un caso como en el otro la integral exterior es la que tiene los límites de integración constantes y los límites de la integral interior depende de la variable que [b]no[/b] es la variable de integración de esa integral.

Arriba a la derecha se ve el triángulo, que es el dominio de integración. Abajo en azul está la gráfica [math]f[/math] de definida en ese dominio.[br][br]Arriba a la izquierda, al marcar la casilla "Integración primero con respecto a x" aparece en el valor de [math]y[/math]indicado por el deslizador verde, a la derecha se ve sobre el triángulo el segmento horizontal rojo a esa altura y en la parte de abajo de la construcción se aprecia la gráfica de la función (con [math]y=y_0[/math] constante y [math]x[/math] entre [math]\psi\left(y_0\right)=a+\frac{b-a}{d-c}\left(y_0-c\right)[/math] y [math]b[/math]). La zona coloreada en rojo que aparece bajo la gráfica de [math]f[/math] representa la integral interior y el plano azul es el plano [math]y=y_0[/math]. Al variar [math]y_0[/math] con el deslizador va variando la zona coloreada y su integral. [br][br]De forma análoga, al marcar la casilla "Integración primero con respecto a y" aparece en el valor de [math]x[/math]indicado por el deslizador rojo, a la derecha el segmento vertical en [math]x=x_0[/math] del dominio, y abajo la gráfica de la función (con [math]x=x_0[/math] constante e [math]y[/math] entre [math]c[/math] y [math]\varphi\left(x_0\right)[/math]). La zona coloreada en verde bajo la gráfica de [math]f[/math] representa la integral interior y el plano azul es ahora el plano [math]x=x_0[/math] . Al variar [math]x[/math] con el deslizador va variando la zona coloreada y su integral. [br][br]La función [math]f[/math], y los valores de [math]a[/math], [math]b[/math], [math]c[/math] y [math]d[/math] se pueden introducir en las casillas de entrada.