%I#10 2016年11月4日03:08:13
%S 2,0,6,5,9,2,2,2,0,2,3,7,6,6,2,1,8,8,8,1,0,4,6,1,2,5,4,1,0,1,8,
%温度4,3,0,0,1,4,2,4,9,8,5,3,1,9,0,6,7,3,2,8,5,1,7,9,1,8,0,5,9,8,
%U 8,5,9,4,9,3,2,6,0,6,7,7,7、7,8,3,5,5,4,5,7,0,2,8,7,1,5,9、8,4,6,7、8,6
%N对数积分li(z)在C中不动点的虚部的十进制展开式。
%C真实部分以及详细注释和链接请参见A276762。
%H Stanislav Sykora,n的表格,n=1..2000的a(n)</a>
%电子邮箱:2.0659222023706621889881046112541084301424985319006732838579118。。。
%t RealDigits[Im[z/.FindRoot[LogIntegral[z]==z,{z,2+I},WorkingPrecision->100]][1](*_Vaclav Kotesovec_,2016年10月30日*)
%o(PARI)\\z可以是t_INT、t_REAL或t_COMPLEX,但0或1除外
%o li(z)=
%o个{
%o my(sgn=(-1)^if(实(z)<1,imag(z)<0,imag;
%o sgn*Pi*I-eint1(-log(z));
%o}(o)
%o默认值(realprecision,2100);\\执行:
%o Eps_=4.0*10.0^(-default(realprecision));
%o z=1+I;zlast=0;\\初始化和迭代
%o表示(k=1,1e6,z=li(z));如果(abs(z-zlast)<Eps_,break);zlast=z);
%o imag(z)\\显示结果
%Y参见A276762(实部),A070769。
%K nonn,cons公司
%O 1,1号机组
%A _Stanislav Sykora_,2016年10月28日