GfT = np.zeros((n + 2) * (n + 1))