pythontexgodel

\setcounter{MaxMatrixCols}{20}

\begin{pycode}
import numpy as np

def Table_Layer(n):
    A=np.zeros((2**n+1,2**n+1))
    for i in range(2**n):
        A[i+1,1]=(i+2)%2**n
    for j in range(2**n):
        A[2**n,j+1]=(j+1)%2**n
    for i in range(2**n+1):
        for j in range(2**n+1):
            if A[i,j]==0:
                A[i,j]=2**n
    for i in range(2**n+1):
        for j in range(2**n+1):
            if (i<2**n) and (i>0):
                if j>1:
                    k=int(A[2**n-i,j-1])
                    p=int(A[2**n-i,1])
                    ci=int(2**n-i)
                    A[ci,j]=A[k,p]
    B=np.zeros((2**n,2**n))
    for i in range(2**n):
        for j in range(2**n):
            B[i,j]=int(A[i+1,j+1])
    print(r"\begin{pmatrix}")
    for i in range(2**n) :
        for j in range(2**n -1) :
            print(f"{int(B[i,j])}",'& ',end='')
        print(f"{int(B[i,2**n-1])}",'\\\\')
    print(r"\end{pmatrix}")
\end{pycode}

Pour $n=3$ :

\[\pyc{Table_Layer(3)}\]

Pour $n=4$ :

\[\pyc{Table_Layer(4)}\]