octave:2> randdisc([1/6,1/3,1/2]) ans = 3 octave:3> randdisc([1/6,1/3,1/2]) ans = 2 octave:4> randdisc([1/6,1/3,1/2]) ans = 1 octave:5> randdisc([1/6,1/3,1/2]) ans = 3 octave:6> randdisc([1/6,1/3,1/2]) ans = 3 octave:7> randdisc([1/6,1/3,1/2]) ans = 3 octave:8> randdisc([1/6,1/3,1/2]) ans = 3 octave:9> randdisc([1/6,1/3,1/2]) ans = 2 octave:10> randdisc([1/6,1/3,1/2]) ans = 1 octave:11> randdisc([1/6,1/3,1/2]) ans = 2 octave:12> Q = toeplitz((4:-1:1).^2); octave:13> Q Q = 16 9 4 1 9 16 9 4 4 9 16 9 1 4 9 16 octave:14> sumQ2 = sum(Q,2) sumQ2 = 30 38 38 30 octave:15> Q = Q+diag(max(sumQ2)-sumQ2) Q = 24 9 4 1 9 16 9 4 4 9 16 9 1 4 9 24 octave:16> sum(Q,2) ans = 38 38 38 38 octave:17> Q = Q/max(sumQ2) Q = 0.631579 0.236842 0.105263 0.026316 0.236842 0.421053 0.236842 0.105263 0.105263 0.236842 0.421053 0.236842 0.026316 0.105263 0.236842 0.631579 octave:18> sum(Q) ans = 1.00000 1.00000 1.00000 1.00000 octave:19> sa2(@(i) H(i,100)) counterstati = 159 ans = -30.371 octave:20> sa2(@(i) H(i,200)) counterstati = 86 ans = -35.453 octave:21> sa2(@(i) H(i,500)) counterstati = 163 ans = -38.076 octave:22> sa2(@(i) H(i,1000)) counterstati = 91 ans = -38.704 octave:23> diary off