%r

attach(faithful)
long <- eruptions[eruptions > 3]
par(pty="s")   # square figure
qqnorm(long)
qqline(long)

12 + 30

show(icosahedron(color='green', opacity=.5, mesh=3), spin=1)

stats.TimeSeries(1000).randomize('normal').sums().plot()