octave:2> P=randstoc(4);
octave:3> invariantPotenze(P,1e-6,100)
ans =

   0.27652   0.22941   0.24143   0.25264

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

   0.28000   0.24000   0.29000   0.19000

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

   0.29500   0.21000   0.24600   0.24900

octave:7> harry
v =

   0.45000   0.15000   0.40000

ans =  29.500
ans =

   0.45110   0.15210   0.39680

ans =  29.457
octave:8> diary off