CoCalc
Sharedsupport / 2015-03-18-120814-sqrtn-mod-1.sagewsOpen in CoCalc

Examples for support purposes...

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)