Skip to main content

Section 5.3 Integrating a two-form over a region in \(\mathbb{R}^2\)

Step 2: we define the integral of a two-form on \(U \subset \mathbb{R}^2\) along a closed bounded domain in \(U\) with a choice of orientation. The definition is in terms of standard double integrals from calculus. We show that the resulting integral is invariant under orientation-preserving reparametrizations, and changes sign under orientation-reversing reparametrizations. This is closely connected to the transformation formula for double integrals.

Subsection 5.3.1 The integral of a two-form over an oriented closed bounded region in \(\mathbb{R}^2\)

We define the integral of a two-form over an oriented closed bounded region in \(\mathbb{R}^2\text{.}\)

Definition 5.3.1. Integral of a two-form over an oriented closed bounded region in \(\mathbb{R}^2\).

Let \(D \subset \mathbb{R}^2\) be a closed bounded region. Let \((x,y)\) be coordinates on \(\mathbb{R}^2\text{.}\) The canonical (\(+\)) orientation on \(D\) is described by the ordered basis \(\{ \mathbf{e}_1, \mathbf{e}_2 \}\) for \(\mathbb{R}^2\text{,}\) with \(\mathbf{e}_1=(1,0)\) and \(\mathbf{e}_2=(0,1)\) the unit vectors pointing in the positive \(x\)- and \(y\)-directions. Let \(\omega = f\ dx \wedge dy\) be a two-form on an open subset \(U \subseteq \mathbb{R}^2\) such that \(D \subset U\text{.}\) Then, we define the integral of \(\omega\) over \(D\) with canonical orientation as:

\begin{equation*} \int_{D_+} \omega = \iint_D f\ dA, \end{equation*}

where on the right-hand-side we mean the standard double integral from calculus of the function \(f\) over the region \(D\text{.}\) If \(D\) is given the opposite orientation, we define the integral of \(\omega\) over \(D_-\) as:

\begin{equation*} \int_{D_-} \omega = - \iint_D f\ dA. \end{equation*}

Remark 5.3.2.

A bit of notation: we will always use double or triple integral signs to denote the standard double and triple integrals from calculus, while we will use only one integral sign when we are integrating a differential form.

It is probably worth recalling here how double integrals are defined, from your previous calculus course. If \(D\) is a rectangular domain, that is,

\begin{equation*} D = \{ (x,y) \in \mathbb{R}^2 ~|~ x \in [a,b], y \in [c,d] \}, \end{equation*}

then the double integral is defined as an iterated integral:

\begin{equation*} \iint_D f\ dA = \int_c^d \int_a^b f\ dx dy. \end{equation*}

The notation here means that the inner integral is an integral with respect to \(x\) (while keeping \(y\) fixed), while the outer integral is with respect to \(y\text{.}\) 1  We recall Fubini's theorem, which states that the order of integration does not matter:

\begin{equation*} \int_c^d \int_a^b f\ dx dy = \int_a^b \int_c^d f\ dy dx, \end{equation*}

that is, it does not matter whether you integrate in \(x\) first and then in \(y\) or the other way around.

Note that on the right-hand-side here there is no wedge product: this is not the integral of a two-form, it is an interated integral in \(x\) and \(y\) as you have seen in calculus. The inner integral is with respect to \(x\) (keeping \(y\) fixed), while the outer integral is with respect to \(y\text{.}\)

To integrate over a more general closed bounded region \(D\text{,}\) we proceed as follows. Since \(D\) is bounded, we can take it to be inside a rectangular region. We can then extend the function \(f\) to the rectangular region by setting it to zero everywhere outside \(D\text{.}\) The double integral over \(D\) is then defined to be the integral of the extended function over the rectangular region, which can be written as an iterated itegral as above.

There are two types of regions that give rise to nice iterated integrals. If \(D\) can be written as follows:

\begin{equation*} D = \{(x,y) \in \mathbb{R}^2~|~ x \in [a,b], u(x) \leq y \leq v(x) \}, \end{equation*}

with \(u,v : \mathbb{R} \to \mathbb{R}\) continuous functions, we say that the region is \(x\)-supported (or of type I ). In this case, one can show that the double integral can be written as the following iterated integral:

\begin{equation*} \iint_D f\ dA = \int_a^b \int_{u(x)}^{v(x)} f\ dy dx, \end{equation*}

where the inner integral is with respect to \(x\) (keeping \(y\) fixed), while the outer integral is with respect to \(y\text{.}\)

If instead \(D\) can be written as:

\begin{equation*} D = \{(x,y) \in \mathbb{R}^2~|~ y \in [c,d], u(y) \leq x \leq v(y) \}, \end{equation*}

we say that \(D\) is \(y\)-supported (or of type II). The double integral is then the iterated integral:

\begin{equation*} \iint_D f\ dA = \int_c^d \int_{u(y)}^{v(y)} f\ dx dy, \end{equation*}

with the inner integral being with respect to \(y\) (keeping \(x\) fixed), and the outer integral with respect to \(y\text{.}\)

Note that rectangular regions are particular cases of both \(x\)-supported and \(y\)-supported regions. If a region \(D\) is either \(x\)-supported or \(y\)-supported, we say that it is recursively supported. Most of the regions that we will deal with will be either recursively supported regions, or regions that can be expressed as unions of recursively supported regions.

Consider the two-form \(\omega =x y\ dx \wedge dy\) on \(\mathbb{R}^2\text{,}\) and the closed bounded region

\begin{equation*} D = \{(x,y)\in\mathbb{R}^2~|~0 \leq x \leq 2, -1 \leq y \leq 3 \}, \end{equation*}

equipped with the canonical orientation. Then

\begin{align*} \int_{D_+} \omega =\amp \iint_D x y\ dA\\ =\amp \int_{-1}^3 \int_0^2 x y\ dx dy\\ =\amp \int_{-1}^3 y \left[ \frac{x^2}{2} \right]_{x=0}^{x=2} \ dy\\ =\amp 2 \int_{-1}^3 y\ dy\\ =\amp y^2 \Big|_{-1}^3\\ =\amp 8. \end{align*}

To evaluate the double integral, we used the standard procedure for evaluating double integrals over rectangular regions as iterated integrals, with the inner integral being with respect to \(x\) (keeping \(y\) constant), and the outer integral being with respect to \(y\text{.}\)

Consider the two-form \(\omega = x e^{y}\ dx \wedge dy\text{,}\) and the region

\begin{equation*} D = \{(x,y) \in \mathbb{R}^2~|~x \in [1,2], \ln(x) \leq y \leq \ln(2 x)\}. \end{equation*}

\(D\) is an \(x\)-supported region. The integral then reads:

\begin{align*} \int_{D_+} \omega =\amp \int_1^2 \int_{\ln(x)}^{\ln(2 x)} x e^{y}\ dy dx\\ =\amp \int_1^2 x \left[e^y \right]_{\ln(x)}^{\ln(2 x)} dx\\ =\amp \int_1^2 \left(x e^{\ln(2 x)} - x e^{\ln(x)} \right) dx\\ =\amp \int_1^2 \left( 2 x^2 - x^2 \right) dx\\ =\amp \int_1^2 x^2\ dx\\ =\amp \frac{8}{3} - \frac{1}{3}\\ =\amp \frac{7}{3}. \end{align*}

Remark 5.3.5.

It is worth pointing out here that Definition 5.3.1 is actually quite subtle. We defined the integral of the two-form \(\omega = f \ dx \wedge dy\) over \(D_+\) as

\begin{equation*} \int_{D_+} f \ dx \wedge dy = \iint_D f\ dA, \end{equation*}

where the right-hand-side is the standard double integral of a function in calculus. The subtelty is that the right-hand-side is not an oriented integral, while the left-hand-side is. Indeed, suppose for simplicity that \(D = [a,b] \times [c,d]\) is a rectangular region, as in Example 5.3.3. Then, we can interpret the right-hand-side as an iterated integral:

\begin{equation*} \iint_D f\ dA = \int_c^d \int_a^b f\ dx dy. \end{equation*}

But then, by Fubini's theorem, we know that we can exchange the order of integration without issue. That is,

\begin{equation*} \iint_D f\ dA = \int_c^d \int_a^b f\ dx dy = \int_a^b \int_c^d f\ dy dx. \end{equation*}

At first sight, this may appear problematic, as one could be tempted to reinterpret the last integral as

\begin{equation*} \int_{D_+} f\ dy \wedge dx, \end{equation*}

but since \(dx \wedge dy = - dy \wedge dx\text{,}\) this would be minus the integral we started with! But this is incorrect. The subtlety is in the choice of orientation.

The key is that in Definition 5.3.1, we started by choosing coordinates \((x,y)\) on \(\mathbb{R}^2\text{,}\) and then we wrote the one-form \(\omega = f\ dx \wedge dy\) using the basic two-form \(dx \wedge dy\) in which the differentials \(dx\) and \(dy\) appear in the same order as the coordinates \((x,y)\text{.}\) This is important, as \(dy \wedge dx = - dx \wedge dy\text{.}\) So while we can exchange the order of the iterated integrals once we have written everything in terms of double integrals, when we write the integral in terms of differential forms, we must use the correct choice of basic two-form \(dx \wedge dy\) with the differentials appearing in the same order as the coordinates of \(\mathbb{R}^2\text{.}\) That's because integrals of two-forms are oriented, while double integrals of functions are not.

Note that there is nothing special about the variables \(x\) and \(y\text{.}\) We could name the coordinates of \(\mathbb{R}^2\) anything. For instance, if we choose coordinates \((u,v)\) on \(\mathbb{R}^2\text{,}\) then the integral of a two-form \(\omega = f \ du \wedge dv\) over a region \(D\) is equal to the double integral of \(f\) over \(D\) (with a positive sign in front) when \(D\) is endowed with the canonical orientation described by the ordered basis \(\{ \mathbf{e}_1, \mathbf{e}_2 \} \text{,}\) with the basis vectors \(\mathbf{e}_1=(1,0)\) and \(\mathbf{e}_2=(0,1)\) pointing in the positive directions of the coordinates \(u\) and \(v\) respectively. 2 

In a more careful treatment, we would allow to write the two-form \(\omega\) in terms of any basic two-form \(\alpha\text{.}\) We would then define the orientation that is compatible with the basic two-form as corresponding to choices of oriented bases \(\{ \mathbf{u}, \mathbf{v} \}\) such that \(\alpha(\mathbf{u}, \mathbf{v}) > 0\text{.}\) Our definition would then state that the integral of a two-form written in terms of this basic two-form over a region \(D\) with the compatible orientation would be given by the double integral of the function over that region, while the integral over the region with opposite orientation would be given by minus the double integral.

Subsection 5.3.2 Integrals of two-forms over regions in \(\mathbb{R}^2\) are oriented and reparametrization-invariant

We already mentioned that integrals of two-forms oriented. Let us now be a little more precise, and show that, with our definition, integrals of two-forms are invariant under orientation-preserving reparametrizations, and change sign under orientation-reversing reparametrizations.

Let us first define what we mean by orientation-preserving and orientation-reversing reparametrizations. We state the definition in \(\mathbb{R}^2\text{,}\) but it naturally generalizes to \(\mathbb{R}^n\text{.}\)

Definition 5.3.6. Orientation-preserving reparametrizations of regions in \(\mathbb{R}^2\).

Let \(D_1, D_2 \subset \mathbb{R}^2\) be recursively supported regions, 3  and \(\phi: D_2 \to D_1\) a bijective function that can be extended to a \(C^1\)-function on an open subset \(U \subseteq \mathbb{R}^2\) that contains \(D_2\text{.}\) We assume that \(\phi\) is invertible (except possibly on the boundary of \(D_2\)), i.e. the determinant of its Jacobian \(\det J_\phi\) is non-zero on the interior of \(D_2\) (recall the definition of the Jacobian in Definition 4.7.5). We say that \(\phi\) is an orientation-preserving reparametrization if \(\det J_\phi > 0\) for all points in the interior of \(D_2\text{,}\) and that it is an orientation-reserving reparametrization if \(\det J_\phi \lt 0\) for all points in the interiori of \(D_2\text{.}\)

We could extend the statement more generally to closed bounded regions by taking unions of recursively supported regions.

With this definition, we can now show that our theory of integration for two-forms over regions in \(\mathbb{R}^2\) is oriented and reparametrization-invariant, our two guiding principles.

The key is to use Lemma 4.7.7 to calculate \(\phi^* \omega\text{.}\) Since we are pulling back a two-form \(\omega\) on an open subset \(U \subseteq \mathbb{R}^2\text{,}\) and that \(\phi: V \to U \) with \(V \subseteq \mathbb{R}^2\text{,}\) we can apply the result of Lemma 4.7.7. We write \(\omega = f\ dx \wedge dy\text{,}\) and \(\phi(u,v) = (x(u,v), y(u,v) )\text{.}\) The jacobian of the transformation is

\begin{equation*} J_\phi = \begin{pmatrix} \frac{\partial x}{\partial u} \amp \frac{\partial x}{\partial v} \\ \frac{\partial y}{\partial u} \amp \frac{\partial y}{\partial v} \end{pmatrix}. \end{equation*}

Then Lemma 4.7.7 tells us that

\begin{equation*} \phi^* \omega = f(\phi(u,v) ) (\det J_\phi)\ du \wedge dv. \end{equation*}

This means that

\begin{equation*} \int_{D_2} \phi^* \omega = \int_{D_2} f(\phi(u,v) ) (\det J_\phi)\ du \wedge dv. \end{equation*}

Using our definition of integration in Definition 5.3.1, we know that

\begin{equation*} \int_{D_1} \omega = \iint_{D_1} f(x,y)\ dx dy, \qquad \int_{D_2} \phi^* \omega = \iint_{D_2} f(\phi(u,v) ) (\det J_\phi)\ du dv, \end{equation*}

where on the right-hand-side of each equation we mean the double integral for the recursively supported regions \(D_1\) in the \(xy\)-plane and \(D_2\) in the \(uv\)-plane.

But recall from your previous calculus course that double integrals satisfy a “transformation formula”, or “change of variables formula”, which is the natural generalization of the substitution formula for definite integrals. The transformation formula states that

\begin{equation*} \iint_{D_1} f(x,y)\ dx dy = \iint_{D_2} f(\phi(u,v) ) |\det J_\phi| du dv. \end{equation*}

Note that there is now an absolute value around the determinant of the Jacobian. Thus, what this means is that if our transformation is such that \(\det J_\phi >0\text{,}\) then \(|\det J_\phi| = \det J_\phi\text{,}\) and

\begin{equation*} \int_{D_1} \omega = \iint_{D_1} f(x,y)\ dx dy= \int_{D_2} f(\phi(u,v) ) (\det J_\phi)\ du \wedge dv =\int_{D_2} \phi^* \omega, \end{equation*}

while if \(\det J_\phi \lt 0\text{,}\) then \(|\det J_\phi| = - \det J_\phi\text{,}\) and

\begin{equation*} \int_{D_1} \omega = \iint_{D_1} f(x,y)\ dx dy= - \int_{D_2} f(\phi(u,v) ) (\det J_\phi)\ du \wedge dv = -\int_{D_2} \phi^* \omega. \end{equation*}

This is the statement of the lemma: integrals of two-forms are oriented and reparametrization-invariant!

What is particularly nice with the proof of the lemma is that the transformation (or change of variables) formula for double integrals is simply the statement that integrals of two-forms over regions in \(\mathbb{R}^2\) are invariant under orientation-preserving reparametrizations! Isn't that cool? It explains why the determinant of the Jacobian appears; it comes from pulling back the two-form under the change of variables.

The fact that double integrals involve the absolute value of the determinant of the Jacobian, while our integrals do not (and change signs under orientation-reversing reparametrizations), is also interesting. As alluded to above, the reason is that our integrals are oriented, while standard double integrals in calculus are not.

Consider the basic two-form \(\omega = dx \wedge dy\) on \(\mathbb{R}^2\text{.}\) Let us define the following \(x\)-supported domain:

\begin{equation*} D = \{ (x,y) \in \mathbb{R}^2 ~|~ x \in [-1,1], -\sqrt{1-x^2} \leq y \leq \sqrt{1-x^2} \}. \end{equation*}

It is easy to see that \(D\) is a closed disk of radius one centered at the origin. The integral of the basic two-form \(\omega\) over \(D\) should give us the area of the disk, namely \(\pi\text{.}\) We calculate:

\begin{align*} \int_D \omega =\amp \iint_D \ dA \\ =\amp \int_{-1}^1 \int_{-\sqrt{1-x^2}}^{\sqrt{1-x^2}} \ dy dx\\ =\amp \int_{-1}^1 \left( \sqrt{1-x^2} + \sqrt{1-x^2} \right)\ dx\\ =\amp 2 \int_{-1}^1 \sqrt{1-x^2}\ dx. \end{align*}

To evaluate this integral, we do a trigonometric substitution, \(x = \sin(\theta)\text{:}\)

\begin{align*} 2 \int_{-1}^1 \sqrt{1-x^2}\ dx =\amp 2 \int_{-\pi/2}^{\pi/2} \sqrt{1-\sin^2(\theta)} \cos(\theta)\ d \theta\\ =\amp 2 \int_{-\pi/2}^{\pi/2} \cos^2 (\theta)\ d \theta\\ =\amp \int_{-\pi/2}^{\pi/2} (1 + \cos(2 \theta))\ d \theta\\ =\amp \left( \frac{\pi}{2} + \frac{\pi}{2} \right) + \frac{1}{2} \left( \sin(\pi) - \sin(-\pi) \right)\\ =\amp \pi. \end{align*}

We could have instead use a change of variables to evaluate this integral: polar coordinates. Define the map \(\phi: \mathbb{R}^2 \to \mathbb{R}^2\) with

\begin{equation*} \phi(r, \theta) = (r \cos(\theta), r \sin(\theta) ). \end{equation*}

Then, if we define the region

\begin{equation*} D_2 =\{ (r, \theta) \in \mathbb{R}^2~|~ r \in [0,1], \theta \in [0,2\pi ] \}, \end{equation*}

the map \(\phi: D_2 \to D\) is bijective and invertible in the interior of \(D_2\text{.}\) The determinant of the Jacobian is

\begin{equation*} \det J_\phi = \det \begin{pmatrix} \cos(\theta) \amp - r \sin(\theta) \\ \sin(\theta) \amp r \cos(\theta) \end{pmatrix} = r \cos^2(\theta) + r \sin^2(\theta) = r, \end{equation*}

which is positive for all points in the interior of \(D_2\text{.}\) Thus \(\phi\) is an orientation-preserving reparametrization, so the integral of the pullback \(\phi^* \omega = r\ dr \wedge d \theta\) over \(D_2\) should give us \(\pi\) again. Indeed, we get:

\begin{align*} \int_{D_2} \phi^* \omega =\amp \int_{0}^{2 \pi} \int_0^1 (\det J_\phi )\ dr d\theta\\ =\amp \int_{0}^{2 \pi} \int_0^1 r\ dr d\theta\\ =\amp \frac{1}{2} \int_0^{2 \pi}d \theta\\ =\amp \pi. \end{align*}

Notice how easier the integral was! That's of course because polar coordinates are well suited for evaluating integrals over regions that have circular symmetry.

Exercises 5.3.3 Exercises

1.

Evaluate the integral of the two-form \(\omega = \frac{1}{(1+x+y)^2} dx \wedge dy\) over the rectangular region

\begin{equation*} D = \{ (x,y) \in \mathbb{R}^2\ | \ x \in [1,2], y \in [2,3] \} \end{equation*}

with canonical orientation.

Solution.

By definition of the integral, we have:

\begin{align*} \int_D \omega =\amp \int_D \frac{1}{(1+x+y)^2}\ dx \wedge dy\\ =\amp \int_2^3 \int_1^2 \frac{1}{(1+x+y)^2}\ dx dy\\ =\amp \int_2^3 \left[-\frac{1}{1+x+y} \right]_{x=1}^{x=2}\ dy\\ =\amp \int_2^3 \left( \frac{1}{2+y} - \frac{1}{3+y} \right)\ dy \\ =\amp \left[ \ln(2+y) - \ln(3+y) \right]_{y=2}^{y=3}\\ =\amp \ln(5)-\ln(6)-\ln(4)+\ln(5)\\ =\amp \ln \left(\frac{25}{24} \right). \end{align*}

2.

Evaluate the integral of the two-form \(\omega = \frac{\ln(y)}{x y}\ dx \wedge dy\) over the rectangular region

\begin{equation*} D = \{(x,y) \in \mathbb{R}^2\ | \ x \in [1,4], y \in [1,5] \} \end{equation*}

with clockwise orientation.

Solution.

By definition of the integral, we have (we add a minus sign since we are evaluating the integral with clockwise orientation):

\begin{align*} \int_D \omega =\amp - \int_D \frac{\ln(y)}{x y}\ dx \wedge dy\\ =\amp -\int_1^5 \int_1^4 \frac{\ln(y)}{x y}\ dx dy\\ =\amp -\int_1^5 \frac{\ln(y)}{y} \left[\ln(x) \right]_{x=1}^{x=4}\ dy\\ =\amp -\ln(4) \int_1^5 \frac{\ln(y)}{y}\ dy. \end{align*}

To evaluate the remaining definite integral, we do a substitution \(u = \ln(y)\text{,}\) \(du = \frac{1}{y}\ dy\text{,}\) and \(y=1\) goes to \(u=\ln(1)=0\text{,}\) while \(y=5\) goes to \(u = \ln(5)\text{.}\) We get:

\begin{align*} \int_D \omega =\amp -\ln(4) \int_0^{\ln(5)} u\ du\\ =\amp -\ln(4) \left[ \frac{u^2}{2} \right]_{u=0}^{u=\ln(5)}\\ =\amp -\frac{\ln(4) (\ln(5))^2}{2}. \end{align*}

3.

Evaluate the integral of the two-form \(\omega = \cos(y^3)\ dx \wedge dy\) over the region \(D\) bounded by the parabola \(x=y^2\) and the lines \(x=0\text{,}\) \(y=0\) and \(y=1\text{,}\) with canonical orientation.

Solution.

We first describe the region explicitly. It is shown in the figure below:

The region.
Figure 5.3.9. The region \(D\) is the region bounded by the two curves shown above, the \(y\)-axis, and the \(x\)-axis. The green curve is the upper half of the curve \(x = y^2\text{,}\) while the orange line is the line \(y=1\text{.}\)
We can describe this region as an \(x\)-supported or \(y\)-supported region. If we write it as a \(y\)-supported region, we would write \(0 \leq y \leq 1\text{,}\) and \(0 \leq x \leq y^2\text{.}\) Thus

\begin{equation*} D = \{(x,y) \in \mathbb{R}^2\ | \ y \in [0,1], 0 \leq x \leq y^2 \}. \end{equation*}

If we write it as an \(x\)-supported region, we would write \(0 \leq x \leq 1\text{,}\) and \(\sqrt{x} \leq y \leq 1\text{,}\) that is,

\begin{equation*} D = \{(x,y) \in \mathbb{R}^2\ | \ x \in [0,1], \sqrt{x} \leq y \leq 1 \}. \end{equation*}

Which of the two descriptions should we choose? If we use the second description (\(x\)-supported), we first need to integrate in \(y\text{.}\) But the function that we have to integrate is \(\cos(y^3)\text{,}\) which has no elementary antiderivative. So this is problematic. We will not run into this problem if we use the first description (\(y\)-supported); with this description, we first integrate in \(x\text{,}\) which is fine. And the next integral in \(y\) will then be easy, as we will see.

So we use the first description of the region as \(x\)-supported. The integral then becomes

\begin{align*} \int_D \omega =\amp \int_D \cos(y^3)\ dx \wedge dy\\ =\amp \int_0^1 \int_0^{y^2} \cos(y^3)\ dx dy\\ =\amp \int_0^1 \cos(y^3) \left[x \right]_{x=0}^{x=y^2}\ dy\\ =\amp \int_0^1 y^2 \cos(y^3)\ dy. \end{align*}

To evaluate the remaining definite integral, we do the substitution \(u = y^3\text{,}\) \(du = 3 y^2\ dy\text{.}\) Then \(y=0\) becomes \(u=0\text{,}\) and \(y=1\) becomes \(u=1\text{.}\) We get:

\begin{align*} \int_D \omega =\amp \frac{1}{3} \int_0^1 \cos(u)\ du\\ =\amp \frac{1}{3} \sin(u) \Big|_{u=0}^{u=1}\\ =\amp \frac{1}{3} \sin(1). \end{align*}

4.

Consider the two-form \(\omega = e^{(x^2+y^2)^2}\ dx \wedge dy\text{,}\) and let \(D\) be the disk of radius one with counterclockwise orientation. Show that

\begin{equation*} \int_D \omega = 2 \pi \int_0^1 r e^{r^4}\ dr. \end{equation*}
Solution.

To show this, we change coordinates from Cartesian coordinates to polar coordinates. We define the function \(\phi: D_2 \to D\text{,}\) with

\begin{equation*} \phi(r,\theta) = (r \cos(\theta), r \sin(\theta) ), \end{equation*}

and

\begin{equation*} D_2 = \{ (r,\theta) \in \mathbb{R}^2 \ | \ r \in [0,1], \theta \in [0, 2 \pi] \}. \end{equation*}

It maps the rectangular region \(D_2\) to the unit disk \(D\text{.}\) The determinant of the Jacobian of \(\phi\) is

\begin{align*} \det J_\phi =\amp \det \begin{pmatrix} \frac{\partial x}{\partial r} \amp \frac{\partial x }{\partial \theta} \\ \frac{\partial y}{\partial r} \amp \frac{\partial y}{\partial \theta} \end{pmatrix}\\ =\amp \det \begin{pmatrix} \cos(\theta) \amp -r \sin(\theta) \\ \sin(\theta) \amp r \cos(\theta) \end{pmatrix}\\ =\amp r. \end{align*}

As \(r \in [0,1]\text{,}\) this is positive on the interior of \(D_2\text{,}\) and thus we know that the integral will be invariant under pullback, that is,

\begin{equation*} \int_{D_2} \phi^* \omega = \int_D \omega. \end{equation*}

We then calculate the pullback two-form \(\phi^* \omega\text{:}\)

\begin{align*} \phi^* \omega =\amp e^{r^4} (\det J_\phi) dr \wedge d \theta\\ =\amp r e^{r^4}\ dr \wedge d\theta. \end{align*}

Therefore,

\begin{align*} \int_D \omega =\amp \int_{D_2} r e^{r^4}\ dr \wedge d \theta\\ =\amp \int_0^1 \int_0^{2 \pi} r e^{r^4}\ d \theta d r\\ =\amp \int_0^1 r e^{r^4} \left[ \theta \right]_{\theta=0}^{\theta=2 \pi} dr\\ =\amp 2 \pi \int_0^1 r e^{r^4}\ dr. \end{align*}

5.

Consider the two-form \(\omega = x y\ dx \wedge dy\text{,}\) and let \(D\) be the region bounded by the four curves \(y=\frac{1}{x}\text{,}\) \(y=\frac{4}{x}\text{,}\) \(y=x\) and \(y = 4 x\text{,}\) with canonical orientation. Evaluate the integral of \(\omega\) over \(D\) using the change of variables

\begin{equation*} (x,y) = \left(\frac{u}{v}, v \right). \end{equation*}
Solution.

The region \(D\) is shown in the figure below:

The region.
Figure 5.3.10. The region \(D\) is the region bounded by the three curves shown above. The blue curve is the curve \(y=1/x\text{,}\) the orange curve is \(y=4/x\text{,}\) the green line is \(y=x\text{,}\) and the red line is \(y=4 x\text{.}\)
We could evaluate the integral of \(\omega\) over \(D\) by splitting the region \(D\) into two sub-regions, and then realizing these sub-regions as \(x\)-supported or \(y\)-supported. Or, we can do a change of variables, as specified in the question.

We consider the change of variables \(\phi: D_2 \to D\) with

\begin{equation*} \phi(u,v) = \left( \frac{u}{v}, v \right). \end{equation*}

What is the domain \(D_2\) in the \((u,v)\)-plane such that \(\phi(D_2) = D\text{?}\) In the variables \(u\) and \(v\text{,}\) the four curves \(y=1/x\text{,}\) \(y=4/x\text{,}\) \(y=x\) and \(y= 4x\) become

\begin{equation*} v = \frac{v}{u}, \qquad v= \frac{4 v}{u}, \qquad, v = \frac{u}{v}, \qquad v = 4 \frac{u}{v}. \end{equation*}

Those four bounding equations can be rewritten as

\begin{equation*} u = 1,\qquad u=4, \qquad v= \sqrt{u}, \qquad v = 2 \sqrt{u}. \end{equation*}

Here we used the fact that \(v\) is positive, since \(y=v\) is positive. So we can describe the region \(D_2\) in the \((u,v)\)-plane as the \(v\)-supported region

\begin{equation*} D_2 = \{(u,v) \in \mathbb{R}^2\ |\ u \in [1,4], \sqrt{u} \leq v \leq 2 \sqrt{u} \}. \end{equation*}

This region is shown in the figure below:

The region.
Figure 5.3.11. The region \(D_2\) in the \((u,v)\)-plane is the region bounded by the four curves shown above. The blue curve is the curve \(v = \sqrt{u}\text{,}\) the orange curve is \(v = 2 \sqrt{u}\text{,}\) and the two grey vertical lines are \(u=1\) and \(u=4\text{.}\)

Next, we calculate the determinant of the Jacobian of \(\phi\text{.}\) We get:

\begin{align*} \det J_\phi =\amp \det \begin{pmatrix} \frac{\partial x}{\partial u} \amp \frac{\partial x}{\partial v} \\ \frac{\partial y}{\partial u} \amp \frac{\partial y}{\partial v} \end{pmatrix} \\ =\amp \det \begin{pmatrix} \frac{1}{v} \amp - \frac{u}{v^2} \\ 0 \amp 1 \end{pmatrix}\\ =\amp \frac{1}{v}. \end{align*}

As \(v > 0\) on \(D_2\text{,}\) the determinant of the Jacobian is positive on \(D_2\text{.}\) So we know that the integral is invariant under pullback:

\begin{equation*} \int_D \omega = \int_{D_2} \phi^* \omega. \end{equation*}

The pullback is:

\begin{align*} \phi^* \omega =\amp u (\det J_\phi)\ du \wedge dv\\ =\amp \frac{u}{v}\ du \wedge dv. \end{align*}

We then integrate:

\begin{align*} \int_{D_2} \phi^* \omega =\amp \int_{D_2} \frac{u}{v}\ du \wedge dv\\ =\amp \int_1^4 \int_{\sqrt{u}}^{2 \sqrt{u}} \frac{u}{v} \ dv du\\ =\amp \int_1^4 u \left[ \ln(v) \right]_{v=\sqrt{u}}^{v=2\sqrt{u}}\ du\\ =\amp \int_1^4 u \left( \ln(2 \sqrt{u}) - \ln(\sqrt{u}) \right)\ du\\ =\amp \ln(2) \int_1^4 u\ du\\ =\amp \ln(2) \left( \frac{16}{2} - \frac{1}{2} \right)\\ =\amp \frac{15 \ln(2)}{2}. \end{align*}

For fun, let us show that we would get the same thing by evaluating the integral of \(\omega = x y\ dx \wedge dy\) directly by splitting the region into two recursively supported sub-regions. Looking at the graph in Figure 5.3.10, one choice is to split \(D\) into the two \(x\)-supported regions \(D_1\) and \(D_2\text{,}\) with

\begin{align*} D_1 =\amp \{(x,y) \in \mathbb{R}^2\ | \ x \in [1/2, 1], 1/x \leq y \leq 4 x \},\\ D_2 =\amp \{(x,y) \in \mathbb{R}^2\ | \ x \in [1,2], x \leq y \leq 4/x \}. \end{align*}

We can then evaluate both integrals separately and add them up to get the result. For \(D_1\text{,}\) we calculate:

\begin{align*} \int_{D_1} \omega =\amp \int_{1/2}^1 \int_{1/x}^{4 x} x y\ dy dx\\ =\amp \int_{1/2}^1 x \left[ \frac{y^2}{2} \right]_{y=1/x}^{y= 4x}\ dx\\ =\amp \int_{1/2}^1 x \left(8 x^2 - \frac{1}{2 x^2} \right)\ dx\\ =\amp \left[ 2 x^4 - \frac{1}{2} \ln(x) \right]_{x=1/2}^{x=1}\\ =\amp \frac{15}{8} - \frac{1}{2} \ln(2). \end{align*}

As for \(D_2\text{,}\) we get:

\begin{align*} \int_{D_2} \omega =\amp \int_{1}^2 \int_{x}^{4/x} x y\ dy dx\\ =\amp \int_{1}^1 x \left[ \frac{y^2}{2} \right]_{y=x}^{y= 4/x}\ dx\\ =\amp \int_{1}^1 x \left(\frac{8}{x^2} - \frac{x^2}{2} \right)\ dx\\ =\amp \left[ 8 \ln(x) - \frac{x^4}{8} \ln(x) \right]_{x=1}^{x=2}\\ =\amp 8 \ln(2) - \frac{15}{8}. \end{align*}

Adding those two integrals, we get:

\begin{align*} \int_D \omega =\amp \int_{D_1} \omega + \int_{D_2} \omega\\ =\amp \frac{15}{8} - \frac{1}{2} \ln(2) + 8 \ln(2) - \frac{15}{8}\\ =\amp \frac{15 \ln(2)}{2}. \end{align*}

This is the same answer that we obtained previously via our change of variables. Great!