octave:2> simulatedannealingsave
Elapsed time is 4.55774 seconds.
i =  87
ans = -30.786
valutazioni =  41
octave:3> simulatedannealingsave
Elapsed time is 3.25 seconds.
i =  90
ans = -33.380
valutazioni =  29
octave:4> simulatedannealingsave
Elapsed time is 3.93 seconds.
i =  88
ans = -32.295
valutazioni =  36
octave:5> simulatedannealingsave
Elapsed time is 2.39 seconds.
i =  88
ans = -32.295
valutazioni =  21
octave:6> simulatedannealingsave
Elapsed time is 3.91 seconds.
i =  180
ans = -33.372
valutazioni =  35
octave:7> simulatedannealingsave
Elapsed time is 8.46 seconds.
i =  177
ans = -32.865
valutazioni =  77
octave:8> simulatedannealingsave
Elapsed time is 3.48 seconds.
i =  184
ans = -31.627
valutazioni =  31
octave:9> v=rand(1,4)
v =

   0.83484   0.18557   0.69736   0.31499

octave:10> v=v/sum(v)
v =

   0.410690   0.091291   0.343062   0.154957

octave:11> v*ones(4)
ans =

   1   1   1   1

octave:12> P=[0 1/3 2/3;1/3 0 2/3;1 0 0]
P =

   0.00000   0.33333   0.66667
   0.33333   0.00000   0.66667
   1.00000   0.00000   0.00000

octave:13> ones(1,3)/(eye(3)-P+ones(3))
ans =

   0.45000   0.15000   0.40000

octave:14> simulazione(P,200,@f)
ans =  1
octave:15> [i,s]=simulazione(P,200,@f)
i =  2
s =  0.37500
octave:16> [i,s]=simulazione(P,200,@f)
i =  3
s =  0.42000
octave:17> [i,s]=simulazione(P,200,@f)
i =  1
s =  0.40500
octave:18> [i,s]=simulazione(P,200,@f)
i =  2
s =  0.40000
octave:19> [i,s]=simulazione(P,200,@f)
i =  1
s =  0.39500
octave:20> [i,s]=simulazione(P,200,@f)
i =  3
s =  0.39500
octave:21> [i,s]=simulazione(P,200,@f)
i =  1
s =  0.37500
octave:22> [i,s]=simulazione(P,200,@f)
i =  1
s =  0.37500
octave:23> [i,s]=simulazione(P,200,@f)
i =  1
s =  0.40000
octave:24> [i,s]=simulazione(P,500,@f)
i =  1
s =  0.40000
octave:25> [i,s]=simulazione(P,500,@f)
i =  3
s =  0.39000
octave:26> [i,s]=simulazione(P,500,@f)
i =  1
s =  0.42000
octave:27> [i,s]=simulazione(P,500,@f)
i =  1
s =  0.38400
octave:28> [i,s]=simulazione(P,500,@f)
i =  2
s =  0.39800
octave:29> [i,s]=simulazione(P,500,@f)
i =  2
s =  0.40600
octave:30> [i,s]=simulazione(P,500,@f)
i =  3
s =  0.40000
octave:31> [i,s]=simulazione(P,500,@f)
i =  1
s =  0.39400
octave:32> [i,s]=simulazione(P,500,@f)
i =  3
s =  0.38800
octave:33> [i,s]=simulazione(P,500,@f)
i =  2
s =  0.41000
octave:34> [i,s]=simulazione(P,500,@f)
i =  1
s =  0.40400
octave:35> v
v =

   0.410690   0.091291   0.343062   0.154957

octave:36> v=ones(1,3)/(eye(3)-P+ones(3))
v =

   0.45000   0.15000   0.40000

octave:37> f([1 2 3])
ans =

   20   30   40

octave:38> f([1 2 3])*v'
ans =  29.500
octave:39> [i,s]=simulazione(P,500,@f)
i =  3
s =  29.660
octave:40> [i,s]=simulazione(P,500,@f)
i =  1
s =  29.680
octave:41> [i,s]=simulazione(P,500,@f)
i =  3
s =  29.540
octave:42> [i,s]=simulazione(P,500,@f)
i =  1
s =  29.580
octave:43> [i,s]=simulazione(P,500,@f)
i =  3
s =  29.560
octave:44> [i,s]=simulazione(P,500,@f)
i =  2
s =  29.520
octave:45> lambda=2.5
lambda =  2.5000
octave:46> k=[0:20];
octave:47> a=exp(-lambda)*lambda.^k./factorial(k);
octave:48> bar(k,a)
octave:49> a
a =

 Columns 1 through 5:

   8.2085e-02   2.0521e-01   2.5652e-01   2.1376e-01   1.3360e-01

 Columns 6 through 10:

   6.6801e-02   2.7834e-02   9.9406e-03   3.1064e-03   8.6290e-04

 Columns 11 through 15:

   2.1573e-04   4.9028e-05   1.0214e-05   1.9643e-06   3.5076e-07

 Columns 16 through 20:

   5.8461e-08   9.1345e-09   1.3433e-09   1.8657e-10   2.4549e-11

 Column 21:

   3.0686e-12

octave:50> diary off