In multivariable calculus, an iterated integral is the result of applying integrals to a function of more than one variable (for example ${\displaystyle f(x,y)}$ or ${\displaystyle f(x,y,z)}$) in such a way that each of the integrals considers some of the variables as given constants. For example, the function ${\displaystyle f(x,y)}$, if ${\displaystyle y}$ is considered a given parameter, can be integrated with respect to ${\displaystyle x}$, ${\textstyle \int f(x,y)\,dx}$. The result is a function of ${\displaystyle y}$ and therefore its integral can be considered. If this is done, the result is the iterated integral

${\displaystyle \int \left(\int f(x,y)\,dx\right)\,dy.}$

It is key for the notion of iterated integrals that this is different, in principle, from the multiple integral

${\displaystyle \iint f(x,y)\,dx\,dy.}$

In general, although these two can be different, Fubini's theorem states that under specific conditions, they are equivalent.

The alternative notation for iterated integrals

${\displaystyle \int dy\int dx\,f(x,y)}$

is also used.

In the notation that uses parentheses, iterated integrals are computed following the operational order indicated by the parentheses starting from the most inner integral outside. In the alternative notation, writing ${\textstyle \int dy\,\int dx\,f(x,y)}$, the innermost integrand is computed first.

Examples

A simple computation

For the iterated integral

${\displaystyle \int \left(\int (x+y)\,dx\right)\,dy}$

the integral

${\displaystyle \int (x+y)\,dx={\frac {x^{2)){2))+yx}$

is computed first and then the result is used to compute the integral with respect to y.

${\displaystyle \int \left({\frac {x^{2)){2))+yx\right)\,dy={\frac {yx^{2)){2))+{\frac {xy^{2)){2))}$

This example omits the constants of integration. After the first integration with respect to x, we would rigorously need to introduce a "constant" function of y. That is, If we were to differentiate this function with respect to x, any terms containing only y would vanish, leaving the original integrand. Similarly for the second integral, we would introduce a "constant" function of x, because we have integrated with respect to y. In this way, indefinite integration does not make very much sense for functions of several variables.

The order is important

The order in which the integrals are computed is important in iterated integrals, particularly when the integrand is not continuous on the domain of integration. Examples in which the different orders lead to different results are usually for complicated functions as the one that follows.

Define the sequence ${\displaystyle a_{0}=0<a_{1}<a_{2}<\cdots }$ such that ${\displaystyle a_{n}\to 1}$. Let ${\displaystyle g_{n))$ be a sequence of continuous functions not vanishing in the interval ${\displaystyle (a_{n},a_{n+1})}$ and zero elsewhere, such that ${\textstyle \int _{0}^{1}g_{n}=1}$ for every ${\displaystyle n}$. Define

${\displaystyle f(x,y)=\sum _{n=0}^{\infty }\left(g_{n}(x)-g_{n+1}(x)\right)g_{n}(y).}$

In the previous sum, at each specific ${\displaystyle (x,y)}$, at most one term is different from zero. For this function it happens that^[1]

${\displaystyle \int _{0}^{1}\left(\int _{0}^{1}f(x,y)\,dy\right)\,dx=\int _{0}^{a_{1))\left(\int _{0}^{a_{1))g_{0}(x)g_{0}(y)\,dy\right)\,dx=1\neq 0=\int _{0}^{1}0\,dy=\int _{0}^{1}\left(\int _{0}^{1}f(x,y)\,dx\right)\,dy}$