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]]