def constrmat(n,m):
    Z=[]
    for i in range(n):
        Z+=[[]]
        for j in range(m):
            Z[i]+=[i*len(Z[0])+j+1]
    return Z


print constrmat(6,9)