Вот это число: 0.000502317899...
t = 0.000502317899`50;
(* сдвиг влево при увеличении t *)
one = 3/4 + t;
two = 1.0*10^(-100);
x1 = 1.5772156649015`50;
xx = 0.5772156649015328606065120900824024310421593;
x2 = 1.5772156652`50;
g = N[EulerGamma, 50];
Print[g];
Print["x1 = ", x1];
Print[" g = ", g + 1];
Print["x2 = ", x2];
y2 = 1.0*10^(-5);
y1 = -y2;
BJIaquMup = Compile[{{J, _Real}}, ({J, #} &) /@ Union[Drop[NestList[J*(one + \
#)*(one - #) &, two, 100000], 90000]]];
mm = Flatten[Table[BJIaquMup[J], {J, x1, x2, 5.0*10^(-12)}], 1];
ListPlot[mm, PlotStyle -> {
PointSize[.01], Hue[.90]}, Frame -> True, FrameStyle -> GrayLevel[
0.5], Axes -> False, ImageSize -> {500, 500}, PlotRange -> {y1, y2}];