Re: question about the inverse li function
- To: mathgroup at smc.vnet.net
- Subject: [mg65379] Re: [mg65348] question about the inverse li function
- From: "Carl K. Woll" <carlw at wolfram.com>
- Date: Wed, 29 Mar 2006 06:34:18 -0500 (EST)
- References: <200603280905.EAA26284@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Capet Arthur wrote:
> the question...
>
> i've measured C(x), wich is equal to li(I(x)), where li denotes the
> logarithmic integral function. I would like to compute I(x)= li^(-1)
> (C(x))
>
> how can i compute the inverse of the logarithmic integral function ?
> Is there a function Inverse[_function] ?
>
> thanx a lot
>
>
> Arthur Capet, ULG, Belgium
>
>
Arthur,
There is the function InverseFunction. For example:
In[1]:=
InverseFunction[Sin]
Out[1]=
ArcSin
Unfortunuately, Mathematica does not have a built in function for the
inverse of LogIntegral:
In[2]:=
InverseFunction[LogIntegral]
Out[2]=
InverseFunction[LogIntegral]
So, we'll have to work harder to figure out the inverse. This question
has come up on mathgroup before, with 3 basic possibilities:
1. Plot LogIntegral[x] and then extract the plot points. Use
Interpolation to construct an InterpolatingFunction using these plot points.
2. Use FindRoot to find the inverse at a bunch of points, and then use
Interpolation.
3. Use homotopic continuation. That is, find the inverse of LogIntegral
at some point x, and then construct an ODE to determine how the inverse
varies as x varies.
I am usually the person advocating method 3. I'll show you how it works
in your example. First, form the equation for the inverse x[y]:
eqn = y == LogIntegral[ x[y] ]
Now, construct an ODE satisfied by the inverse, x:
In[4]:=
D[eqn,y]
Out[4]=
x'[y]
1 == ---------
Log[x[y]]
The next step is to use NDSolve on this ODE, and so we need an initial
value. However, the inverse is multivalued for negative input, since
LogIntegral ranges from 0 to -Infinity as x goes from 0 to 1, and then
LogIntegral ranges from -Infinity to Infinity as x goes from 1 to
Infinity. So, we'll need to construct two inverses: one for 0<x<1 and
one for 1<x<Infinity. For 0<x<1 we'll use the initial value
x[LogIntegral[1/2]]==1/2
and for 1<x<Infinity we'll use
x[LogIntegral[2]]==2
Now we are ready to use NDSolve.
low = x /.
NDSolve[{D[eqn, y], x[LogIntegral[1/2]] == 1/2}, x, {y, -5, 5}][[1]];
high = x /.
NDSolve[{D[eqn, y], x[LogIntegral[2]] == 2}, x, {y, -5, 5}][[1]];
Some tests:
In[7]:=
LogIntegral[high[2]]
LogIntegral[high[-2]]
LogIntegral[low[-2]]
Out[7]=
2.
Out[8]=
-2.
Out[9]=
-2.
As you may have noticed, I used the range -5<y<5. You can change this
range to whatever you need. Also, as I pointed out earlier, there are no
real values of x on the low branch for which LogIntegral[x] is positive.
There are complex values, but the branch chosen by my implementation
doesn't agree with the branch used by Mathematica for LogIntegral:
In[10]:=
LogIntegral[low[1]]
Out[10]=
1. + 3.14159 I
Finally, we can use low or high to answer your question:
i[x_] := low[ C[x] ]
or
i[x_] := high[ C[x] ]
Carl Woll
Wolfram Research
- References:
- question about the inverse li function
- From: "Capet Arthur" <Arthur.Capet@student.ulg.ac.be>
- question about the inverse li function