И эту псевдоальфу 0.005675 даёт всё та же псевдопарабола (x+o)*(x-o)
ae = 0.00567504`50;
a = 1.18370401243`50;
o = a*ae;
two = 1.0`*10^(-50);
x1 = 176.2101`50;
c = 0.00002`50;
x2 = 1/ae + c;
Print["x2 = ", x2];
Print["o = ", o];
y2 = 3.0*10^(-6);
y1 = -y2;
BJIaquMup = Compile[{{J, _Real}}, ({J, #} &) /@ Union[Drop[
NestList[J*(o + #)*(o - #) &, two, 7000], 6000]]];
mm = Flatten[Table[BJIaquMup[J], {J, x1, x2, 2.0*10^(-6)}], 1];
ListPlot[mm, PlotStyle -> {PointSize[.01],
Hue[.55]}, Frame ->
True, FrameStyle -> GrayLevel[0.5], Axes -> False,
ImageSize -> {500, 500}, PlotRange -> {y1, y2}]