| Hosted by CoCalc | Download
n=3 sn = RealField(200)(sqrt(3)) v = [float(sn^a-int(sn^a)) for a in range(10000)]
v.sort() t = stats.TimeSeries(v)
d = t.diffs() d = d.clip_remove(0.01) d.plot_histogram(bins=100)