octave:2> randstocsym(5) ans = 0.1708817 0.2257497 0.2670181 0.0572432 0.2791073 0.2257497 0.5071308 0.0020589 0.0815106 0.1835499 0.2670181 0.0020589 0.2338201 0.1548995 0.3422033 0.0572432 0.0815106 0.1548995 0.5162181 0.1901287 0.2791073 0.1835499 0.3422033 0.1901287 0.0050108 octave:3> sum(ans,2) ans = 1.00000 1.00000 1.00000 1.00000 1.00000 octave:4> randstocsym(5) ans = 0.1310165 0.2333485 0.2308360 0.2189799 0.1858191 0.2333485 0.0527982 0.3089200 0.3503769 0.0545564 0.2308360 0.3089200 0.1286444 0.1632953 0.1683043 0.2189799 0.3503769 0.1632953 0.0015755 0.2657724 0.1858191 0.0545564 0.1683043 0.2657724 0.3255479 octave:5> randstocsym(5) ans = 0.276338 0.374315 0.010939 0.135126 0.203282 0.374315 0.292798 0.054905 0.116256 0.161726 0.010939 0.054905 0.501041 0.344570 0.088545 0.135126 0.116256 0.344570 0.187600 0.216449 0.203282 0.161726 0.088545 0.216449 0.329997 octave:6> randstocsym(5) ans = 0.144311 0.057467 0.153189 0.305432 0.339601 0.057467 0.259882 0.393787 0.071011 0.217853 0.153189 0.393787 0.256389 0.027885 0.168749 0.305432 0.071011 0.027885 0.099845 0.495827 0.339601 0.217853 0.168749 0.495827 -0.222030 octave:7> randstocsym(5) ans = 0.22745 0.22589 0.13326 0.22687 0.18654 0.22589 0.28225 0.19050 0.12536 0.17601 0.13326 0.19050 0.29687 0.10668 0.27268 0.22687 0.12536 0.10668 0.38512 0.15597 0.18654 0.17601 0.27268 0.15597 0.20880 octave:8> randstocsym(10) ans = Columns 1 through 6: 0.0220823 0.0468183 0.1792240 0.1456906 0.1508865 0.1041417 0.0468183 0.0886300 0.1472745 0.0938051 0.0604299 0.1976627 0.1792240 0.1472745 0.0900226 0.0925768 0.1335440 0.1186326 0.1456906 0.0938051 0.0925768 0.1073825 0.0241177 0.1024324 0.1508865 0.0604299 0.1335440 0.0241177 0.0784054 0.1633600 0.1041417 0.1976627 0.1186326 0.1024324 0.1633600 0.0355424 0.0322584 0.0170627 0.0445252 0.1271772 0.0969637 0.0752508 0.0766382 0.0943342 0.1331789 0.0703596 0.1317747 0.0036444 0.1953693 0.0971673 0.0468959 0.1359451 0.0959453 0.0958771 0.0468907 0.1568153 0.0141254 0.1005131 0.0645729 0.1034561 Columns 7 through 10: 0.0322584 0.0766382 0.1953693 0.0468907 0.0170627 0.0943342 0.0971673 0.1568153 0.0445252 0.1331789 0.0468959 0.0141254 0.1271772 0.0703596 0.1359451 0.1005131 0.0969637 0.1317747 0.0959453 0.0645729 0.0752508 0.0036444 0.0958771 0.1034561 0.0723323 0.0223753 0.2786339 0.2334205 0.0223753 0.1864700 0.0390373 0.2421875 0.2786339 0.0390373 0.0075979 0.0075308 0.2334205 0.2421875 0.0075308 0.0304878 octave:9> octave:9> octave:9> octave:9> octave:9> octave:9> octave:9> octave:9> octave:9> octave:9> octave:9> octave:9> octave:9> octave:9> randstocsym(10) ans = Columns 1 through 5: 1.3490e-01 2.5264e-02 3.6541e-02 1.4631e-01 1.0771e-01 2.5264e-02 1.0423e-01 1.6391e-01 8.2599e-02 8.1776e-02 3.6541e-02 1.6391e-01 3.1474e-02 1.3947e-01 3.2741e-02 1.4631e-01 8.2599e-02 1.3947e-01 2.7751e-01 1.3600e-02 1.0771e-01 8.1776e-02 3.2741e-02 1.3600e-02 7.5114e-02 2.0791e-01 1.5298e-01 1.9038e-01 1.7759e-01 3.4928e-02 2.3539e-02 6.7890e-02 1.0505e-01 2.0724e-02 1.2519e-04 1.0569e-01 3.4504e-02 2.0652e-01 7.3928e-03 3.3659e-01 1.2579e-01 1.3892e-01 6.0356e-02 1.2317e-01 1.2633e-01 8.6337e-02 1.4793e-01 3.3551e-02 1.1621e-02 1.9108e-01 Columns 6 through 10: 2.0791e-01 2.3539e-02 1.0569e-01 1.2579e-01 8.6337e-02 1.5298e-01 6.7890e-02 3.4504e-02 1.3892e-01 1.4793e-01 1.9038e-01 1.0505e-01 2.0652e-01 6.0356e-02 3.3551e-02 1.7759e-01 2.0724e-02 7.3928e-03 1.2317e-01 1.1621e-02 3.4928e-02 1.2519e-04 3.3659e-01 1.2633e-01 1.9108e-01 7.5535e-02 2.9257e-02 3.7009e-02 7.3160e-03 8.7092e-02 2.9257e-02 1.8305e-01 1.7794e-01 1.7679e-01 2.1564e-01 3.7009e-02 1.7794e-01 5.6434e-02 2.8065e-02 9.8533e-03 7.3160e-03 1.7679e-01 2.8065e-02 6.4608e-02 1.4865e-01 8.7092e-02 2.1564e-01 9.8533e-03 1.4865e-01 6.8251e-02 octave:10> octave:10> octave:10> octave:10> octave:10> octave:10> octave:10> octave:10> octave:10> octave:10> randstocsym(10) ans = Columns 1 through 6: 0.0639439 0.0060146 0.2016976 0.0447752 0.1913612 0.0980120 0.0060146 0.0892980 0.0298701 0.1353567 0.0395566 0.1183474 0.2016976 0.0298701 0.0727607 0.1530070 0.1435932 0.1250573 0.0447752 0.1353567 0.1530070 0.0153006 0.1577987 0.1651250 0.1913612 0.0395566 0.1435932 0.1577987 0.0124816 0.1718946 0.0980120 0.1183474 0.1250573 0.1651250 0.1718946 0.1295979 0.0831825 0.1742468 0.0130291 0.0805825 0.0515382 0.0591633 0.0292217 0.1008451 0.0194263 0.1159117 0.0719372 0.0656357 0.1494212 0.1361423 0.1463888 0.0118150 0.0065289 0.0448517 0.1323703 0.1703222 0.0951700 0.1203277 0.1533098 0.0223150 Columns 7 through 10: 0.0831825 0.0292217 0.1494212 0.1323703 0.1742468 0.1008451 0.1361423 0.1703222 0.0130291 0.0194263 0.1463888 0.0951700 0.0805825 0.1159117 0.0118150 0.1203277 0.0515382 0.0719372 0.0065289 0.1533098 0.0591633 0.0656357 0.0448517 0.0223150 0.1858146 0.1660169 0.1760443 0.0103818 0.1660169 0.2031350 0.2199013 0.0079692 0.1760443 0.2199013 0.0895284 0.0193782 0.0103818 0.0079692 0.0193782 0.2684558 octave:11> octave:11> octave:11> octave:11> octave:11> octave:11> octave:11> octave:11> octave:11> octave:11> randstocsym(5) ans = 0.048818 0.241391 0.207336 0.076150 0.426306 0.241391 0.265934 0.339543 0.073517 0.079614 0.207336 0.339543 0.116466 0.243337 0.093318 0.076150 0.073517 0.243337 0.235240 0.371755 0.426306 0.079614 0.093318 0.371755 0.029006 octave:12> randstocsym(5) ans = 0.031125 0.172184 0.114881 0.354322 0.327488 0.172184 0.232578 0.251365 0.119920 0.223953 0.114881 0.251365 0.254518 0.306975 0.072261 0.354322 0.119920 0.306975 0.033291 0.185491 0.327488 0.223953 0.072261 0.185491 0.190807 octave:13> randstocsym(5) ans = 0.042627 0.154682 0.228817 0.289161 0.284712 0.154682 0.309369 0.183526 0.185096 0.167327 0.228817 0.183526 0.102157 0.194027 0.291473 0.289161 0.185096 0.194027 0.030760 0.300956 0.284712 0.167327 0.291473 0.300956 -0.044468 octave:14> any(ans)<0 ans = 0 0 0 0 0 octave:15> P=randstocsym(5) P = 0.281886 0.229887 0.245657 0.075628 0.166942 0.229887 0.215141 0.307067 0.088581 0.159324 0.245657 0.307067 0.063401 0.125788 0.258086 0.075628 0.088581 0.125788 0.066898 0.643105 0.166942 0.159324 0.258086 0.643105 -0.227458 octave:16> any(P) ans = 1 1 1 1 1 octave:17> any(P(:)<0) ans = 1 octave:18> P(:)<0 ans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 octave:19> octave:19> octave:19> octave:19> octave:19> octave:19> octave:19> octave:19> octave:19> octave:19> v=[1,2,3,4,5,6] v = 1 2 3 4 5 6 octave:20> any(v==3) ans = 1 octave:21> v==3 ans = 0 0 1 0 0 0 octave:22> find(v==3) ans = 3 octave:23> [r,idx]=find(v==3) r = 1 idx = 3 octave:24> find(v==3) ans = 3 octave:25> P=randstocsym(5) P = 0.242098 0.141221 0.184338 0.168991 0.263353 0.141221 0.220237 0.230446 0.198448 0.209648 0.184338 0.230446 0.195603 0.181137 0.208476 0.168991 0.198448 0.181137 0.191920 0.259505 0.263353 0.209648 0.208476 0.259505 0.059018 octave:26> P=randstocsym(5) P = 0.157682 0.076334 0.381679 0.237484 0.146820 0.076334 0.114026 0.258643 0.517724 0.033273 0.381679 0.258643 0.170438 0.111919 0.077321 0.237484 0.517724 0.111919 0.065671 0.067202 0.146820 0.033273 0.077321 0.067202 0.675384 octave:27> P=randstocsym(5) P = 0.127804 0.166384 0.048045 0.280515 0.377251 0.166384 0.193608 0.251495 0.267092 0.121422 0.048045 0.251495 0.275492 0.204783 0.220184 0.280515 0.267092 0.204783 0.059806 0.187804 0.377251 0.121422 0.220184 0.187804 0.093340 octave:28> P=randstocsym(5) P = 0.237808 0.202625 0.232784 0.175721 0.151062 0.202625 0.123758 0.188580 0.233283 0.251754 0.232784 0.188580 0.168839 0.121936 0.287861 0.175721 0.233283 0.121936 0.220583 0.248476 0.151062 0.251754 0.287861 0.248476 0.060847 octave:29> P=randstocsym(5) P = 0.162924 0.186851 0.015284 0.530477 0.104464 0.186851 0.262460 0.220594 0.190163 0.139932 0.015284 0.220594 0.067512 0.110537 0.586073 0.530477 0.190163 0.110537 0.136242 0.032581 0.104464 0.139932 0.586073 0.032581 0.136950 octave:30> sum(P,2) ans = 1.00000 1.00000 1.00000 1.00000 1.00000 octave:31> v=rand(1,4) v = 0.3156068 0.5224611 0.0074561 0.0029091 octave:32> V=repmat(v,4,1) V = 0.3156068 0.5224611 0.0074561 0.0029091 0.3156068 0.5224611 0.0074561 0.0029091 0.3156068 0.5224611 0.0074561 0.0029091 0.3156068 0.5224611 0.0074561 0.0029091 octave:33> W=V' W = 0.3156068 0.3156068 0.3156068 0.3156068 0.5224611 0.5224611 0.5224611 0.5224611 0.0074561 0.0074561 0.0074561 0.0074561 0.0029091 0.0029091 0.0029091 0.0029091 octave:34> repmat(v',1,4) ans = 0.3156068 0.3156068 0.3156068 0.3156068 0.5224611 0.5224611 0.5224611 0.5224611 0.0074561 0.0074561 0.0074561 0.0074561 0.0029091 0.0029091 0.0029091 0.0029091 octave:35> V>=W ans = 1 1 0 0 0 1 0 0 1 1 1 0 1 1 1 1 octave:36> Q=randstocsym(4) Q = 0.27284 0.32300 0.14345 0.26071 0.32300 0.17956 0.26176 0.23568 0.14345 0.26176 0.29883 0.29595 0.26071 0.23568 0.29595 0.20766 octave:37> Q(V>=W) ans = 0.27284 0.14345 0.26071 0.32300 0.17956 0.26176 0.23568 0.29883 0.29595 0.20766 octave:38> sum(Q) ans = 1.0000 1.0000 1.0000 1.0000 octave:39> P=zeros(4) P = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 octave:40> P(V>=W)=Q(V>=W) P = 0.27284 0.32300 0.00000 0.00000 0.00000 0.17956 0.00000 0.00000 0.14345 0.26176 0.29883 0.00000 0.26071 0.23568 0.29595 0.20766 octave:41> Q Q = 0.27284 0.32300 0.14345 0.26071 0.32300 0.17956 0.26176 0.23568 0.14345 0.26176 0.29883 0.29595 0.26071 0.23568 0.29595 0.20766 octave:42> R(V>=W)=Q(V>=W) R = Columns 1 through 7: 0.27284 0.00000 0.14345 0.26071 0.32300 0.17956 0.26176 Columns 8 through 14: 0.23568 0.00000 0.00000 0.29883 0.29595 0.00000 0.00000 Columns 15 and 16: 0.00000 0.20766 octave:43> Q(V>=W) ans = 0.27284 0.14345 0.26071 0.32300 0.17956 0.26176 0.23568 0.29883 0.29595 0.20766 octave:44> V V = 0.3156068 0.5224611 0.0074561 0.0029091 0.3156068 0.5224611 0.0074561 0.0029091 0.3156068 0.5224611 0.0074561 0.0029091 0.3156068 0.5224611 0.0074561 0.0029091 octave:45> W W = 0.3156068 0.3156068 0.3156068 0.3156068 0.5224611 0.5224611 0.5224611 0.5224611 0.0074561 0.0074561 0.0074561 0.0074561 0.0029091 0.0029091 0.0029091 0.0029091 octave:46> V>=W ans = 1 1 0 0 0 1 0 0 1 1 1 0 1 1 1 1 octave:47> (V>=W) - eye(4) ans = 0 1 0 0 0 0 0 0 1 1 0 0 1 1 1 0 octave:48> P(ans) error: subscript indices must be either positive integers less than 2^31 or logicals octave:48> (V>=W) - logical(4) ans = 0 0 -1 -1 -1 0 -1 -1 0 0 0 -1 0 0 0 0 octave:49> logical(4) ans = 1 octave:50> !eye(4) ans = 0 1 1 1 1 0 1 1 1 1 0 1 1 1 1 0 octave:51> whos Variables in the current scope: Attr Name Size Bytes Class ==== ==== ==== ===== ===== P 4x4 128 double Q 4x4 128 double R 1x16 128 double V 4x4 128 double W 4x4 128 double ans 4x4 16 logical idx 1x1 8 double r 1x1 8 double v 1x4 32 double Total is 102 elements using 704 bytes octave:52> I=eye(4) I = Diagonal Matrix 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 octave:53> whos Variables in the current scope: Attr Name Size Bytes Class ==== ==== ==== ===== ===== I 4x4 32 double P 4x4 128 double Q 4x4 128 double R 1x16 128 double V 4x4 128 double W 4x4 128 double ans 4x4 16 logical idx 1x1 8 double r 1x1 8 double v 1x4 32 double Total is 118 elements using 736 bytes octave:54> I=!eye(4) I = 0 1 1 1 1 0 1 1 1 1 0 1 1 1 1 0 octave:55> whos Variables in the current scope: Attr Name Size Bytes Class ==== ==== ==== ===== ===== I 4x4 16 logical P 4x4 128 double Q 4x4 128 double R 1x16 128 double V 4x4 128 double W 4x4 128 double ans 4x4 16 logical idx 1x1 8 double r 1x1 8 double v 1x4 32 double Total is 118 elements using 720 bytes octave:56> clear al octave:57> clear all octave:58> I=eye(4) I = Diagonal Matrix 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 octave:59> whos Variables in the current scope: Attr Name Size Bytes Class ==== ==== ==== ===== ===== I 4x4 32 double Total is 16 elements using 32 bytes octave:60> I=!eye(4) I = 0 1 1 1 1 0 1 1 1 1 0 1 1 1 1 0 octave:61> whos Variables in the current scope: Attr Name Size Bytes Class ==== ==== ==== ===== ===== I 4x4 16 logical Total is 16 elements using 16 bytes octave:62> logical(eye(4)) ans = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 octave:63> whos Variables in the current scope: Attr Name Size Bytes Class ==== ==== ==== ===== ===== I 4x4 16 logical ans 4x4 16 logical Total is 32 elements using 32 bytes octave:64> (V>=W)-logical(eye(4)) error: 'V' undefined near line 1 column 2 octave:64> v=rand(1,4) v = 0.010823 0.450081 0.013137 0.780784 octave:65> V=repmat(v,4,1) V = 0.010823 0.450081 0.013137 0.780784 0.010823 0.450081 0.013137 0.780784 0.010823 0.450081 0.013137 0.780784 0.010823 0.450081 0.013137 0.780784 octave:66> W=V' W = 0.010823 0.010823 0.010823 0.010823 0.450081 0.450081 0.450081 0.450081 0.013137 0.013137 0.013137 0.013137 0.780784 0.780784 0.780784 0.780784 octave:67> (V>=W)-logical(eye(4)) ans = 0 1 1 1 0 0 0 1 0 1 0 1 0 0 0 0 octave:68> whos Variables in the current scope: Attr Name Size Bytes Class ==== ==== ==== ===== ===== I 4x4 16 logical V 4x4 128 double W 4x4 128 double ans 4x4 128 double v 1x4 32 double Total is 68 elements using 432 bytes octave:69> (V-diag(diag(V))>=W) ans = 0 1 1 1 0 0 0 1 0 1 0 1 0 0 0 0 octave:70> V>=W ans = 1 1 1 1 0 1 0 1 0 1 1 1 0 0 0 1 octave:71> v=[1,0,1,1] v = 1 0 1 1 octave:72> V=repmat(v,4,1) V = 1 0 1 1 1 0 1 1 1 0 1 1 1 0 1 1 octave:73> W=V' W = 1 1 1 1 0 0 0 0 1 1 1 1 1 1 1 1 octave:74> V>=W ans = 1 0 1 1 1 1 1 1 1 0 1 1 1 0 1 1 octave:75> (V-diag(diag(V))>=W) ans = 0 0 1 1 1 1 1 1 1 0 0 1 1 0 1 0 octave:76> (V>=W) & !eye(4) ans = 0 0 1 1 1 0 1 1 1 0 0 1 1 0 1 0 octave:77> V>=W ans = 1 0 1 1 1 1 1 1 1 0 1 1 1 0 1 1 octave:78> (V>=W) & !eye(4) ans = 0 0 1 1 1 0 1 1 1 0 0 1 1 0 1 0 octave:79> (V>=W)&!eye(4) ans = 0 0 1 1 1 0 1 1 1 0 0 1 1 0 1 0 octave:80> whos Variables in the current scope: Attr Name Size Bytes Class ==== ==== ==== ===== ===== I 4x4 16 logical V 4x4 128 double W 4x4 128 double ans 4x4 16 logical v 1x4 32 double Total is 68 elements using 320 bytes octave:81> !eye(4) ans = 0 1 1 1 1 0 1 1 1 1 0 1 1 1 1 0 octave:82> (V>=W) & !eye(4) ans = 0 0 1 1 1 0 1 1 1 0 0 1 1 0 1 0 octave:83> v v = 1 0 1 1 octave:84> v=rand(1,4) v = 0.96777 0.13399 0.66388 0.33515 octave:85> P=metropolis(v) P = 0.6092276 0.0198951 0.3652923 0.0055850 0.1436993 0.3189833 0.1713833 0.3659341 0.5325017 0.0345892 0.3631964 0.0697127 0.0161271 0.1462946 0.1380908 0.6994875 octave:86> invariantPotenze(P,1e-6,100) ans = 0.460669 0.063780 0.316016 0.159536 octave:87> v/sum(v) ans = 0.460670 0.063779 0.316016 0.159535 octave:88> invariantPotenze(P,1e-8,100) ans = 0.460670 0.063779 0.316016 0.159535 octave:89> exp(-1000) ans = 0 octave:90> exp(-100) ans = 3.7201e-44 octave:91> exp(-500) ans = 7.1246e-218 octave:92> exp(-800) ans = 0 octave:93> exp(-700) ans = 9.8597e-305 octave:94> realmin ans = 2.2251e-308 octave:95> format long e octave:96> realmin ans = 2.22507385850720e-308 octave:97> 2^-1024 ans = 5.56268464626800e-309 octave:98> 2^-1023 ans = 1.11253692925360e-308 octave:99> realmax ans = 1.79769313486232e+308 octave:100> H(1)=-1 H = -1.00000000000000e+00 octave:101> H(2)=1 H = -1.00000000000000e+00 1.00000000000000e+00 octave:102> epsilon=1e-3 epsilon = 1.00000000000000e-03 octave:103> v=exp(-H/epsilon) v = Inf 0.00000000000000e+00 octave:104> v = 1./sum(exp(-(repmat(H',1,2)-repmat(H,2,1))/epsilon)) v = 1.00000000000000e+00 0.00000000000000e+00 octave:105> H=[0,-1,1] H = 0.00000000000000e+00 -1.00000000000000e+00 1.00000000000000e+00 octave:106> epsilon=1 epsilon = 1.00000000000000e+00 octave:107> v = 1./sum(exp(-(repmat(H',1,3)-repmat(H,3,1))/epsilon)) v = 2.44728471054798e-01 6.65240955774822e-01 9.00305731703805e-02 octave:108> format octave:109> v = 1./sum(exp(-(repmat(H',1,3)-repmat(H,3,1))/epsilon)) v = 0.244728 0.665241 0.090031 octave:110> epsilon=1e-1 epsilon = 0.10000 octave:111> v = 1./sum(exp(-(repmat(H',1,3)-repmat(H,3,1))/epsilon)) v = 4.5398e-05 9.9995e-01 2.0611e-09 octave:112> epsilon=1e-2 epsilon = 0.010000 octave:113> v = 1./sum(exp(-(repmat(H',1,3)-repmat(H,3,1))/epsilon)) v = 3.7201e-44 1.0000e+00 1.3839e-87 octave:114> H H = 0 -1 1 octave:115> epsmet ans = Columns 1 through 5: 2.6616e-01 8.8394e-02 6.2747e-02 8.8394e-02 1.0000e-01 1.0000e-01 1.9607e-01 7.0985e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.5359e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 7.0985e-02 1.9607e-01 1.0000e-01 5.7792e-02 5.1085e-02 3.6263e-02 5.1085e-02 5.3369e-01 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 2.1540e-02 1.9040e-02 1.3515e-02 1.9040e-02 3.7271e-02 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 Columns 6 through 10: 8.8394e-02 2.9121e-02 8.8394e-02 1.0000e-01 8.8394e-02 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 4.6411e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 5.1085e-02 1.6830e-02 5.1085e-02 1.0000e-01 5.1085e-02 1.9607e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 3.2945e-02 1.9607e-01 1.0000e-01 1.0000e-01 1.9040e-02 6.2726e-03 1.9040e-02 8.2620e-01 1.9040e-02 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.9607e-01 guarda P(5,5) e P(9,9) ans = Columns 1 through 5: 3.6148e-01 7.8134e-02 3.9371e-02 7.8134e-02 1.0000e-01 1.0000e-01 2.3876e-01 5.0389e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.7846e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 5.0389e-02 2.3876e-01 1.0000e-01 3.3400e-02 2.6097e-02 1.3150e-02 2.6097e-02 7.2013e-01 1.0000e-01 1.0000e-01 5.0389e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 5.0389e-02 1.0000e-01 1.0000e-01 4.6396e-03 3.6251e-03 1.8267e-03 3.6251e-03 1.3891e-02 1.0000e-01 1.0000e-01 5.0389e-02 1.0000e-01 1.0000e-01 Columns 6 through 10: 7.8134e-02 8.4805e-03 7.8134e-02 1.0000e-01 7.8134e-02 1.0000e-01 1.0854e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.1540e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0854e-02 1.0000e-01 1.0000e-01 1.0000e-01 2.6097e-02 2.8325e-03 2.6097e-02 1.0000e-01 2.6097e-02 2.3876e-01 1.0854e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0854e-02 2.3876e-01 1.0000e-01 1.0000e-01 3.6251e-03 3.9346e-04 3.6251e-03 9.6112e-01 3.6251e-03 1.0000e-01 1.0854e-02 1.0000e-01 1.0000e-01 2.3876e-01 guarda P(5,5) e P(9,9) ans = Columns 1 through 5: 7.5759e-01 8.4805e-03 8.9498e-06 8.4805e-03 1.0000e-01 1.0000e-01 2.9989e-01 1.0553e-04 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0553e-04 2.9989e-01 1.0000e-01 1.7275e-06 1.4650e-07 1.5461e-10 1.4650e-07 9.0000e-01 1.0000e-01 1.0000e-01 1.0553e-04 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0553e-04 1.0000e-01 1.0000e-01 4.6217e-15 3.9194e-16 4.1363e-19 3.9194e-16 2.6753e-10 1.0000e-01 1.0000e-01 1.0553e-04 1.0000e-01 1.0000e-01 Columns 6 through 10: 8.4805e-03 1.9240e-12 8.4805e-03 1.0000e-01 8.4805e-03 1.0000e-01 2.2688e-11 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.1498e-08 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.2688e-11 1.0000e-01 1.0000e-01 1.0000e-01 1.4650e-07 3.3238e-17 1.4650e-07 1.0000e-01 1.4650e-07 2.9989e-01 2.2688e-11 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.2688e-11 2.9989e-01 1.0000e-01 1.0000e-01 3.9194e-16 8.8923e-26 3.9194e-16 1.0000e+00 3.9194e-16 1.0000e-01 2.2688e-11 1.0000e-01 1.0000e-01 2.9989e-01 guarda P(5,5) e P(9,9) ans = Columns 1 through 5: 8.0000e-01 1.9240e-12 3.2971e-42 1.9240e-12 1.0000e-01 1.0000e-01 3.0000e-01 1.7137e-31 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.7137e-31 3.0000e-01 1.0000e-01 2.3674e-49 4.5551e-60 7.8058e-90 4.5551e-60 9.0000e-01 1.0000e-01 1.0000e-01 1.7137e-31 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.7137e-31 1.0000e-01 1.0000e-01 4.4462e-135 8.5546e-146 1.4660e-175 8.5546e-146 1.8780e-87 1.0000e-01 1.0000e-01 1.7137e-31 1.0000e-01 1.0000e-01 Columns 6 through 10: 1.9240e-12 6.9524e-109 1.9240e-12 1.0000e-01 1.9240e-12 1.0000e-01 3.6134e-98 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.1086e-68 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 3.6134e-98 1.0000e-01 1.0000e-01 1.0000e-01 4.5551e-60 1.6459e-156 4.5551e-60 1.0000e-01 4.5551e-60 3.0000e-01 3.6134e-98 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 3.6134e-98 3.0000e-01 1.0000e-01 1.0000e-01 8.5546e-146 3.0911e-242 8.5546e-146 1.0000e+00 8.5546e-146 1.0000e-01 3.6134e-98 1.0000e-01 1.0000e-01 3.0000e-01 guarda P(5,5) e P(9,9) ans = Columns 1 through 5: 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 Columns 6 through 10: 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 0.0000e+00 0.0000e+00 0.0000e+00 1.0000e+00 0.0000e+00 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 guarda P(5,5) e P(9,9) octave:116> octave:116> octave:116> epsmet ans = Columns 1 through 5: 2.6616e-01 8.8394e-02 6.2747e-02 8.8394e-02 1.0000e-01 1.0000e-01 1.9607e-01 7.0985e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.5359e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 7.0985e-02 1.9607e-01 1.0000e-01 5.7792e-02 5.1085e-02 3.6263e-02 5.1085e-02 5.3369e-01 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 2.1540e-02 1.9040e-02 1.3515e-02 1.9040e-02 3.7271e-02 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 Columns 6 through 10: 8.8394e-02 2.9121e-02 8.8394e-02 1.0000e-01 8.8394e-02 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 4.6411e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 5.1085e-02 1.6830e-02 5.1085e-02 1.0000e-01 5.1085e-02 1.9607e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 3.2945e-02 1.9607e-01 1.0000e-01 1.0000e-01 1.9040e-02 6.2726e-03 1.9040e-02 8.2620e-01 1.9040e-02 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.9607e-01 guarda P(5,5) e P(9,9) ans = Columns 1 through 5: 3.6148e-01 7.8134e-02 3.9371e-02 7.8134e-02 1.0000e-01 1.0000e-01 2.3876e-01 5.0389e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.7846e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 5.0389e-02 2.3876e-01 1.0000e-01 3.3400e-02 2.6097e-02 1.3150e-02 2.6097e-02 7.2013e-01 1.0000e-01 1.0000e-01 5.0389e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 5.0389e-02 1.0000e-01 1.0000e-01 4.6396e-03 3.6251e-03 1.8267e-03 3.6251e-03 1.3891e-02 1.0000e-01 1.0000e-01 5.0389e-02 1.0000e-01 1.0000e-01 Columns 6 through 10: 7.8134e-02 8.4805e-03 7.8134e-02 1.0000e-01 7.8134e-02 1.0000e-01 1.0854e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.1540e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0854e-02 1.0000e-01 1.0000e-01 1.0000e-01 2.6097e-02 2.8325e-03 2.6097e-02 1.0000e-01 2.6097e-02 2.3876e-01 1.0854e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0854e-02 2.3876e-01 1.0000e-01 1.0000e-01 3.6251e-03 3.9346e-04 3.6251e-03 9.6112e-01 3.6251e-03 1.0000e-01 1.0854e-02 1.0000e-01 1.0000e-01 2.3876e-01 guarda P(5,5) e P(9,9) ans = Columns 1 through 5: 7.5759e-01 8.4805e-03 8.9498e-06 8.4805e-03 1.0000e-01 1.0000e-01 2.9989e-01 1.0553e-04 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0553e-04 2.9989e-01 1.0000e-01 1.7275e-06 1.4650e-07 1.5461e-10 1.4650e-07 9.0000e-01 1.0000e-01 1.0000e-01 1.0553e-04 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0553e-04 1.0000e-01 1.0000e-01 4.6217e-15 3.9194e-16 4.1363e-19 3.9194e-16 2.6753e-10 1.0000e-01 1.0000e-01 1.0553e-04 1.0000e-01 1.0000e-01 Columns 6 through 10: 8.4805e-03 1.9240e-12 8.4805e-03 1.0000e-01 8.4805e-03 1.0000e-01 2.2688e-11 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.1498e-08 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.2688e-11 1.0000e-01 1.0000e-01 1.0000e-01 1.4650e-07 3.3238e-17 1.4650e-07 1.0000e-01 1.4650e-07 2.9989e-01 2.2688e-11 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.2688e-11 2.9989e-01 1.0000e-01 1.0000e-01 3.9194e-16 8.8923e-26 3.9194e-16 1.0000e+00 3.9194e-16 1.0000e-01 2.2688e-11 1.0000e-01 1.0000e-01 2.9989e-01 guarda P(5,5) e P(9,9) octave:116> octave:116> octave:116> diary off