[UPDATE: 25/09/2014 with the iterative method] This kind of transcendental equation is often encountered in physics. Undergraduate students are usually shown (or asked to draw) the graphical solution:
The numerical solutions are easily found by an iterative method using a scientific calculator (see below), but how far can one go with only pen and paper?
Expansion
Aside from the trivial solution
x0=0, one clearly has
xk≃(2k+1)π2 (
k≥1), so we can write:
xk=(2k+1)π2−εk,withεk<1 One would like to do an expansion in
εk, but of course this will not work for the tangent around its divergence points. We can however use the cotangent, since
tan(xk)=xk⇒cot(xk)=1/xk. Using standard substitution formulas for the sine and cosine yields:
cot[(2k+1)π2−εk]=tan(εk)≃εk≃2(2k+1)π where in the last equality we neglected
εk in the denominator. One can include it for a more rigorous treatment. Finally, we have:
xk≃(2k+1)π2−2(2k+1)π,fork≥1, giving for the first three solutions 4.5002, 7.7267, and 10.9046, to be compared with the "exact" values 4.4934…, 7.7253…, and 10.9041…. The quality of the approximation increases with the order
k, since
εk decreases (the intersections are closer and closer to the vertical asymptotes).
Iteration
Let us rewrite the initial equation by applying the arctangent to both members:
x=arctan(x)
For the
k-th solution, the initial estimate is:
x0k=(2k+1)π2. Let us plug it in the right-hand side of Eq. (1) to obtain the first order estimate
x1k and then iterate. Note that the arctangent is a multi-valued function, and the standard implementation reduces it to the first branch (the one going through the origin). We are looking for the solution sitting on the
k-th branch, so we need to add
kπ each time:
xi+1k=arctan(xik)+kπ
For the first non-trivial solution (
k=1), the sequence is: 4.71239, 4.50328, 4.49387, 4.49343, 4.49341,... with the second iteration already reaching an excellent precision!