Necesito ayuda con la generación de una secuencia a partir de 

    P = np.array([[0,0],
                  [1,0],
                  [0.5,np.sqrt(3)/2]])

`P` consta de tres entradas, cada una es un par ordenado que indica el vértice de un triángulo unitario

Con este triángulo necesito generar tres sub-triángulos, los vértices están en la lista `T` por eso obtengo una lista con 3^2 = 9 array. La cuestión es que quiero generar los siguientes sub-triángulos, lo que implica aplicar las mismas transformaciones a la nueva lista `T` para generar los nuevos 3^3 =27 vértices.

Este es el código que tengo actualmente:

    P = np.array([[0,0],[1,0],[0.5,np.sqrt(3)/2]])

    m = []
    for p in P:
        m.append(p*0.5)

    n=[]
    for p in P:
        n.append(p*0.5 + np.array([0.5,0]))

    k=[]
    for p in P:
        k.append(p*0.5 + np.array([0.25,np.sqrt(3)/4]))
    
    T = np.array([m,n,k])

donde `T` es de la forma 

    T = [ 
        array([[0.,  0.], [0.5, 0.], [0.25, 0.4330127]]), 
        array([[0.5,  0.], [1, 0.], [0.75, 0.4330127]]), 
        array([[0.25, 0.4330127], [0.75, 0.4330127],[0.5, 0.8660254]]) ]

El objetivo es obtener el triángulo de Sierpinski y graficarlo utilizando la función de matplolib siguiente: 

    mpatches.Polygon(T[i], fc="y")

Saludos!