illustrations for RH book
%load code.sage v = prime_gap_distribution(10^7)[:50] b = bar_chart(v)
b
b.show(fontsize=18, ticks=[[2,4,6,8]+[10,20,..,40]+[48],[20000,60000,90000]], tick_formatter=['latex','latex'])
%load code.sage draw('primegapdist')
%load code.sage draw('sieves')
race_mult_parity(10)
@interact def _(s='prime_pi'): %load code.sage for t in s.split(): print t draw(t)