Hilbert curve iteration

hilbert.py

from pylab import *

def hilbert(pts):
    s = 1 + (pts[-1][0] - pts[0][0]) / 2
    return ( [(y-s, x-s) for x,y in pts]
           + [(x-s, y+s) for x,y in pts]
           + [(x+s, y+s) for x,y in pts]
           + [(s-y,-x-s) for x,y in pts] )

pts = [(0, 0)]
for i in range(8):
    pts = hilbert(pts)
    plot(*zip(*pts))
    show()