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