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) be defined on the domain, D, where D is a bounded region. D 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 D, with a rectangle, let’s call this rectangle for R.
We define a function, f~, on R as the following:
f~(x,y)={f(x,y)0(x,y)∈D(x,y)∈/D
Definition 1 (Double integral over a general region)
The double integral of, f over D is:
∬Df(x,y) dA=∬Rf~(x,y) dA
Definition 2 (Equivalent definition)
Recall our original definition of the double integral, we can also divide the domain, D, into smaller pieces.
Therefore, the double integral of, f, over D is also:
∬Df(x,y) dA=Δx→0Δy→0limj=1∑mi=1∑nf(xi\*,yj\*)(ΔA)ij
Since everything outside D, is 0, they contribute nothing to the overall volume. Therefore, the volume under f~ is the same as the volume under f.
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, D, is of type I, if it lies between two graphs of two continuous functions of x.
Meaning:
D={(x,y)∣a≤x≤b,g1(x)≤y≤g2(x)}
Let’s see how we can compute this now. Let D be of type I. Let’s enclose this region within a rectangle as well:
∬f(x,y) dA=∬Rf~(x,y) dA=∫ab∫cdf~(x,y) dy dx
If we fix x, then f~ becomes:
f~(x,y)=⎩⎨⎧0f(x,y)0c≤y≤g1(x)g1(x)≤y≤g2(x)g2(x)≤y≤d
Which means:
∫ab∫g1(x)g2(x)f(x,y) dy dx
Proposition 1 (Integration over a type I region)
If D is of type I, then:
∬Df(x,y) dA=∫ab∫g1(x)g2(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)≤x≤h2(y),c≤y≤d}
We perform the same procedure as type I, instead of fixing x, we fix y.
Proposition 2 (Integration over a type II region)
If D is of type II, then:
∬Df(x,y) dA=∫cd∫h1(y)h2(y)f(x,y) dx dy
Before we head into doing a few examples, it’s a good practice always draw where our D lies. So we easily get a geometrical sense of how it lies.
Example 1
Find ∬Dy dA, where D is the region bounded by, y=2x2,y=1+x2.
So, our region is bounded by two parabolas, we have two functions of x, therefore, our region is of type I.
Let’s firstly find where these graphs intersect:
2x2=1+x2x2=1x=±1Graphing these we also see that, y=1+x2 is the upper graph, and the lower is y=2x2.
∫−11∫2x21+x2y dy dx21∫−11y2y=2x2y=1+x2 dx21∫−11(1+x2)2−(2x2)2 dx21∫−111+2x2+x4−4x4 dx21∫−111+2x2−3x4 dx21[x+32x3−53x5x=−1x=1]21[(1+32−53)−(−1−32+53)]21[1+32−53+1+32−53]1+32−53
Example 2
Find ∬Dxy dA, where D is the region bounded by, y=x−1,y2=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=y2−1.
Let’s find where these intersect:
y+1=y2−1y2−y−2=0y1=…=−1y2=…=2Which means:
∫−12∫y2−1y+1xy dx dy∫−122y[x2]x=y2−1x=y+1 dy∫−122y[(y+1)2−(y2−1)2] dy∫−122y[y2+2y+1−y4+2y2−1] dy∫−122y[3y2+2y−y4] dy∫−1223y3+y2−2y5 dy83y4+3y3−12y6y=−1y=2=…=827
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
Computing this entire integral will not yield a singular number, it will yield a function of x instead.
However, if D 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, where D is the region bounded by, x=0,y=1,y=x
Graphing this we can see that this a region of both types.
For, type I:
∬Df(x,y) dA=∫01∫xysin(y2) dy dxType II:
∬Df(x,y) dA=∫01∫0ysin(y2) dx dyIn this case, one of these isn’t even integrable in terms of elementary functions. So yeah…
Solving the second integral:
∫01∫0ysin(y2) dx dy∫01x⋅sin(y2)x=0x=y dy∫01y⋅sin(y2) dyUsing normal u-sub, we obtain:
−21cos(y2)y=0y=1(−21cos(1))−(−21cos(0))−21cos(1)+2121(−cos(1)+1)