v=rand(1,4)

v =

    0.8147    0.9058    0.1270    0.9134

v=v/sum(v)

v =

    0.2951    0.3281    0.0460    0.3308

u=rand

u =

    0.6324

cumsum(v)

ans =

    0.2951    0.6232    0.6692    1.0000

u<cumsum(v)

ans =

     0     0     1     1

find(u<cumsum(v))

ans =

     3     4

find(u<cumsum(v),1)

ans =

     3

v

v =

    0.2951    0.3281    0.0460    0.3308

randdisc(v)

ans =

     1

randdisc(v)

ans =

     1

randdisc(v)

ans =

     2

randdisc(v)

ans =

     4

randdisc(v)

ans =

     4

randdisc(v)

ans =

     1

testsimula
testsimula
m=10;
x=linspace(pi/2,2*pi,m);
plot(x,H(x,1:m))
hold off
plot(x,H(x,1:m))
plot(x,H(x,1:m),'*')
m=100;
x=linspace(pi/2,2*pi,m);
plot(x,H(x,1:m),'*')
sa(100)

counterstati =

    66

Elapsed time is 0.021127 seconds.

ans =

  -32.8900

sa(100)

counterstati =

    57

Elapsed time is 0.014559 seconds.

ans =

  -32.8900

sa(100)

counterstati =

    28

Elapsed time is 0.006138 seconds.

ans =

  -33.3802

sa(100)

counterstati =

    23

Elapsed time is 0.006076 seconds.

ans =

  -33.3802

sa(100)

counterstati =

    24

Elapsed time is 0.006990 seconds.

ans =

  -32.2947

sa(100)

counterstati =

    46

Elapsed time is 0.007439 seconds.

ans =

  -32.2947

sa(200)

counterstati =

    42

Elapsed time is 0.014086 seconds.

ans =

  -31.6268

sa(200)

counterstati =

    61

Elapsed time is 0.014877 seconds.

ans =

  -33.3723

sa(200)

counterstati =

    28

Elapsed time is 0.015405 seconds.

ans =

  -33.3724

sa(200)

counterstati =

    41

Elapsed time is 0.014804 seconds.

ans =

  -32.3660

sa(200)

counterstati =

    25

Elapsed time is 0.012778 seconds.

ans =

  -33.2021

sa(400)

counterstati =

    48

Elapsed time is 0.038655 seconds.

ans =

  -33.2159

sa(400)

counterstati =

    41

Elapsed time is 0.036173 seconds.

ans =

  -33.3678

sa(400)

counterstati =

    24

Elapsed time is 0.036437 seconds.

ans =

  -30.0196

sa(400)

counterstati =

    33

Elapsed time is 0.036753 seconds.

ans =

  -29.6626

sa(400)

counterstati =

    35

Elapsed time is 0.034804 seconds.

ans =

  -31.7554

sa(400)

counterstati =

    22

Elapsed time is 0.033241 seconds.

ans =

  -32.5879

sa(400)

counterstati =

    52

Elapsed time is 0.035522 seconds.

ans =

  -32.8312

sa(400)

counterstati =

    25

Elapsed time is 0.038361 seconds.

ans =

  -27.0521

sa(800)

counterstati =

    44

Elapsed time is 0.105658 seconds.

ans =

  -33.2944

sa(800)

counterstati =

    35

Elapsed time is 0.102667 seconds.

ans =

  -31.9521

sa(800)

counterstati =

    26

Elapsed time is 0.098106 seconds.

ans =

  -23.3247

sa(800)

counterstati =

    42

Elapsed time is 0.104982 seconds.

ans =

  -33.3353