A(X) = RX + XT

Laplacian2D = kron(speye(n),D) + kron(D, speye(n));