import matplotlib.pyplot as plt
import numpy as np
from itertools import combinations

def center(points):
	return sum(points)/len(points)

def draw_triangle(points):
	for ps in combinations(points, 2):
		xs = [ps[0][0], ps[1][0]]
		ys = [ps[0][1], ps[1][1]]        
		plt.plot(xs, ys, 'b-', alpha=0.1)

def mesh(points, depth):
	if depth > 0:
		c = center(points)
		for pair in combinations(points, 2):
			pts = [pair[0], pair[1], c]
			draw_triangle(pts)
			mesh(pts, depth-1)
			
points = [
	np.array([0, 1]),
	np.array([-0.866, -0.5]),
	np.array([ 0.866, -0.5]) ]

mesh(points, 3)
plt.axis("off")
plt.gca().set_aspect("equal")
plt.show()