Definition of Double Integrals

Section 5.1 Definition of Double Integrals

Before we introduce double integrals let us recall how the integral of a function of one variable, say

t,

was defined. Geometrically, the definite integral of a function

f

over an interval

I = [a, b]

represents the area of the region between the graph and the

t

-axis. To find an approximation of this area the idea is to use a partition of

I

given by

a = t_{0} < t_{1} < t_{2} < \dots < t_{n} = b .

Then we choose a point

t_{i}^{*}

in every interval

[t_{i - 1}, t_{i}] .

Now we look at the rectangle of width

Δ t_{i} := t_{i} - t_{i - 1}

and height

f (t_{i}^{*}) .

Such a rectangle is highlighted in Figure 5.1. The area of this rectangle is

f (t_{i}^{*}) Δ t_{i} .

Hence the total area of all rectangles as shown in Figure 5.1 is

🔗

\sum_{i = 1}^{n} f (t_{i}^{*}) Δ t_{i} .

This total area is an approximation for the area under the graph lightly shaded in Figure 5.1. The above sum is called a Riemann sum of

f

over

I .

Taking a finer partition, that is,

max Δ t_{i}

smaller, we expect that the corresponding Riemann sum provides a better approximation. Indeed, for continuous functions the limit

🔗

\int_{a}^{b} f (t) d t := lim_{max Δ t_{i} \to 0} \sum_{i = 1}^{n} f (t_{i}^{*}) Δ t_{i}

exists, and is independent of the intermediate points

t_{i}^{*}

chosen. The limit is called the Riemann integral or simply the integral of

f

over the interval

[a, b] .

f (t_{i}^{*})

is negative then the area is considered negative.

🔗

Figure 5.1. Approximation of the area under the graph of a function.
🔗

The idea to define integrals of functions of two variables (double integrals) is very similar to the above. The area under the graph of a function of one variable corresponds to the volume of the region between the graph of a function of two variables and the

x y

-plane. To get an approximation of that volume we replace the rectangles by rectangular boxes, and then pass to the limit by taking smaller and smaller rectangular boxes. Two such approximations are shown in Figure 5.2, together with the limiting case.

🔗

Figure 5.2. Approximation of the volume under the graph by rectangular boxes
🔗

To be more precise let us look at a function,

f,

defined on a rectangle

[a, b] \times [c, d]

of the plane. By choosing partitions

a = x_{0} < x_{1} < \dots < x_{m} = b

and

c = y_{0} < y_{2} < \dots < y_{n} = d

of the intervals

[a, b]

and

[c, d]

we get a subdivision of the rectangle

[a, b] \times [c, d]

into small rectangles

[x_{i - 1}, x_{i}] \times [y_{j - 1}, y_{j}],

one of which is highlighted in Figure 5.3.

🔗

The lengths of the edges of such a rectangle is

Δ x_{i}

and

Δ y_{j},

respectively. Now we choose an arbitrary point

x_{i j}^{*} := (x_{i}^{*}, y_{j}^{*})

in that rectangle similarly as done in the one dimensional case with

t_{i}^{*} .

Then we look at the rectangular box with basis

[x_{i - 1}, x_{i}] \times [y_{j - 1}, y_{j}]

and height

f (x_{i}^{*}, y_{j}^{*}) .

Such a box is shown in Figure 5.4.

🔗

Figure 5.4. The graph of a function and the box over a small rectangle.
🔗

Its volume is the product of the lengths of its edges, namely

f (x_{i}^{*}, y_{j}^{*}) Δ x_{i} Δ y_{j} .

f (x_{i}^{*}, y_{j}^{*})

is negative then the volume is considered negative, as in the case of a function of one variable. Hence the total volume of all boxes is

🔗

\begin{matrix} (5.1) & \sum_{i = 1}^{m} \sum_{j = 1}^{n} f (x_{i}^{*}, y_{j}^{*}) Δ x_{i} Δ y_{j} . \end{matrix}

This sum corresponds to the Riemann sums we looked at before, and is also called a Riemann sum. It turns out that, if

f

is continuous on the closed rectangle, then one can pass to the limit as

max Δ x_{i} \to 0

and

max Δ y_{j} \to 0 .

This leads to the following definition.

🔗

Definition 5.5.

Suppose that

f

is continuous on the rectangle

R = [a, b] \times [c, d] .

Then we call

🔗

\begin{aligned} \int_{R} f (x) d x & := \iint_{R} f (x, y) d x d y \\ := lim_{\begin{array}{c} max Δ x_{i} \to 0 \\ max Δ y_{j} \to 0 \end{array}} \sum_{i = 1}^{m} \sum_{j = 1}^{n} f (x_{i}^{*}, y_{j}^{*}) Δ x_{i} Δ y_{j} \end{aligned}

the double integral or simply the integral of

f

over

R

🔗

Remark 5.6.

The notation for integrals used in the literature varies. In older or more applied texts the notation

🔗

\iint_{R} f (x, y) d x d y

is used for double integrals. In the more modern literature one simply writes

🔗

\int_{R} f (x) d x .

In this course we will mainly use the second, more concise form.

🔗

We would like to integrate over more general regions than just rectangles. We can certainly integrate over a union of finitely many rectangles. If we want to integrate over more general sets then we need to approximate the sets by a union of small rectangles as well. To get upper and lower bounds we approximate from the inside and the outside as shown in Figure 5.7.

🔗

Figure 5.7. Approximation of a domain from the in- and outside by rectangles.
🔗

We then form a sum similar to (5.1) over the relevant rectangles. If the set is reasonably `nice,’ that is, its boundary is not too wild, then the limits of the Riemann sums from the inside and the outside are the same, and we can define the double integral over such regions. It is quite difficult to characterise domains for which this is possible, and we omit such a discussion. All domains we encounter in this course have the required properties.

🔗

We next collect some properties of the integral. They follow directly from the definition as limits of Riemann sums.

🔗

Theorem 5.8.

Suppose that

D, D_{1}, D_{2} \subset R^{2}

are domains on which the integral of every continuous function is well defined. If

f, g

are continuous functions on

D

and

α \in R

then

🔗

$\int_{D} (f (x) + g (x)) d x = \int_{D} f (x) d x + \int_{D} g (x) d x$
$\int_{D} α f (x) d x = α \int_{D} f (x) d x$
If $D = D_{1} \cup D_{2}$ and $D_{1} \cap D_{2} = \emptyset$ then

$\int_{D} f (x) d x = \int_{D_{1}} f (x) d x + \int_{D_{2}} f (x) d x .$

🔗

Remark 5.9.

We introduced the double integral as the volume of the region between the graph of a function and the

x y

-plane. Note, however, that this is not the only possible interpretation. If we choose

f

to be the constant function with value one, that is,

f (x, y) = 1,

then the integral represents the area of the domain:

🔗

\begin{matrix} (5.2) & Area of the domain D = \int_{D} 1 d x . \end{matrix}

🔗

Another physical interpretation is the following. Let

D

represent the area occupied by a thin plate made of some material. At every point

(x, y)

on that plate let

ϱ (x, y)

denote the mass density of the material. Then the

🔗

Total mass of the plate occupying D = \int_{D} ϱ (x) d x .

🔗

The theory developed in this section does not allow us to compute the double integral of a given function. We will see how to do this in the next section.

🔗

Prev Top Next