octave:2> P=randstoc(4)
P =

   0.25308   0.16481   0.12034   0.46177
   0.19110   0.10095   0.28960   0.41834
   0.16877   0.26009   0.31528   0.25586
   0.13469   0.19339   0.31957   0.35235

octave:3> [zeros(1,4),1]/[(P-eye(4)),ones(4,1)]
ans =

   0.17563   0.18939   0.27771   0.35727

octave:4> rigaP = @(i) P(i,:);
octave:5> rigaP(3)
ans =

   0.16877   0.26009   0.31528   0.25586

octave:6> invariantSIM(rigaP,1000)
ans =

   0.19300   0.17800   0.27500   0.35400

octave:7> invariantSIM(rigaP,10000)
ans =

   0.16980   0.19430   0.27560   0.36030

octave:8> 
octave:8> harry
v =

   0.45000   0.15000   0.40000

ans =  29.500
ans =

   0.45070   0.15240   0.39690

ans =  29.462
octave:9> diary off