octave:3> octave:3> ls octave:4> pwd ans = /home/accounts/altri/caliari/aa0809/sistemi_stocastici octave:5> cd metropolis/ octave:6> ls octave:7> epsilonmetropolis epsilon = 10 ans = Columns 1 through 5: 3.2932e-01 1.0005e-01 3.0084e-02 1.0026e-01 8.3043e-02 1.2805e-01 2.4554e-01 5.1523e-02 1.2286e-01 6.9027e-02 7.6412e-02 1.0225e-01 1.4728e-01 1.2875e-01 1.2294e-01 1.2832e-01 1.2286e-01 6.4875e-02 2.4302e-01 9.1547e-02 2.7736e-02 1.8014e-02 1.6167e-02 2.3891e-02 7.3793e-01 8.2581e-02 9.1796e-02 6.2204e-02 8.2898e-02 1.4048e-01 9.3085e-02 8.4546e-02 9.3604e-02 1.0117e-01 8.8846e-02 1.2144e-01 7.0342e-02 3.4385e-02 8.6766e-02 1.1116e-01 4.7467e-03 3.2498e-03 1.9491e-03 3.6540e-03 1.0848e-02 1.1215e-01 1.2204e-01 5.2315e-02 6.7940e-02 1.1489e-01 Columns 6 through 10: 6.4524e-02 7.8941e-03 9.4884e-02 1.0231e-01 8.7629e-02 9.1796e-02 9.1764e-03 7.0342e-02 8.9647e-02 1.2204e-01 1.2345e-01 2.0162e-02 6.8239e-02 1.0670e-01 1.0382e-01 8.2898e-02 1.0981e-02 8.6766e-02 1.0080e-01 6.7940e-02 3.6661e-02 2.5165e-03 2.9009e-02 7.8090e-02 2.9984e-02 2.6926e-01 1.1606e-02 8.5808e-02 7.9532e-02 9.3835e-02 1.0693e-01 7.4179e-02 7.6435e-02 1.3195e-01 1.4926e-01 8.5808e-02 8.2961e-03 1.9848e-01 1.5487e-01 1.2846e-01 2.8831e-03 5.1918e-04 5.6141e-03 9.6386e-01 2.6763e-03 9.3835e-02 1.6200e-02 1.2846e-01 7.3827e-02 2.1834e-01 guarda P(5,5) epsilon = 0.10000 ans = Columns 1 through 5: 8.1465e-01 2.4638e-12 2.5194e-42 2.4689e-12 8.3043e-02 1.2805e-01 3.0624e-01 1.7522e-31 1.2286e-01 6.9027e-02 7.6412e-02 1.0225e-01 1.6744e-01 1.2875e-01 1.2294e-01 1.2832e-01 1.2286e-01 2.2063e-31 3.1887e-01 9.1547e-02 1.9660e-49 3.1442e-60 9.5965e-90 4.1700e-60 9.2191e-01 8.2581e-02 9.1796e-02 2.1155e-31 8.2898e-02 1.4048e-01 9.3085e-02 8.4546e-02 9.3604e-02 1.0117e-01 8.8846e-02 1.2144e-01 7.0342e-02 1.1694e-31 8.6766e-02 1.1116e-01 4.5488e-135 7.6689e-146 1.5642e-175 8.6227e-146 1.4666e-87 1.1215e-01 1.2204e-01 1.7792e-31 6.7940e-02 1.1489e-01 Columns 6 through 10: 1.5889e-12 6.4716e-109 2.3365e-12 1.0231e-01 2.1578e-12 9.1796e-02 3.0550e-98 7.0342e-02 8.9647e-02 1.2204e-01 1.2345e-01 1.9737e-68 6.8239e-02 1.0670e-01 1.0382e-01 8.2898e-02 3.6557e-98 8.6766e-02 1.0080e-01 6.7940e-02 6.3991e-60 1.4623e-156 5.0634e-60 7.8090e-02 5.2335e-60 3.4307e-01 3.8638e-98 8.5808e-02 7.9532e-02 9.3835e-02 1.0693e-01 7.4179e-02 7.6435e-02 1.3195e-01 1.4926e-01 8.5808e-02 2.7619e-98 2.4116e-01 1.5487e-01 1.2846e-01 6.8037e-146 4.0788e-242 1.3248e-145 1.0000e+00 6.3156e-146 9.3835e-02 5.3932e-98 1.2846e-01 7.3827e-02 2.8686e-01 guarda P(5,5) epsilon = 0.010000 ans = Columns 1 through 7: 0.81465 0.00000 0.00000 0.00000 0.08304 0.00000 0.00000 0.12805 0.30624 0.00000 0.12286 0.06903 0.09180 0.00000 0.07641 0.10225 0.16744 0.12875 0.12294 0.12345 0.00000 0.12832 0.12286 0.00000 0.31887 0.09155 0.08290 0.00000 0.00000 0.00000 0.00000 0.00000 0.92191 0.00000 0.00000 0.08258 0.09180 0.00000 0.08290 0.14048 0.34307 0.00000 0.09309 0.08455 0.09360 0.10117 0.08885 0.10693 0.07418 0.12144 0.07034 0.00000 0.08677 0.11116 0.08581 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.11215 0.12204 0.00000 0.06794 0.11489 0.09383 0.00000 Columns 8 through 10: 0.00000 0.10231 0.00000 0.07034 0.08965 0.12204 0.06824 0.10670 0.10382 0.08677 0.10080 0.06794 0.00000 0.07809 0.00000 0.08581 0.07953 0.09383 0.07644 0.13195 0.14926 0.24116 0.15487 0.12846 0.00000 1.00000 0.00000 0.12846 0.07383 0.28686 guarda P(5,5) octave:8> v=[1 2 3] v = 1 2 3 octave:9> cumsum(v) ans = 1 3 6 octave:10> v=v/sum(v) v = 0.16667 0.33333 0.50000 octave:11> v=[1/2 1/2] v = 0.50000 0.50000 octave:12> ls octave:13> cd .. octave:14> cd simulazione/ octave:15> v v = 0.50000 0.50000 octave:16> cumsum(v) ans = 0.50000 1.00000 octave:17> vargen(v) ans = 1 octave:18> vargen(v) ans = 2 octave:19> vargen(v) ans = 1 octave:20> vargen(v) ans = 1 octave:21> vargen(v) ans = 1 octave:22> vargen(v) ans = 1 octave:23> vargen(v) ans = 2 octave:24> vargen(v) ans = 2 octave:25> vargen(v) ans = 1 octave:26> vargen(v) ans = 1 octave:27> vargen(v) ans = 2 octave:28> vargen(v) ans = 2 octave:29> vargen(v) ans = 1 octave:30> v=[1/3 2/3] v = 0.33333 0.66667 octave:31> for i=1:1000 > y(i)=vargen(v); > end octave:32> y y = Columns 1 through 18: 2 2 2 2 2 1 2 2 2 2 1 2 1 2 2 1 1 2 Columns 19 through 36: 1 2 1 1 1 1 2 1 1 1 2 2 2 1 1 2 2 2 Columns 37 through 54: 1 1 2 2 2 2 1 1 2 1 2 2 1 2 2 2 1 2 Columns 55 through 72: 2 2 2 2 2 1 2 2 2 2 1 2 2 2 1 2 2 2 Columns 73 through 90: 1 2 2 1 2 1 2 2 1 2 1 2 2 1 1 2 2 1 Columns 91 through 108: 2 2 2 2 2 2 1 2 2 1 2 2 2 2 2 2 1 1 Columns 109 through 126: 2 1 1 1 1 2 2 2 1 2 2 1 2 2 2 2 2 2 Columns 127 through 144: 2 2 2 1 1 2 2 2 2 2 2 2 2 1 2 1 2 2 Columns 145 through 162: 1 2 2 2 2 1 2 2 2 2 1 2 2 2 2 1 1 2 Columns 163 through 180: 2 1 2 1 1 1 2 1 1 1 1 1 2 2 2 1 1 2 Columns 181 through 198: 2 2 2 1 2 1 2 1 2 2 1 2 1 2 1 2 2 2 Columns 199 through 216: 1 2 2 1 2 2 1 2 2 2 1 2 1 2 2 2 1 2 Columns 217 through 234: 2 2 2 1 2 1 2 2 1 2 2 1 1 2 1 1 2 2 Columns 235 through 252: 2 2 1 2 1 2 2 1 2 2 2 2 2 2 1 2 2 1 Columns 253 through 270: 1 2 1 1 1 1 2 2 2 2 2 2 2 2 1 2 1 2 Columns 271 through 288: 2 2 2 1 2 2 2 2 1 1 1 2 2 2 2 2 1 2 Columns 289 through 306: 2 1 2 1 2 2 2 2 2 1 2 2 2 1 2 2 1 2 Columns 307 through 324: 2 2 2 2 2 2 1 1 2 2 1 2 1 2 1 2 1 1 Columns 325 through 342: 1 1 1 1 2 1 2 2 1 2 1 2 2 1 2 1 2 2 Columns 343 through 360: 2 2 1 2 2 2 2 2 2 2 2 1 2 2 2 1 2 2 Columns 361 through 378: 2 2 2 2 2 2 2 1 2 2 2 2 2 2 2 2 1 2 Columns 379 through 396: 2 1 1 2 2 2 2 2 1 2 2 2 2 2 2 2 2 1 Columns 397 through 414: 2 1 2 2 2 1 1 2 2 1 1 2 1 1 2 2 1 2 Columns 415 through 432: 2 2 1 2 2 2 2 2 1 1 2 2 1 2 2 2 1 2 Columns 433 through 450: 1 2 1 1 1 1 1 2 1 2 1 2 2 2 2 2 1 2 Columns 451 through 468: 2 1 1 2 2 2 2 2 2 2 2 2 1 1 2 1 1 2 Columns 469 through 486: 1 1 2 1 2 2 2 1 2 2 2 2 2 1 2 2 2 1 Columns 487 through 504: 1 1 2 2 1 2 1 2 1 1 1 1 2 2 2 2 2 2 Columns 505 through 522: 2 1 2 2 2 1 1 2 2 2 2 2 1 2 2 2 2 1 Columns 523 through 540: 1 2 2 2 1 1 1 2 2 1 2 2 1 2 2 2 2 2 Columns 541 through 558: 2 2 1 1 2 2 2 2 2 1 2 2 2 2 1 2 1 2 Columns 559 through 576: 2 2 1 2 2 1 2 2 1 2 1 2 1 2 1 2 1 2 Columns 577 through 594: 1 2 1 1 1 2 2 1 1 2 2 2 2 2 2 2 1 1 Columns 595 through 612: 1 2 1 2 2 1 2 1 2 1 2 2 2 1 2 2 1 1 Columns 613 through 630: 2 1 1 2 1 2 2 1 1 1 2 2 2 1 2 2 2 2 Columns 631 through 648: 1 1 2 1 2 2 1 2 2 1 1 2 2 1 2 1 1 2 Columns 649 through 666: 2 2 2 2 2 2 2 1 1 1 2 2 2 1 2 1 1 2 Columns 667 through 684: 2 1 1 2 1 2 2 2 2 2 2 2 2 2 2 2 2 2 Columns 685 through 702: 2 2 2 2 2 2 2 2 1 1 2 2 1 2 2 2 2 2 Columns 703 through 720: 1 1 2 2 1 2 2 2 2 2 2 1 1 2 2 1 2 2 Columns 721 through 738: 2 2 1 2 1 2 1 1 2 2 2 1 1 2 1 2 2 1 Columns 739 through 756: 2 2 1 2 2 2 2 1 2 1 2 1 2 1 2 2 2 2 Columns 757 through 774: 2 2 2 2 1 2 2 2 2 2 2 1 1 2 2 2 1 2 Columns 775 through 792: 2 2 2 2 2 1 2 1 2 1 1 2 2 1 1 2 1 2 Columns 793 through 810: 1 1 2 2 1 2 2 2 2 1 2 2 2 2 2 1 2 2 Columns 811 through 828: 2 2 1 2 2 1 2 2 2 1 2 1 1 2 2 1 2 2 Columns 829 through 846: 1 2 1 2 2 2 2 2 2 2 2 2 2 2 2 2 1 2 Columns 847 through 864: 2 2 2 2 1 2 2 2 2 1 1 1 1 2 2 1 2 2 Columns 865 through 882: 2 1 2 1 1 2 2 2 1 2 1 1 2 2 1 1 1 2 Columns 883 through 900: 2 2 2 2 2 2 2 2 1 2 1 1 2 2 2 1 1 1 Columns 901 through 918: 2 2 1 2 2 2 2 2 1 1 2 2 2 2 1 2 2 1 Columns 919 through 936: 2 2 2 2 1 2 1 1 1 2 2 1 2 1 2 1 1 2 Columns 937 through 954: 2 2 2 2 2 2 1 1 2 1 2 2 1 1 2 1 2 1 Columns 955 through 972: 1 1 2 2 2 1 1 1 1 2 1 2 2 2 2 1 2 2 Columns 973 through 990: 1 2 2 1 2 2 2 2 2 1 2 1 2 1 2 2 2 2 Columns 991 through 1000: 2 2 2 2 1 2 2 1 1 2 octave:33> find(y==1) ans = Columns 1 through 12: 6 11 13 16 17 19 21 22 23 24 26 27 Columns 13 through 24: 28 32 33 37 38 43 44 46 49 53 60 65 Columns 25 through 36: 69 73 76 78 81 83 86 87 90 97 100 107 Columns 37 through 48: 108 110 111 112 113 117 120 130 131 140 142 145 Columns 49 through 60: 150 155 160 161 164 166 167 168 170 171 172 173 Columns 61 through 72: 174 178 179 184 186 188 191 193 195 199 202 205 Columns 73 through 84: 209 211 215 220 222 225 228 229 231 232 237 239 Columns 85 through 96: 242 249 252 253 255 256 257 258 267 269 274 279 Columns 97 through 108: 280 281 287 290 292 298 302 305 313 314 317 319 Columns 109 through 120: 321 323 324 325 326 327 328 330 333 335 338 340 Columns 121 through 132: 345 354 358 368 377 380 381 387 396 398 402 403 Columns 133 through 144: 406 407 409 410 413 417 423 424 427 431 433 435 Columns 145 through 156: 436 437 438 439 441 443 449 452 453 463 464 466 Columns 157 through 168: 467 469 470 472 476 482 486 487 488 491 493 495 Columns 169 through 180: 496 497 498 506 510 511 517 522 523 527 528 529 Columns 181 through 192: 532 535 543 544 550 555 557 561 564 567 569 571 Columns 193 through 204: 573 575 577 579 580 581 584 585 593 594 595 597 Columns 205 through 216: 600 602 604 608 611 612 614 615 617 620 621 622 Columns 217 through 228: 626 631 632 634 637 640 641 644 646 647 656 657 Columns 229 through 240: 658 662 664 665 668 669 671 693 694 697 703 704 Columns 241 through 252: 707 714 715 718 723 725 727 728 732 733 735 738 Columns 253 through 264: 741 746 748 750 752 761 768 769 773 780 782 784 Columns 265 through 276: 785 788 789 791 793 794 797 802 808 813 816 820 Columns 277 through 288: 822 823 826 829 831 845 851 856 857 858 859 862 Columns 289 through 300: 866 868 869 873 875 876 879 880 881 891 893 894 Columns 301 through 312: 898 899 900 903 909 910 915 918 923 925 926 927 Columns 313 through 324: 930 932 934 935 943 944 946 949 950 952 954 955 Columns 325 through 336: 956 960 961 962 963 965 970 973 976 982 984 986 Columns 337 through 339: 995 998 999 octave:34> length(find(y==1)) ans = 339 octave:35> length(find(y==2)) ans = 661 octave:36> v=rand(1,10); octave:37> v=v/sum(v) v = Columns 1 through 6: 0.041852 0.014427 0.133991 0.121041 0.072111 0.098408 Columns 7 through 10: 0.176337 0.158095 0.121414 0.062324 octave:38> bar(v) octave:39> simulatedannealing Elapsed time is 10.6119 seconds. iminimo = 90 Hminimo = -33.380