Part 10 - General regions (1)

Introduction

In this part we’ll cover how to compute the double integral over general regions.

Double integrals over general regions

Let f(x,y)f(x, y) be defined on the domain, DD, where DD is a bounded region. DD is a general region, meaning it can have any geometrical shape.

However, if we reduce our problem into something we know, a rectangular region, we know how to compute the double integral.

Therefore, let’s enclose DD, with a rectangle, let’s call this rectangle for RR.

We define a function, f~\tilde{f}, on RR as the following:

f~(x,y)={f(x,y)(x,y)D0(x,y)D\tilde{f}(x, y) = \begin{cases} f(x, y) & (x, y) \in D \newline 0 & (x, y) \notin D \end{cases}
Definition 1 (Double integral over a general region)

The double integral of, ff over DD is:

Df(x,y) dA=Rf~(x,y) dA\iint_D f(x, y)\ dA = \iint_R \tilde{f}(x, y)\ dA
Definition 2 (Equivalent definition)

Recall our original definition of the double integral, we can also divide the domain, DD, into smaller pieces.

Therefore, the double integral of, ff, over DD is also:

Df(x,y) dA=limΔx0Δy0j=1mi=1nf(xi\*,yj\*)(ΔA)ij\iint_{D} f(x, y)\ dA = \lim_{\substack{\Delta x \to 0 \newline \Delta y \to 0}} \sum_{j = 1}^{m} \sum_{i = 1}^{n} f(x_{i}^{\*}, y_{j}^{\*}) (\Delta A)_{ij}

Since everything outside DD, is 0, they contribute nothing to the overall volume. Therefore, the volume under f~\tilde{f} is the same as the volume under ff.

How to compute general double integrals

To compute double integrals over a general region, we firstly need to introduce a new concept.

Definition 3 (Type I region)

A domain, DD, is of type I, if it lies between two graphs of two continuous functions of x.

Meaning:

D={(x,y)axb,g1(x)yg2(x)}D = \{(x, y) | a \leq x \leq b, g_1(x) \leq y \leq g_2(x)\}

Let’s see how we can compute this now. Let DD be of type I. Let’s enclose this region within a rectangle as well:

f(x,y) dA=Rf~(x,y) dA=abcdf~(x,y) dy dx\iint f(x, y)\ dA = \iint_R \tilde{f}(x, y)\ dA = \int_a^b \int_c^d \tilde{f}(x, y)\ dy\ dx

If we fix xx, then f~\tilde{f} becomes:

f~(x,y)={0cyg1(x)f(x,y)g1(x)yg2(x)0g2(x)yd\tilde{f}(x, y) = \begin{cases} 0 & c \leq y \leq g_1(x) \newline f(x, y) & g_1(x) \leq y \leq g_2(x) \newline 0 & g_2(x) \leq y \leq d \end{cases}

Which means:

abg1(x)g2(x)f(x,y) dy dx\boxed{\int_a^b \int_{g_1(x)}^{g_2(x)} f(x, y)\ dy\ dx}
Proposition 1 (Integration over a type I region)

If DD is of type I, then:

Df(x,y) dA=abg1(x)g2(x)f(x,y) dy dx\iint_D f(x, y)\ dA = \int_a^b \int_{g_1(x)}^{g_2(x)} f(x, y)\ dy\ dx
Definition 4 (Type II region)

D is of type II, if it lies between two graphs of two continuous functions of y.

D={(x,y)h1(y)xh2(y),cyd}D = \{(x, y) | h_1(y) \leq x \leq h_2(y), c \leq y \leq d\}

We perform the same procedure as type I, instead of fixing xx, we fix yy.

Proposition 2 (Integration over a type II region)

If DD is of type II, then:

Df(x,y) dA=cdh1(y)h2(y)f(x,y) dx dy\iint_D f(x, y)\ dA = \int_c^d \int_{h_1(y)}^{h_2(y)} f(x, y)\ dx\ dy

Before we head into doing a few examples, it’s a good practice always draw where our DD lies. So we easily get a geometrical sense of how it lies.

Example 1

Find Dy dA\iint_D y\ dA, where DD is the region bounded by, y=2x2,y=1+x2y = 2x^2, y = 1 + x^2.

So, our region is bounded by two parabolas, we have two functions of xx, therefore, our region is of type I.

Let’s firstly find where these graphs intersect:

2x2=1+x2x2=1x=±12x^2 = 1 + x^2 \newline x^2 = 1 \newline x = \pm 1

Graphing these we also see that, y=1+x2y = 1 + x^2 is the upper graph, and the lower is y=2x2y = 2x^2.

112x21+x2y dy dx\int_{-1}^1 \int_{2x^2}^{1 + x^2} y\ dy\ dx1211y2y=2x2y=1+x2 dx\dfrac{1}{2} \int_{-1}^1 y^2 \bigg\rvert_{y = 2x^2}^{y = 1 + x^2}\ dx1211(1+x2)2(2x2)2 dx\dfrac{1}{2} \int_{-1}^1 (1 + x^2)^2 - (2x^2)^2\ dx12111+2x2+x44x4 dx\dfrac{1}{2} \int_{-1}^1 1 + 2x^2 + x^4 - 4x^4 \ dx12111+2x23x4 dx\dfrac{1}{2} \int_{-1}^1 1 + 2x^2 - 3x^4\ dx12[x+23x335x5x=1x=1]\dfrac{1}{2} \left[x + \dfrac{2}{3} x^3 - \dfrac{3}{5} x^5 \bigg\rvert_{x = -1}^{x = 1} \right]12[(1+2335)(123+35)]\dfrac{1}{2} \left[\left(1 + \dfrac{2}{3} - \dfrac{3}{5}\right) - \left(-1 - \dfrac{2}{3} + \dfrac{3}{5}\right) \right]12[1+2335+1+2335]\dfrac{1}{2} \left[1 + \dfrac{2}{3} - \dfrac{3}{5} + 1 + \dfrac{2}{3} - \dfrac{3}{5} \right]1+2335\boxed{1 + \dfrac{2}{3} - \dfrac{3}{5}}
Example 2

Find Dxy dA\iint_D xy\ dA, where DD is the region bounded by, y=x1,y2=x+1y = x -1, y^2 = x+ 1.

Now, this region is technically both of type I and II. However, one of the choices will make our lives easier.

Let’s rewrite our bounds as, x=y+1,x=y21x = y + 1, x = y^2 -1.

Let’s find where these intersect:

y+1=y21y2y2=0y1==1y2==2y + 1 = y^2 - 1 \newline y^2 -y -2 = 0 \newline y_1 = \ldots = -1 \newline y_2 = \ldots = 2 \newline

Which means:

12y21y+1xy dx dy\int_{-1}^{2} \int_{y^2 - 1}^{y + 1} xy\ dx\ dy12y2[x2]x=y21x=y+1 dy\int_{-1}^{2} \dfrac{y}{2} \left[x^2\right] \bigg\rvert_{x = y^2 - 1}^{x = y + 1}\ dy12y2[(y+1)2(y21)2] dy\int_{-1}^{2} \dfrac{y}{2} \left[(y + 1)^2 - (y^2 - 1)^2\right]\ dy12y2[y2+2y+1y4+2y21] dy\int_{-1}^{2} \dfrac{y}{2} \left[y^2 + 2y + 1 - y^4 + 2y^2 - 1 \right]\ dy12y2[3y2+2yy4] dy\int_{-1}^{2} \dfrac{y}{2} \left[3y^2 + 2y - y^4 \right]\ dy1232y3+y2y52 dy\int_{-1}^{2} \dfrac{3}{2} y^3 + y^2 - \dfrac{y^5}{2}\ dy38y4+y33y612y=1y=2==278\dfrac{3}{8} y^4 + \dfrac{y^3}{3} - \dfrac{y^6}{12} \bigg\rvert_{y = -1}^{y = 2} = \ldots = \boxed{\dfrac{27}{8}}

Changing order of integral

One might now think that we can change the order of integration, using Fubini’s Theorem.

That is not the case.

Say we have:

g1(x)g2(x)abf(x,y) dx dy\int_{g_1(x)}^{g_2(x)} \int_a^b f(x, y)\ dx\ dy

Computing this entire integral will not yield a singular number, it will yield a function of xx instead.

However, if DD is both type I and II, we can change order of integration, but we need to change our limits accordingly.

Example 3

Find Dsin(y2) dA\iint_D sin(y^2)\ dA, where DD is the region bounded by, x=0,y=1,y=xx = 0, y = 1, y = x

Graphing this we can see that this a region of both types.

For, type I:

Df(x,y) dA=01xysin(y2) dy dx\iint_D f(x, y)\ dA = \int_0^1 \int_x^y sin(y^2)\ dy\ dx

Type II:

Df(x,y) dA=010ysin(y2) dx dy\iint_D f(x, y)\ dA = \int_0^1 \int_0^y sin(y^2)\ dx\ dy

In this case, one of these isn’t even integrable in terms of elementary functions. So yeah…

Solving the second integral:

010ysin(y2) dx dy\int_0^1 \int_0^y sin(y^2)\ dx\ dy01xsin(y2)x=0x=y dy\int_0^1 x \cdot sin(y^2) \bigg\rvert_{x = 0}^{x = y}\ dy01ysin(y2) dy\int_0^1 y \cdot sin(y^2) \ dy

Using normal u-sub, we obtain:

12cos(y2)y=0y=1-\dfrac{1}{2} cos(y^2) \bigg\rvert_{y = 0}^{y = 1}(12cos(1))(12cos(0))\left( -\dfrac{1}{2} cos(1)\right) - \left( -\dfrac{1}{2} cos(0)\right)12cos(1)+12-\dfrac{1}{2} cos(1) + \dfrac{1}{2}12(cos(1)+1)\boxed{\dfrac{1}{2} \left(-cos(1) + 1 \right)}