Calculus/Polar Integration

Introduction

Integrating a polar equation requires a different approach than integration under the Cartesian system, hence yielding a different formula, which is not as straight forward as integrating the function f(x).

Proof

In creating the concept of integration, we used riemann sums on rectangles to approximate the area under the curve. However, with polar graphs, one must use triangles that start from the origin, and have a radius ending on the curve. If you don't mind skipping the proof, this is the form to use to integrate polar equations:

$\frac{1}{2}\underset{a}{\overset{b}{\int }}r(\theta {)}^{2}d\theta$

Whereas the $(a,r(a))$ and $(b,r(b))$ are the ends of the curve that you wish to integrate.

Integral calculus

The integration region R is bounded by the curve r(θ) and the rays θ = a and θ = b.

Let R denote the region enclosed by a curve r(θ) and the rays θ = a and θ = b, where 0 < b − a < 2π. Then, the area of R is

$\frac{1}{2}\underset{a}{\overset{b}{\int }}r(\theta {)}^{2}d\theta .$

The region R is approximated by n sectors (here, n = 5).

This result can be found as follows. First, the interval [a, b] is divided into n subintervals, where n is an arbitrary positive integer. Thus Δθ, the length of each subinterval, is equal to b − a (the total length of the interval), divided by n, the number of subintervals. For each subinterval i = 1, 2, …, n, let θ_i be the midpoint of the subinterval, and construct a circular sector with the center at the pole, radius r(θ_i), central angle Δθ and arc length $r({\theta }_i)\mathrm{\Delta }\theta$. The area of each constructed sector is therefore equal to $\frac{1}{2}r({\theta }_i{)}^2\mathrm{\Delta }\theta$. Hence, the total area of all of the sectors is

$\sum _{i=1}^n\frac{1}{2}r({\theta }_i{)}^2\mathrm{\Delta }\theta .$

As the number of subintervals n is increased, the approximation of the area continues to improve. In the limit as n → ∞, the sum becomes the Riemann sum for the above integral.

Generalization

Using Cartesian coordinates, an infinitesimal area element can be calculated as dA = dx dy. The substitution rule for multiple integrals states that, when using other coordinates, the Jacobian determinant of the coordinate conversion formula has to be considered:

$J=\det \frac{\partial (x,y)}{\partial (r,\theta )}=\left|\begin{array}{cc}\frac{\partial x}{\partial r}& \frac{\partial x}{\partial \theta }\\ \frac{\partial y}{\partial r}& \frac{\partial y}{\partial \theta }\end{array}\right|=\left|\begin{array}{cc}\cos \theta & -r\sin \theta \\ \sin \theta & r\cos \theta \end{array}\right|=r{\cos }^2\theta +r{\sin }^2\theta =r.$

Hence, an area element in polar coordinates can be written as

$dA=Jdrd\theta =rdrd\theta .$

Now, a function that is given in polar coordinates can be integrated as follows:

$\underset{R}{\iint }f(r,\theta )dA=\underset{a}{\overset{b}{\int }}\underset{0}{\overset{r(\theta )}{\int }}f(r,\theta )rdrd\theta .$

Here, R is the same region as above, namely, the region enclosed by a curve r(θ) and the rays θ = a and θ = b.

The formula for the area of R mentioned above is retrieved by taking f identically equal to 1.

Applications

Polar integration is often useful when the corresponding integral is either difficult or impossible to do with the Cartesian coordinates. For example, let's try to find the area of the closed unit circle. That is, the area of the region enclosed by $x^2+y^2=1$.

In Cartesian

$\underset{-1}{\overset{1}{\int }}\underset{-\sqrt{1-x^2}}{\overset{\sqrt{1-x^2}}{\int }}dydx=2\underset{-1}{\overset{1}{\int }}\sqrt{1-x^2}dx$

In order to evaluate this, one usually uses trigonometric substitution. By setting $sin(\theta )=x$, we get both $cos(\theta )=\sqrt{1-x^2}$ and $cos(\theta )d\theta =dx$.

$\int \sqrt{1-x^2}dx=\int co{s}^2(\theta )d\theta =\int (\frac{1}{2}+\frac{1}{2}cos(2\theta ))d\theta =\frac{\theta }{2}+\frac{1}{4}sin(2\theta )=\frac{\theta }{2}+\frac{1}{2}sin(\theta )cos(\theta )=\frac{arcsin(x)}{2}+\frac{x\sqrt{1-x^2}}{2}$

Putting this back into the equation, we get

$2\underset{-1}{\overset{1}{\int }}\sqrt{1-x^2}dx=2[\frac{arcsin(x)}{2}+\frac{x\sqrt{1-x^2}}{2}{]}_{-1}^1=arcsin(1)-arcsin(-1)=\pi$

In Polar

To integrate in polar coordinates, we first realize $r=\sqrt{x^2+y^2}=\sqrt{1}=1$ and in order to include the whole circle, $a=0$ and $b=2\pi$

$\underset{0}{\overset{2\pi }{\int }}\underset{0}{\overset{1}{\int }}rdrd\theta =\underset{0}{\overset{2\pi }{\int }}[\frac{r^2}{2}{]}_0^{1}d\theta =\underset{0}{\overset{2\pi }{\int }}\frac{1}{2}d\theta =[\frac{\theta }{2}{]}_0^{2\pi }=\frac{2\pi }{2}=\pi$

An interesting example

A less intuitive application of polar integration yields the Gaussian integral

$\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{e}^{-x^2}dx=\sqrt{\pi }.$

Try it! (Hint: multiply $\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{e}^{-x^2}dx$ and $\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{e}^{-y^2}dy$)