How to implement a transition matrix as a numpy array and generate a Markov chain from it?

Question

I want to use python to generate a realisation of the Markov chain using the transition matrix given as

T_{ij} =

0.40 & 0.56 & 0.03 & 0.01
0.