login
This site is supported by donations to The OEIS Foundation.

 

Logo

"Email this user" was broken Aug 14 to 9am Aug 16. If you sent someone a message in this period, please send it again.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A276762 Decimal expansion of the real part of a fixed point of the logarithmic integral li(z) in C. 2
1, 8, 7, 8, 8, 8, 1, 7, 4, 7, 9, 0, 8, 1, 2, 3, 0, 9, 1, 9, 6, 9, 4, 8, 6, 4, 9, 7, 3, 2, 6, 8, 8, 3, 4, 6, 9, 6, 8, 9, 5, 9, 4, 9, 5, 2, 8, 5, 5, 6, 3, 6, 6, 0, 6, 8, 9, 4, 5, 2, 7, 5, 1, 5, 0, 8, 5, 2, 1, 2, 5, 5, 4, 2, 6, 4, 7, 2, 6, 6, 9, 8, 8, 5, 0, 4, 8, 4, 5, 1, 2, 7, 0, 8, 9, 4, 4, 8, 5, 3, 5, 0, 8, 0, 6 (list; constant; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

The function li(z) has in C two fixed points, namely z0 = A276762 + A276763*i and its conjugate. For real arguments, z = 0 is also a fixed point but, since it lies on a branch cut and has discontinuous imaginary part, setting li(0) = 0 in C is problematic. This leaves z_0 and its conjugate as the only 'regular' fixed points of li(z). They are both attractors of the mapping. The attraction basin of z0 appears to be the whole upper half of the complex plane (including the real axis, minus the singular points z = 0 and z =1), while that of the conjugate of z0 appears to be the whole lower half (exluding the real axis). The convergence is exponential but rather slow, with the asymptotic convergence factor of 0.756330... per iteration, equal to abs(li'(z0)) = 1/abs(log(z0)).

LINKS

Stanislav Sykora, Table of n, a(n) for n = 1..2000

Eric Weisstein's World of Mathematics, Logarithmic Integral.

Wikipedia, Logarithmic integral function.

FORMULA

Setting z0 = A276762 + A276763*i, li(z0) = z0.

EXAMPLE

1.878881747908123091969486497326883469689594952855636606894527515...

MATHEMATICA

RealDigits[Re[z/.FindRoot[LogIntegral[z] == z, {z, 2+I}, WorkingPrecision -> 100]]][[1]] (* Vaclav Kotesovec, Oct 30 2016 *)

PROG

(PARI) \\ z may be t_INT, t_REAL, or t_COMPLEX except 0 or 1

li(z)=

{

  my(sgn=(-1)^if(real(z)<1, imag(z)<0, imag(z)<=0));

  sgn*Pi*I - eint1(-log(z));

}

default(realprecision, 2100); \\ Execution:

Eps_= 4.0*10.0^(-default(realprecision));

z=1+I; zlast=0; \\ Initialize and iterate

for(k=1, 1e6, z=li(z); if(abs(z-zlast)<Eps_, break); zlast=z);

real(z) \\ Display the result

CROSSREFS

Cf. A276763 (imaginary part), A070769.

Sequence in context: A037077 A094106 A277064 * A256609 A269892 A242968

Adjacent sequences:  A276759 A276760 A276761 * A276763 A276764 A276765

KEYWORD

nonn,cons

AUTHOR

Stanislav Sykora, Oct 28 2016

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy .

Last modified August 17 23:58 EDT 2017. Contains 290682 sequences.