Invert a power series?

$\begingroup$

If one is given a power series of the form:

$$\sum^\infty_{m=1} a_m x^m = a_0 + a_1x + a_2 x^2 + \dots$$

for known $a_m$'s i.e $a_m = f(m)$ and as $m\rightarrow\infty$, $a_m\rightarrow0$ is there anyway to find out what the function is given by? For example, I know that

$$\sum^\infty_{m=1} x^m = {1\over 1-x}.$$

So is there an expression for for non-unity coefficients? Is there a way t0f doing this? Like inverting the Taylor series or something similar?

$\endgroup$ 3

2 Answers

$\begingroup$

There are two questions here, one in the title and one in the body of the text:

(i) Given a power series $$\sum_{k=0}^\infty a_k x^k\tag{1}$$ with $a_0=c\ne0$ in general there is no way to get hold of some functional inverse "in finite terms". Assume that this series defines a monotonic function $f:\ U\to{\mathbb R}$ in some neighborhood $U$ of $0$. The corresponding inverse function would then be defined in the neighborhood of $c$ and maybe have a development of the form $f^{-1}(y)=\sum_{l=0}^\infty b_l(y-c)^l$. Even when this is the case the $b_l$ cannot be obtained from the $a_k$ by some recursive formula involving only finitely many $a_k$ at each step.

If, however, $a_0=0$ and $a_1\ne0$ then the series $(1)$ has a "formal" inverse, as well as a functional inverse. This means the following: There is a sequence $(b_l)_{l\geq1}$ of coefficients $b_l$, defined recursively from the coefficients $a_k$, such that "formally" $$\sum_{l=1}^\infty b_l\left(\sum_{k=1}^\infty a_k x^k\right)^l=x\ .$$ Furthermore, if the series $(1)$ actually converges to some function $f(x)$ in a neighborhood $U$ of $0$, then the series $\sum_{l=1}^\infty b_ly^l$ represents $f^{-1}(y)$ in some neighborhood $V$ of $0$. Mathematica has a command that computes the $b_l$ from the $a_k$, no question of convergence involved.

(ii) A series $(1)$ involves an arbitrary infinite sequence of coefficients $a_k$, whereas any well known special function, as $\exp$, $x\mapsto x/(1-x)^2$ and the like, has a Taylor series whose coefficients follow some rigid law, say, the law of the Fibonacci numbers. Therefore you cannot expect that a "random" series of the form $(1)$ can be identified as the series of a function in your catalogue of "well known functions".

But note that there is some help when the encountered $a_k$ happen to be integers: You can look up an initial segment $(a_0,a_1,\ldots, a_n)$ of the sequence $(a_k)_{k\geq0}$ in the online catalogue OEIS and will maybe find out that it belongs to some "well known function".

$\endgroup$ $\begingroup$

You cannot "invert" a power series (P.S.) by inverting term by term. If you want to find the inverse of a function of which a powerseries representation is established, you got to kind of start all over. For example, there is a P.S. for sinx. But reciprocating these polynomial terms does not get the P.S. for the cosecant function. Instead you have to find a couple of derivatives of the cosecant and established the coefficients by putting in zero for x every time. Other mechanisms are also possible to find it, but still involves some algebra.

Sometimes you can build an inverse function from a given function. For example, a P.S. representation of arcsine can be found by creating √(1-x) in a denominator (or use Binomial Series) and then replacing x by x² etc.

$\endgroup$

Your Answer

Sign up or log in

Sign up using Google Sign up using Facebook Sign up using Email and Password

Post as a guest

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

You Might Also Like