This interpretation of double integrals is basically the same. Instead of area, we have volume:
∬Df(x,y)dA=volume under f(x,y) over D
You can already see that unlike in 1D, the region you integrate over can get very complicated. In 1D, you mostly integrated over intervals, but in 2D, D can be a rectangle, a disk, a star, etc. However, despite this, Fubini's theorem tells you that a double integral can be (in a nutshell) calculated as two 1D integrals:
Theorem (Fubini)
∬Df(x,y)dA=∬Df(x,y)dxdy=∬Df(x,y)dydx
In other words, Fubini's theorem tells you that you can calculate a double integral as an iterated integral, and that the order of the variables you integrate in doesn't matter.
Example 1.
Integrate f(x,y)=(x+y)2 over R=[0,2]×[−1,1].
Solution.
When given a rectangle, the first interval always corresponds to x and the second interval always corresponds to y. This gives us the following bounds:
Integrate f(x,y)=ex over the region D bounded by y=x+1, y=x, x=0, and x=1.
Solution.
Unlike the previous example, D is slightly more complicated than a rectangle. For example, for each x in the interval [0,1], the interval that y lives in changes:
x=0x=1⟹0≤y≤1⟹1≤y≤2.
This means that the bounds of y are a function of x:
For each x in [0,1], y will range over [x,x+1], so the integral is