Conceptually, computing double integrals in polar coordinates is the
same as in rectangular coordinates. After all, the idea of an integral
doesn't depend on the coordinate system. If R is a region in the plane
and f(x,y) is a function, then ∬Rf(x,y)dA is what we get when we
Chop R into a bunch of small pieces.
Compute f(x∗,y∗)ΔA for each piece, where ΔA is
the area of the piece.
Add up the contributions of the pieces, and
Take a limit as we chop R into smaller and smaller pieces.
The differences when working with polar coordinates are explained
in the following video:
When working with rectangular coordinates, our pieces are boxes of width
Δx, height Δy, and area ΔA=ΔxΔy.
When working with polar coordinates, our pieces are polar rectangles.
These are regions of the form θ∈[α,β], r∈[a,b].
The approximate
contribution of a single box to our integral is f(r∗i,θ∗j)ΔA, where ΔA is the area of the box. [Here we are thinking
of f as a function of r and θ. If f is a function of x and
y, then strictly speaking we should
write f(rcos(θ),rsin(θ)) instead of
f(r,θ).]
So what is the area of a polar rectangle? The polar rectangle is the
difference of two pie wedges, one with radius b and one with radius a,
and so has area
ΔA=b22(β−α)−a22(β−α)=b2−a22(β−α)=a+b2(b−a)(β−α)=r∗ΔrΔθ,
where r∗=a+b2 is the average value of r,
Δr=b−a is the change in r and Δθ=β−α
is the change in θ.
Another way to understand this is to look at the shape of a small
polar rectangle, like the shaded box in the figure above. This is
almost a rectangular, only rotated and with slightly curved sides. One side
has length Δr. The other had length approximately rΔθ.
So the whole rectangle has approximate area rΔrΔθ.
Now we organize our boxes. Since the contribution of each box is
approximately f(r∗i,θ∗j)r∗iΔrΔθ,
the contribution of the boxes in the same thin pie wedge is approximately
(∫baf(r,θ∗j)rdr)Δθ.
Finally we put all the pie wedges together, to get a total of
∫βα(∫baf(r,θ)rdr)dθ.
The infinitestimal area in rectangular coordinates is dA=dxdy. The
infinitesimal area in polar coordinates is dA=rdrdθ.
Don't forget the factor of r!!
Double integrals in polar coordinates become iterated integrals
∬f(x,y)dA=∫∫f(rcos(θ),rsin(θ))rdrdθ=∫∫f(rcos(θ),rsin(θ))rdθdr,