series inversion
The method of series inversion allows us to derive the power series
of an inverse function f-1 given the power series of f.
Clearly, since we are representing function f and f-1 by power series,
the function f must necessarily be holomorphic; that is, it is differentiable
as a complex-valued function on an open subset of the complex plane
.
(It follows that f-1 must also be holomorphic.)
Holomorphic functions include the elementary functions
studied in calculus
such as sin, cos, tan and exp.
For the method to work smoothly, it is best to assume that we want to invert f at around the origin, and its value there is zero. There is no loss of generality, since if f(b)=a, we can apply series inversion to the function g defined by f(z)=a+g(z-b). Then f-1(w)=g-1(w-a)+b; we obtain the power series for f-1 centred at a.
Also, it must be true that f′(0)≠0, otherwise f will not even
be invertible around the origin.
An example
We explain the method by an example, for f(z)=arctanz. In the following, we will consistently use the notation O(zn) to denote a holomorphic function h(z) whose power series expansion begins with zn. And similarly when z is replaced with the variable w.
First, we start with the well-known power series expansion for w=arctanz:
w=z-z33+z55+O(z7). | (1) |
The number of explicit terms in the power series expansion determines the number of terms that we will be able to obtain in the power series expansion of f-1. So in this case, we are going to seek an expansion of f-1(w)=tanw up to (but excluding) the w7 term.
A simple rearrangement of (1) gives
z=w+z33-z55+O(z7) | (2) |
Now we substitute equation (2) into itself. Of course, usually when we substitute an equation into itself we do not get anything, but here it works because we can perform simplications using the O notation. So for instance, in the following, in second term z3/3 on the right of equation (2), we put in equation (2) simplified to z=w+O(z3). Why we should choose this simplication will be clear in a moment:
z | =w+13(w+O(z3))3-z55+O(z7) | (3) | ||
=w+13(w3+O(z3)(w2+⋯))-z55+O(z7) | (4) | |||
=w+13(w3+O(w3)⋅O(w2))+O(w5). | (5) |
In equation (5) we used the fact that the expansion for z=f-1(w) must begin with a w term, i.e. f-1(0)=0. Also note that we are guaranteed that the w and z terms have non-zero coefficients, because f′(0)≠0. (Otherwise we would not be able to isolate z in equation (2).)
Now equation (5) simplifies to
z=w+w33+O(w5), | (6) |
which is already an achievement, because we have identified exactly what the w3 term must be.
Summary
To summarize the procedure in general, we start with the expansion
w=f(z)=a1z+a2z2+a3z3+⋯,a1≠0, | (10) |
and rearrange it to,
z=b1(w-a2z2-a3z3-⋯-anzn)+O(wn+1),b1=1a1. | (11) |
So we know that z=b1w+O(w2), and we can substitute this into the term z2 of equation (11). At the end we will get an equation of the form z=b1w+b2w2+O(w3), and we can substitute this into the terms z2 and z3 of (11). And what ever results we will substitute back into the terms z2, z3, z4 of equation (11). We can repeat this process until we have all the terms of z=f-1(w) that we need.
We probably should normalize the functions so that a1=b1=1 to make the computations easier.
References
- 1 Lars V. Ahlfors. Complex Analysis. McGraw-Hill, 1979.
Title | series inversion |
---|---|
Canonical name | SeriesInversion |
Date of creation | 2013-03-22 15:39:09 |
Last modified on | 2013-03-22 15:39:09 |
Owner | stevecheng (10074) |
Last modified by | stevecheng (10074) |
Numerical id | 8 |
Author | stevecheng (10074) |
Entry type | Derivation |
Classification | msc 30B10 |