miau

GC, GS = [], [] # canonical and systematic generator matrices (dimension 16x8)
for a in alpha:
    row_GC, row_GS = [], []
    for k in range(8):
        row_GC.append(a^k)
        p = lagrange_polys[k]
        row_GS.append(p(a))
    GC.append(row_GC)
    GS.append(row_GS)
GC, GS = matrix(GC), matrix(GS) # y = G*x, no x*G porque quedaron al reves