← Back to calculus 3

Double Integrals over General Domains

Double Integrals over General Domains

In the previous note we integrated over rectangles. Most regions of interest in applications are not rectangular. This note explains how to integrate over regions bounded by curves.


1. Type I regions

A Type I region DD is bounded on the left and right by x=ax = a and x=bx = b, and vertically by two continuous curves:

D={(x,y):axb,    g1(x)yg2(x)},D = \{\,(x,y) : a \le x \le b,\;\; g_1(x) \le y \le g_2(x)\,\},

where g1(x)g2(x)g_1(x) \le g_2(x) for all x[a,b]x \in [a,b].

Think of a Type I region as one that can be sliced into vertical strips. For each fixed xx, the strip runs from y=g1(x)y = g_1(x) (the lower curve) up to y=g2(x)y = g_2(x) (the upper curve).

The interactive graph below shows the Type I region between y=x2y = x^2 and y=2xy = 2x on [0,2][0, 2]. Adjust the slider to see the vertical strips that make up the region.

Type I: for each fixed x, y ranges from g₁(x) to g₂(x). Each vertical strip is one inner integral.


2. Iterated integral over a Type I region

The double integral over a Type I region is evaluated by integrating yy first (over the strip), then integrating xx:

Df(x,y)dA=ab ⁣g1(x)g2(x)f(x,y)dydx.\boxed{\iint_D f(x,y)\,dA = \int_a^b \!\int_{g_1(x)}^{g_2(x)} f(x,y)\,dy\,dx.}

How to read this: for each fixed xx, compute the inner integral A(x)=g1(x)g2(x)f(x,y)dy\displaystyle A(x) = \int_{g_1(x)}^{g_2(x)} f(x,y)\,dy (one slice of volume), then integrate A(x)A(x) over x[a,b]x \in [a,b].

Key point: The limits of the inner integral depend on xx, so they vary as you move the outer integration variable.


3. Type II regions

A Type II region DD is bounded above and below by y=cy = c and y=dy = d, and horizontally by two continuous curves:

D={(x,y):cyd,    h1(y)xh2(y)},D = \{\,(x,y) : c \le y \le d,\;\; h_1(y) \le x \le h_2(y)\,\},

where h1(y)h2(y)h_1(y) \le h_2(y) for all y[c,d]y \in [c,d].

A Type II region can be sliced into horizontal strips. For each fixed yy, the strip runs from x=h1(y)x = h_1(y) (the left curve) to x=h2(y)x = h_2(y) (the right curve).

Type II: for each fixed y, x ranges from h₁(y) to h₂(y). Each horizontal strip is one inner integral.

The double integral over a Type II region is

Df(x,y)dA=cd ⁣h1(y)h2(y)f(x,y)dxdy.\boxed{\iint_D f(x,y)\,dA = \int_c^d \!\int_{h_1(y)}^{h_2(y)} f(x,y)\,dx\,dy.}


4. Example: Area between y=2x2y = 2x^2 and y=1+x2y = 1 + x^2

Find the area of the region DD between the parabolas y=2x2y = 2x^2 and y=1+x2y = 1 + x^2.

Step 1: Find the intersection.

2x2=1+x2    x2=1    x=±1.2x^2 = 1 + x^2 \implies x^2 = 1 \implies x = \pm 1.

At x=0x = 0: 2x2=02x^2 = 0 and 1+x2=11+x^2 = 1, so 2x21+x22x^2 \le 1+x^2 on [1,1][-1,1]. The lower curve is y=2x2y = 2x^2, the upper curve is y=1+x2y = 1 + x^2.

Step 2: Set up as Type I.

Area=DdA=11 ⁣2x21+x2dydx.\text{Area} = \iint_D dA = \int_{-1}^{1}\!\int_{2x^2}^{1+x^2} dy\,dx.

Step 3: Evaluate the inner integral.

2x21+x2dy=(1+x2)2x2=1x2.\int_{2x^2}^{1+x^2} dy = (1+x^2) - 2x^2 = 1 - x^2.

Step 4: Evaluate the outer integral.

11(1x2)dx=[xx33]11=(113)(1+13)=43.\int_{-1}^{1}(1 - x^2)\,dx = \Bigl[x - \tfrac{x^3}{3}\Bigr]_{-1}^{1} = \Bigl(1 - \tfrac{1}{3}\Bigr) - \Bigl(-1 + \tfrac{1}{3}\Bigr) = \frac{4}{3}.

The area is 43\dfrac{4}{3}.

For a non-trivial integrand, say f(x,y)=x+yf(x,y) = x + y, the same limits apply:

D(x+y)dA=11 ⁣2x21+x2(x+y)dydx.\iint_D (x+y)\,dA = \int_{-1}^{1}\!\int_{2x^2}^{1+x^2}(x+y)\,dy\,dx.

Inner integral:

2x21+x2(x+y)dy=[xy+y22]2x21+x2=x(1x2)+(1+x2)2(2x2)22.\int_{2x^2}^{1+x^2}(x+y)\,dy = \bigl[xy + \tfrac{y^2}{2}\bigr]_{2x^2}^{1+x^2} = x(1-x^2) + \tfrac{(1+x^2)^2 - (2x^2)^2}{2}.

=x(1x2)+1+2x23x42.= x(1-x^2) + \tfrac{1+2x^2-3x^4}{2}.

Since x(1x2)x(1-x^2) is an odd function on [1,1][-1,1], its integral is 00:

111+2x23x42dx=[x2+x333x510]11=2(12+13310)=1615.\int_{-1}^{1} \tfrac{1+2x^2-3x^4}{2}\,dx = \Bigl[\tfrac{x}{2}+\tfrac{x^3}{3}-\tfrac{3x^5}{10}\Bigr]_{-1}^{1} = 2\Bigl(\tfrac{1}{2}+\tfrac{1}{3}-\tfrac{3}{10}\Bigr) = \frac{16}{15}.


5. Switching order: Type I vs Type II

Sometimes it is easier to describe a region as Type II and integrate in the order dxdydx\,dy, even if the natural first description is Type I.

Example. Evaluate 01 ⁣x1sin(y3)dydx\displaystyle\int_0^1\!\int_{\sqrt{x}}^{1} \sin(y^3)\,dy\,dx.

The inner integrand sin(y3)\sin(y^3) has no elementary antiderivative with respect to yy, so the dydy-first approach fails.

The region is {(x,y):0x1,xy1}\{(x,y): 0 \le x \le 1,\,\sqrt{x} \le y \le 1\}, which equals {(x,y):0y1,0xy2}\{(x,y): 0 \le y \le 1,\, 0 \le x \le y^2\} as a Type II region.

Switching order:

01 ⁣0y2sin(y3)dxdy=01y2sin(y3)dy=[cos(y3)3]01=1cos130.153.\int_0^1\!\int_0^{y^2} \sin(y^3)\,dx\,dy = \int_0^1 y^2 \sin(y^3)\,dy = \Bigl[-\tfrac{\cos(y^3)}{3}\Bigr]_0^1 = \frac{1-\cos 1}{3} \approx 0.153.


6. Properties of double integrals

Let DD be a region in the plane and ff, gg be integrable functions on DD.

Linearity:

D(c1f+c2g)dA=c1DfdA+c2DgdA.\iint_D \bigl(c_1 f + c_2 g\bigr)\,dA = c_1\iint_D f\,dA + c_2\iint_D g\,dA.

Additivity: If D=D1D2D = D_1 \cup D_2 with D1D2D_1 \cap D_2 having zero area, then

DfdA=D1fdA+D2fdA.\iint_D f\,dA = \iint_{D_1} f\,dA + \iint_{D_2} f\,dA.

Comparison: If f(x,y)g(x,y)f(x,y) \le g(x,y) on DD, then

DfdADgdA.\iint_D f\,dA \le \iint_D g\,dA.

Area as a special case: Setting f=1f = 1,

D1dA=Area(D).\iint_D 1\,dA = \text{Area}(D).

Bound: If mf(x,y)Mm \le f(x,y) \le M on DD, then

mArea(D)DfdAMArea(D).m \cdot \text{Area}(D) \le \iint_D f\,dA \le M \cdot \text{Area}(D).


7. Exercises

Evaluate DxdA\displaystyle\iint_D x\,dA where DD is bounded by y=xy = x and y=x2y = x^2 (intersecting at (0,0)(0,0) and (1,1)(1,1)).

Solution.

Type I: 0x10 \le x \le 1, x2yxx^2 \le y \le x.

01 ⁣x2xxdydx=01x(xx2)dx=01(x2x3)dx=[x33x44]01=112.\int_0^1\!\int_{x^2}^{x} x\,dy\,dx = \int_0^1 x(x - x^2)\,dx = \int_0^1 (x^2 - x^3)\,dx = \Bigl[\tfrac{x^3}{3} - \tfrac{x^4}{4}\Bigr]_0^1 = \frac{1}{12}.

Evaluate 04 ⁣y2ex3dxdy\displaystyle\int_0^4\!\int_{\sqrt{y}}^{2} e^{x^3}\,dx\,dy.

Solution.

The region {(x,y):0y4,yx2}\{(x,y): 0 \le y \le 4,\,\sqrt{y} \le x \le 2\} equals {(x,y):0x2,0yx2}\{(x,y): 0 \le x \le 2,\, 0 \le y \le x^2\} as a Type I region.

02 ⁣0x2ex3dydx=02x2ex3dx=[ex33]02=e813.\int_0^2\!\int_0^{x^2} e^{x^3}\,dy\,dx = \int_0^2 x^2 e^{x^3}\,dx = \Bigl[\tfrac{e^{x^3}}{3}\Bigr]_0^2 = \frac{e^8 - 1}{3}.

Find the area of the region DD enclosed by y=x3y = x^3 and y=xy = x in the first quadrant.

Solution.

They intersect at (0,0)(0,0) and (1,1)(1,1). For 0x10 \le x \le 1: x3xx^3 \le x, so DD is Type I with x3yxx^3 \le y \le x.

Area=01(xx3)dx=[x22x44]01=1214=14.\text{Area} = \int_0^1 (x - x^3)\,dx = \Bigl[\tfrac{x^2}{2} - \tfrac{x^4}{4}\Bigr]_0^1 = \frac{1}{2} - \frac{1}{4} = \frac{1}{4}.

Evaluate DydA\displaystyle\iint_D y\,dA where D={(x,y):0y1,y2xy}D = \{(x,y): 0 \le y \le 1,\, y^2 \le x \le y\}.

Solution.

Type II with h1(y)=y2h_1(y) = y^2, h2(y)=yh_2(y) = y:

01 ⁣y2yydxdy=01y(yy2)dy=01(y2y3)dy=1314=112.\int_0^1\!\int_{y^2}^{y} y\,dx\,dy = \int_0^1 y(y - y^2)\,dy = \int_0^1 (y^2 - y^3)\,dy = \frac{1}{3} - \frac{1}{4} = \frac{1}{12}.

Without computing, bound D(1+sin2x)dA\displaystyle\iint_D (1 + \sin^2 x)\,dA where D={(x,y):0x1,0yx}D = \{(x,y): 0 \le x \le 1,\, 0 \le y \le x\}.

Solution.

Since 0sin2x10 \le \sin^2 x \le 1, we have 11+sin2x21 \le 1+\sin^2 x \le 2 on DD.

The area of DD (the triangle 0yx10\le y\le x\le 1) is 12\frac{1}{2}.

By the bounding property:

12D(1+sin2x)dA1.\frac{1}{2} \le \iint_D(1+\sin^2 x)\,dA \le 1.