octave:2> A=randn(4) A = 0.97806 -0.74282 -1.57662 -0.53546 -1.14182 -0.25662 -2.40505 -0.88764 0.21792 1.05508 0.28486 0.83630 -0.23240 -1.42744 -0.26754 -1.70392 octave:3> (A*A*A*A*A*A*A-A^7) ans = -3.5527e-15 -1.7764e-15 -7.1054e-15 -5.9952e-15 3.5527e-15 -5.3291e-15 1.3545e-14 7.7716e-15 -4.4409e-15 -1.3323e-14 -1.0658e-14 -1.1546e-14 1.4655e-14 3.0198e-14 2.8422e-14 2.8422e-14 octave:4> A2=A*A; octave:5> A4=A2*A2; octave:6> (A*A2*A4-A^7) ans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 octave:7> A4*(A2*A)-A^7 ans = -3.5527e-15 -3.5527e-15 -1.2434e-14 -6.2172e-15 3.5527e-15 -3.5527e-15 -7.1054e-15 1.7764e-15 -8.8818e-16 2.6645e-15 1.5987e-14 1.0658e-14 5.3291e-15 1.7764e-15 -1.4211e-14 -7.1054e-15 octave:8> (A*A2)*A4-A^7 ans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 octave:9> A*(A2*A4)-A^7 ans = -7.1054e-15 1.7764e-15 -7.1054e-15 -2.6645e-15 1.7764e-15 1.7764e-15 -1.0658e-14 -1.7764e-15 2.6645e-15 -8.8818e-16 8.8818e-15 1.7764e-15 -1.7764e-15 1.7764e-15 -1.4211e-14 -3.5527e-15 octave:10> P=randstoc(4); octave:11> rigaP=@(i) P(i,:); octave:12> invariantPotenze(P,1e-6,100) ans = 0.21106 0.22157 0.31990 0.24747 octave:13> octave:13> v=invariantPotenze(P,1e-6,100) v = 0.21106 0.22157 0.31990 0.24747 octave:14> v*P ans = 0.21106 0.22157 0.31990 0.24747 octave:15> invariantSIM(rigaP,100) ans = 0.18000 0.21000 0.36000 0.25000 octave:16> invariantSIM(rigaP,1000) ans = 0.21800 0.23000 0.32300 0.22900 octave:17> invariantSIM(rigaP,10000) ans = 0.21380 0.22070 0.32910 0.23640 octave:18> invariantSIM(rigaP,10000) ans = 0.21120 0.21640 0.32370 0.24870 octave:19> harry v = 0.45000 0.15000 0.40000 ans = 29.500 ans = 0.44990 0.15210 0.39800 ans = 29.481 octave:20> harry v = 0.45000 0.15000 0.40000 ans = 29.500 ans = 0.45020 0.15120 0.39860 ans = 29.484 octave:21> diary off