from __future__ import division
import matplotlib.pyplot as plt
plt.figure()
plt.subplot(1,1,1)
plt.axis('off')

def drawBlue(p1,p2,p3):
  plt.fill([p1[0],p2[0],p3[0]],[p1[1],p2[1],p3[1]],facecolor='b',edgecolor='none')

def drawWhite(p1,p2,p3):
   plt.fill([p1[0],p2[0],p3[0]],[p1[1],p2[1],p3[1]],facecolor='w',edgecolor='none')

def midpoint(p1,p2):
  return ((p1[0]+p2[0])/2,(p1[1]+p2[1])/2)
def sierpinski(p1,p2,p3,repeat):
  drawWhite(midpoint(p1,p2),midpoint(p2,p3),midpoint(p3,p1))
  if (repeat <= 8):
    sierpinski(p1,midpoint(p1,p2),midpoint(p1,p3),repeat+1)
    sierpinski(p2,midpoint(p1,p2),midpoint(p2,p3),repeat+1)
    sierpinski(p3,midpoint(p1,p3),midpoint(p2,p3),repeat+1)

drawBlue((0,0),(7,0),(3.5,6.0621778265))

sierpinski((0,0),(7,0),(3.5,6.0621778265),0)

plt.show()