Some time ago I had to do a calculation concerning a physical problem. In this calculation some integrals were needed to be solved and the method on how to do this was briefly described. I was able to solve the problem completely and found the method so beautiful that I would like to share it with all of you. Consider the following integrals:(adsbygoogle = window.adsbygoogle || []).push({});

[tex]\int \frac{1}{1+\epsilon \cdot cos \theta}d\theta[/tex]

[tex]\int \frac{1}{\left(1+\epsilon \cdot cos \theta\right)^2}d\theta[/tex]

[tex]\int \frac{1}{\left(1+\epsilon \cdot cos \theta\right)^3}d\theta[/tex]

[tex]\int \frac{cos \theta}{1+\epsilon \cdot cos \theta}d\theta[/tex]

[tex]\int \frac{cos \theta}{\left(1+\epsilon \cdot cos \theta\right)^2}d\theta[/tex]

[tex]\int \frac{cos \theta}{\left(1+\epsilon \cdot cos \theta\right)^3}d\theta[/tex]

[tex]\int \frac{sin \theta}{1+\epsilon \cdot cos \theta}d\theta[/tex]

[tex]\int \frac{sin \theta}{\left(1+\epsilon \cdot cos \theta\right)^2}d\theta[/tex]

[tex]\int \frac{sin \theta}{\left(1+\epsilon \cdot cos \theta\right)^3}d\theta[/tex]

Some of these are straightforward, but some are not. In order to solve the non trivial ones in a systematic way, one can use the following substitution:

[tex]1+\epsilon \cdot cos \theta = \frac{1-\epsilon^2}{1-\epsilon \cdot cos \gamma}[/tex]

With:

[tex]0\leq \theta \leq 2\pi[/tex]

[tex]0\leq \gamma \leq 2\pi[/tex]

The following relations can be obtained:

[tex]cos \theta = \frac{cos \gamma -\epsilon}{1-\epsilon \cdot cos \gamma}[/tex]

[tex]sin \theta = \frac{\sqrt{1-\epsilon^2} sin \gamma}{1-\epsilon \cdot cos \gamma}[/tex]

[tex]cos \gamma = \frac{\epsilon +cos \theta}{1+\epsilon \cdot cos \theta}[/tex]

[tex]sin \gamma = \frac{\sqrt{1-\epsilon^2} sin \theta}{1+\epsilon \cdot cos \theta}[/tex]

[tex]d \theta = \frac{\sqrt{1-\epsilon^2}}{1-\epsilon \cdot cos \gamma}d \gamma[/tex]

[tex]d \gamma = \frac{\sqrt{1-\epsilon^2}}{1+\epsilon \cdot cos \theta}d \theta[/tex]

After rewriting some of the integrals into smaller ones, substituting this and rearranging it is possible to solve them in a fairly easy way. As an example, let's take the one before the last, it was:

[tex]\int \frac{sin \theta}{\left(1+\epsilon \cdot cos \theta\right)^2}d\theta[/tex]

Using the substitution we get:

[tex]\int \frac{sin \gamma}{1-\epsilon^2}d\gamma=\frac{-cos \gamma}{1-\epsilon^2}+C[/tex]

The solution is now obtained by substituting either

[tex]\gamma = arccos \left( \frac{\epsilon +cos \theta}{1+\epsilon \cdot cos \theta} \right)[/tex]

or

[tex]\gamma = arcsin \left( \frac{\sqrt{1-\epsilon^2} sin \theta}{1+\epsilon \cdot cos \theta} \right)[/tex]

This substitution is called the "Sommerfeld substitution" after the "inventor" and is one of the nicest substitutions I've ever encountered for solving integrals in the real with "classical" functions. I hope this is helpful in solving other types of integrals you might be working on.

best regards,

Coomast

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Sharing knowledge

**Physics Forums | Science Articles, Homework Help, Discussion**