octave:2> v=[1,2,3,4,5] v = 1 2 3 4 5 octave:3> find(2.5<4) ans = 1 octave:4> find(2.5 v=[10,20,30,40,50] v = 10 20 30 40 50 octave:6> find(25 find(25 v=rand(1,4) v = 0.97177 0.12629 0.93624 0.69728 octave:9> v=v/sum(v) v = 0.355753 0.046234 0.342748 0.255265 octave:10> cumsum(v) ans = 0.35575 0.40199 0.74473 1.00000 octave:11> find(a=rand a a = 1 1 1 1 octave:13> a=rand a = 0.15563 octave:14> find(a a=rand a = 0.87758 octave:16> find(a randdisc([1/2,1/2]) error: 'randdisc' undefined near line 1 column 1 octave:17> clear all octave:18> randdisc([1/2,1/2]) error: 'randdisc' undefined near line 1 column 1 octave:18> pwd ans = /home/accounts/personale/clrmrc90/aa1617/sistemi_stocastici octave:19> ls 2016-11-17.txt 2016-12-01.txt invariantPotenze.m randstoc.m 2016-11-17.txt~ epsmet.m metropolis.m randstocsym.m 2016-11-24.txt epsmet.m~ randdisc.m octave:20> addpath(.) parse error: syntax error >>> addpath(.) ^ octave:20> addpath('.') octave:21> randdisc(1) error: 'randdisc' undefined near line 1 column 1 octave:21> quit octave:2> randdisc([1/2,1/2]) ans = 1 octave:3> randdisc([1/2,1/2]) ans = 2 octave:4> randdisc([1/2,1/2]) ans = 1 octave:5> randdisc([1/2,1/2]) ans = 2 octave:6> randdisc([1/2,1/2]) ans = 1 octave:7> randdisc([1/2,1/2]) ans = 1 octave:8> randdisc([1/2,1/2]) ans = 2 octave:9> randdisc([1/2,1/2]) ans = 2 octave:10> randdisc([1/2,1/2]) ans = 2 octave:11> randdisc([1/2,1/2]) ans = 1 octave:12> randdisc(ones(1,6)/6) ans = 6 octave:13> randdisc(ones(1,6)/6) ans = 2 octave:14> randdisc(ones(1,6)/6) ans = 2 octave:15> randdisc(ones(1,6)/6) ans = 5 octave:16> randdisc(ones(1,6)/6) ans = 5 octave:17> randdisc(ones(1,6)/6) ans = 3 octave:18> randdisc(ones(1,6)/6) ans = 4 octave:19> randdisc(ones(1,6)/6) ans = 2 octave:20> rand ans = 0.023559 octave:21> format long e octave:22> ans ans = 2.35594686543770e-02 octave:23> randdisc([1/2,1/2]) ans = 1.00000000000000e+00 octave:24> rand ans = 1.70671419216848e-01 octave:25> rand ans = 2.61695707934039e-01 octave:26> rand ans = 2.64055636554661e-01 octave:27> rand('seed',1) octave:28> rand ans = 8.66802453994751e-01 octave:29> rand ans = 9.12651062011719e-01 octave:30> rand('seed',1) octave:31> rand ans = 8.66802453994751e-01 octave:32> rand ans = 9.12651062011719e-01 octave:33> H(0,4) ans = 4.00000000000000e+00 octave:34> H(1,4) ans = -1.24743230832179e+00 octave:35> H(2,4) ans = 8.45319032669073e-01 octave:36> H(3,4) ans = -2.45013467050290e+01 octave:37> H(4,4) ans = -3.22861611843114e+00 octave:38> H(1,4) ans = -1.24743230832179e+00 octave:39> sa(@(i) H(i,100)) counterstati = 1.89000000000000e+02 ans = -3.53906977328493e+01 octave:40> sa(@(i) H(i,100)) counterstati = 7.10000000000000e+01 ans = -4.00247199211813e+01 octave:41> sa(@(i) H(i,200)) counterstati = 9.80000000000000e+01 ans = -3.27590584632349e+01 octave:42> sa(@(i) H(i,1000)) counterstati = 1.72000000000000e+02 ans = -3.85734038425034e+01 octave:43> sa(@(i) H(i,1000)) counterstati = 1.36000000000000e+02 ans = -2.98786489613260e+01 octave:44> sa(@(i) H(i,10000)) counterstati = 5.60000000000000e+01 ans = -3.52644326886715e+01 octave:45> sa(@(i) H(i,100000)) counterstati = 1.17000000000000e+02 ans = -3.35873089783406e+01 octave:46> quit