ls 2013-10-18.txt invariantSIM.m 2013-10-25.txt metropolis.m 2013-11-08.txt randdisc.m 2013-11-15.txt randstoc.m epsmet.m sa.m fraccont.m~ simula.m harry.m testsimula.m H.m testsimula.m~ invariantPotenze.m edit invariantSIM.m v=[0,0,0,0] v = 0 0 0 0 v(1)=v(1)+1 v = 1 0 0 0 v(3)=v(3)+1 v = 1 0 1 0 v(1)=v(1)+1 v = 2 0 1 0 v(5)=v(5)+1 ??? Attempted to access v(5); index out of bounds because numel(v)=4. v(5)=1 v = 2 0 1 0 1 v(4)=v(4)+1 v = 2 0 1 1 1 v(8)=1 v = Columns 1 through 7 2 0 1 1 1 0 0 Column 8 1 P=randstoc(4); v=([P'-eye(4);ones(1,4)]\[zeros(4,1);1])' v = 0.2548 0.1739 0.2826 0.2887 v*P ans = 0.2548 0.1739 0.2826 0.2887 rigaP=@(i) P(i,:); rigaP=@(i) [1-sum(s(1:i)),s(i:-1:1)]; s ??? Undefined function or variable 's'. rigaP=@(i) P(i,:); invarianSIM(rigaP,1000) ??? Undefined function or method 'invarianSIM' for input arguments of type 'function_handle'. invariantSIM(rigaP,1000) ans = 0.2650 0.1810 0.2800 0.2740 v v = 0.2548 0.1739 0.2826 0.2887 P=randstoc(4);,v=([P'-eye(4);ones(1,4)]\[zeros(4,1);1])' v = 0.3245 0.2069 0.1810 0.2876 invariantSIM(rigaP,1000) ans = 0.2480 0.1630 0.2880 0.3010 rigaP=@(i) P(i,:); invariantSIM(rigaP,1000) ans = 0.3180 0.2140 0.1760 0.2920 v v = 0.3245 0.2069 0.1810 0.2876 invariantSIM(rigaP,10000) ans = 0.3245 0.2064 0.1829 0.2862 v v = 0.3245 0.2069 0.1810 0.2876 invariantSIM(rigaP,100000) ans = 0.3244 0.2054 0.1824 0.2878 v v = 0.3245 0.2069 0.1810 0.2876 edit harry.m % Modify expression to add input arguments. % Example: % a = [1 2 3; 4 5 6]; % foo(a); harry v = 0.4500 0.1500 0.4000 ans = 29.5000 ans = 0.4400 0.1700 0.3900 ans = 29.5000 harry v = 0.4500 0.1500 0.4000 ans = 29.5000 ans = 0.4600 0.1500 0.3900 ans = 29.3000 sum(w) ans = 100 sum(w/N) ans = 1 v v = 0.4500 0.1500 0.4000 w/N ans = 0.4600 0.1500 0.3900 sum([0.449,0.149,0.39]) ans = 0.9880 1-sum([0.449,0.149,0.39]) ans = 0.0120 0.45-0.449 ans = 0.0010 sum([0.449,0.149,0.399]) ans = 0.9970 P=randstoc(4);,v=([P'-eye(4);ones(1,4)]\[zeros(4,1);1])' v = 0.2339 0.3380 0.3061 0.1220 v=invariantSIM(rigaP,1000) ??? Undefined function or variable 'rigaP'. rigaP=@(i) P(i,:); v=invariantSIM(rigaP,1000) v = 0.2290 0.3300 0.3190 0.1220 sum(v) ans = 1 v=invariantSIM(rigaP,1000) v = 0.2460 0.3370 0.2860 0.1310 w=[0.229,0.337] w = 0.2290 0.3370 v=invariantSIM(rigaP,1000) v = 0.2310 0.3300 0.3230 0.1160 w=[0.229,0.337,0.323] w = 0.2290 0.3370 0.3230 v=invariantSIM(rigaP,1000) v = 0.2320 0.3180 0.3200 0.1300 w=[0.229,0.337,0.323,0.13] w = 0.2290 0.3370 0.3230 0.1300 sum(w) ans = 1.0190 randn ans = -0.4326 diary off