DtMarkovChain[P_] := Module[{n=Length[P],EE,ee}, EE=Table[1,{n},{n}]-IdentityMatrix[n]; ee=Table[1,{n}]; Return[ee.Inverse[P+EE]//Chop]]