MVA[R_,mu_,K_] := Module[{n=Length[R],EE,ee,la,T,nn},
  EE=Table[1,{n},{n}]-IdentityMatrix[n];
  ee=Table[1,{n}];
  nn=Table[0,{n}];
  la=ee.Inverse[R+EE];
  Do[T=(1+nn)/mu; nn=k*la*T/(la.T), {k,1,K}];
  {nn,T}]