How to do the simulation in MATLAB for this type of transition matrix using inverse?

Hello,

Here is my transition matrix A=[1/3, 0, 1/3, 1/3;1/3, 1/3, 0, 1/3;1/3, 0,1/3,1/3; 0,1/3,1/3,1/3]

N=inv(I-A)

Can anyone show me how to simulate this N (Markov Chain) for 1000 and 5000 times? How to do the simulation in MATLAB for this type of matrix?