Nice. But your plots are line-wrapped. Out[698]= 1 + 2^(-1 - k) Sqrt[2^k + 2 Sqrt[1 + 2^(2 (-1 + k))]] You have, e.g. In[700]:= ArrayPlot[ Partition[RealDigits[%698 /. k -> 63, 2, 63^2][[1]], 63]] Try In[705]:= ArrayPlot[ Partition[ Join[ConstantArray[0, 32], RealDigits[%698 /. k -> 63, 2, 2*63^2][[1]]], 2*63]] --rwg Violating two embargoes at once: Climate change makes planes disappear.