Calculus – Multiple integrals – Double integrals over rectangles

Fubini's theorem

Theorem: If \(f\) is continuous on the rectangle

\[R=[a,b]\times[c,d]=\{(x,y)\,|\,a\leq x\leq b,\;c\leq y\leq d\},\]

then we have:

\[\iint\limits_Rf(x,y)\,dA=\int_a^b\int_c^df(x,y)\,dy\,dx=\int_c^d\int_a^bf(x,y)\,dx\,dy.\]

More generally, this is true if we assume that \(f\) is bounded on \(R\), \(f\) is discontinuous only on a finite number of smooth curves, and the iterated integrals exist.

Stewart §15.1, Example 5
Evaluate \(\displaystyle\iint\limits_R(x-3y^2)\,dA\), where \(R=[0,2]\times[1,2]\).

Solution: We have:

\[\iint\limits_R(x-3y^2)\,dA=\int_0^2\int_1^2(x-3y^2)\,dy\,dx=\int_0^2\bigg[xy-y^3\bigg]_{y=1}^2\,dx=\int_0^2(x-7)\,dx =\bigg[\frac{1}{2}x^2-7x\bigg]_{x=0}^2=-12\]

and

\[\iint\limits_R(x-3y^2)\,dA=\int_1^2\int_0^2(x-3y^2)\,dx\,dy=\int_1^2\bigg[\frac{1}{2}x^2-3xy^2\bigg]_{x=0}^2\,dy =\int_1^2\left(2-6y^2\right)\,dy=\bigg[2y-2y^3\bigg]_{y=1}^2=-12.\]

Stewart §15.1, Example 6
Evaluate \(\displaystyle\iint\limits_Ry\sin(xy)\,dA\), where \(R=[1,2]\times[0,\pi]\).

Solution: We have:

\[\iint\limits_Ry\sin(xy)\,dA=\int_0^{\pi}\int_1^2y\sin(xy)\,dx\,dy=\int_0^{\pi}\bigg[-\cos(xy)\bigg]_{x=1}^2\,dy =\int_0^{\pi}\left(\cos(x)-\cos(2x)\right)\,dx=\bigg[\sin(x)-\frac{1}{2}\sin(2x)\bigg]_{x=0}^{\pi}=0.\]

Here the other order is much more difficult. Using integration by parts we have

\[\int y\sin(xy)\,dy=-\frac{1}{x}\int y\,d\cos(xy)=-\frac{y}{x}\cos(xy)+\frac{1}{x}\int\cos(xy)\,dy=-\frac{y}{x}\cos(xy)+\frac{1}{x^2}\sin(xy)+C.\]

Hence we have:

\[\iint\limits_Ry\sin(xy)\,dA=\int_1^2\int_0^{\pi}y\sin(xy)\,dy\,dx=\int_1^2\bigg[-\frac{y}{x}\cos(xy)+\frac{1}{x^2}\sin(xy)\bigg]_{y=0}^{\pi} =\int_1^2\left(-\frac{\pi}{x}\cos(\pi x)+\frac{1}{x^2}\sin(\pi x)\right)\,dx.\]

Note that \(-\displaystyle\frac{1}{x}\sin(\pi x)\) is an antiderivative, since using the chain rule we have: \(\displaystyle\frac{d}{dx}\left(-\frac{1}{x}\sin(\pi x)\right)=-\frac{\pi}{x}\cos(\pi x)+\frac{1}{x^2}\sin(\pi x)\). Hence:

\[\iint\limits_Ry\sin(xy)\,dA=\int_1^2\int_0^{\pi}y\sin(xy)\,dy\,dx=\int_1^2\left(-\frac{\pi}{x}\cos(\pi x)+\frac{1}{x^2}\sin(\pi x)\right)\,dx =\bigg[-\frac{1}{x}\sin(\pi x)\bigg]_{x=1}^2=\sin(\pi)-\frac{1}{2}\sin(2\pi)=0.\]

Average value

The average value of a function \(f\) of one variable defined on an interval \([a,b]\) is \(\displaystyle\frac{1}{b-a}\int_a^bf(x)\,dx\).

Similarly, we have: the average value of a function \(f\) of two variables defined on a rectangle \(R\) is \(\displaystyle\frac{1}{\textrm{area}(R)}\iint\limits_Rf(x,y)\,dA\).

Catalan's constant

Earlier we considered Catalan's constant \(G\) which is defined as the sum of the series \(\displaystyle\sum_{n=0}^{\infty}\frac{(-1)^n}{(2n+1)^2}\).

We obtained several (integral) representations for \(G\), such as: \(\displaystyle\int_0^1\frac{\arctan(x)}{x}\,dx=G\).

Now this can be used to prove a nice representation as a "double integral": \(\displaystyle\int_0^1\int_0^1\frac{dx\,dy}{1+x^2y^2}=G\). The proof is straightforward:

\[\int_0^1\int_0^1\frac{dx\,dy}{1+x^2y^2}=\int_0^1\int_0^1\frac{dy\,dx}{1+x^2y^2}=\int_0^1\bigg[\frac{\arctan(xy)}{x}\bigg]_{y=0}^1\,dx=\int_0^1\frac{\arctan(x)}{x}\,dx=G.\]

Alternatively, we might use \(\displaystyle\frac{1}{1-t}=\sum_{n=0}^{\infty}t^n\) for \(|t|<1\) to find that

\[\int_0^1\int_0^1\frac{dx\,dy}{1+x^2y^2}=\sum_{n=0}^{\infty}(-1)^n\int_0^1\int_0^1(xy)^{2n}\,dx\,dy =\sum_{n=0}^{\infty}(-1)^n\left(\int_0^1t^{2n}\,dt\right)^2=\sum_{n=0}^{\infty}\frac{(-1)^n}{(2n+1)^2}=G.\]

---

Last modified on September 27, 2021