Necro-post ... Since those servers are gone, anyone else searching may find
this useful. Based on William's original interact and modified to have the
Sage-wrapped mpl histograms. Still amazing :)
>
>
def sqrt2():
PI = float(pi)
return plot(lambda x: (2/PI)*math.sqrt(1-x^2), -1,1, plot_points=200,
rgbcolor=(0.3,0.1,0.1), thickness=2)
delta = delta_qexp(10^5)
@interact
def delta_dist(bins=(20,[10..150]), number = [500,1000,..,delta.prec()]):
D = delta[:number]
w = [float(D[p])/(2*float(p)^(5.5)) for p in prime_range(number + 1)]
show(histogram(w, bins=bins, normed=True) + sqrt2(), frame=True,
gridlines=True)
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at https://groups.google.com/group/sage-support.
For more options, visit https://groups.google.com/d/optout.
