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