Tworzenie macierzy w numpy

Cześć.
Dopiero zaczynam swoją przygodę z Pythonem, więc proszę o wyrozumiałość :)
Moim celem jest stworzenie macierzy NxM przy czym moje N to macierz Nx1, a M to liczba całkowita. Dodatkowo elementy macierzy wynikowej są podnoszone odpowiednio do potęgi M.
Ja to robię tak:

def my_matrix(N, M):
    matrix_N = np.array(N)
    matrix = np.zeros(shape=(matrix_N.size, (M + 1)))
    for i in range(M + 1):
        for j in range(matrix_N.size):
            matrix[j, i] = matrix_N[j]**i
    return matrix

Teraz zastanawiam się czy jest jakiś prostszy sposób na wyliczenie takiej macierzy? Nie bardzo jeszcze znam API numpy.

Co tu dużo pisać, jak Zaczynasz przygodę, to ten tutorial najlepiej najpierw przerobić do spodu.