def METHOD_NAME(a, m, n):
    f = lambda i, j: a[i] ** j
    A = np.fromfunction(np.vectorize(f), (m - 1, n), dtype=int)
    f = lambda i, j: (1 if j == (n - 1) else 0)
    z = np.fromfunction(np.vectorize(f), (1, n), dtype=int)
    return np.append(A, z, axis=0)