octave:2> randstocsym(5) ans = 0.468345 0.170304 0.069024 0.237861 0.054467 0.170304 0.245893 0.245005 0.151076 0.187722 0.069024 0.245005 0.345254 0.284389 0.056329 0.237861 0.151076 0.284389 0.258787 0.067888 0.054467 0.187722 0.056329 0.067888 0.633594 octave:3> Q=randstocsym(5),sum(Q),sum(Q) Q = 0.222985 0.090601 0.261966 0.315183 0.109265 0.090601 0.052445 0.206398 0.345829 0.304726 0.261966 0.206398 0.360696 0.033751 0.137189 0.315183 0.345829 0.033751 0.244434 0.060804 0.109265 0.304726 0.137189 0.060804 0.388016 ans = 1 1 1 1 1 ans = 1 1 1 1 1 octave:4> Q=randstocsym(5),sum(Q),sum(Q) Q = 0.1960887 0.2732016 0.2237157 0.0466381 0.2603559 0.2732016 0.1730529 0.1771933 0.3691955 0.0073567 0.2237157 0.1771933 0.1715016 0.1773107 0.2502787 0.0466381 0.3691955 0.1773107 0.0608418 0.3460139 0.2603559 0.0073567 0.2502787 0.3460139 0.1359948 ans = 1.00000 1.00000 1.00000 1.00000 1.00000 ans = 1.00000 1.00000 1.00000 1.00000 1.00000 octave:5> Q=randstocsym(5),sum(Q),sum(Q) Q = 0.2733644 0.1048214 0.2066349 0.1649334 0.2502459 0.1048214 0.0096099 0.3208609 0.3468579 0.2178498 0.2066349 0.3208609 0.0163409 0.3003727 0.1557906 0.1649334 0.3468579 0.3003727 0.0675995 0.1202364 0.2502459 0.2178498 0.1557906 0.1202364 0.2558772 ans = 1.00000 1.00000 1.00000 1.00000 1.00000 ans = 1.00000 1.00000 1.00000 1.00000 1.00000 octave:6> Q=randstocsym(5),sum(Q),sum(Q) Q = 0.093775 0.087402 0.519329 0.064065 0.235428 0.087402 0.199060 0.196439 0.283780 0.233319 0.519329 0.196439 0.014230 0.189407 0.080595 0.064065 0.283780 0.189407 0.203819 0.258928 0.235428 0.233319 0.080595 0.258928 0.191730 ans = 1 1 1 1 1 ans = 1 1 1 1 1 octave:7> Q=randstocsym(5),sum(Q),sum(Q) Q = 0.024416 0.358123 0.169098 0.385263 0.063100 0.358123 0.172198 0.044175 0.050233 0.375272 0.169098 0.044175 0.347188 0.351920 0.087619 0.385263 0.050233 0.351920 0.159992 0.052592 0.063100 0.375272 0.087619 0.052592 0.421416 ans = 1.0000 1.0000 1.0000 1.0000 1.0000 ans = 1.0000 1.0000 1.0000 1.0000 1.0000 octave:8> Q=randstocsym(5),sum(Q),sum(Q) Q = 0.357621 0.344190 0.022488 0.157773 0.117928 0.344190 0.040470 0.251733 0.257423 0.106184 0.022488 0.251733 0.392597 0.218609 0.114573 0.157773 0.257423 0.218609 0.109176 0.257020 0.117928 0.106184 0.114573 0.257020 0.404295 ans = 1.00000 1.00000 1.00000 1.00000 1.00000 ans = 1.00000 1.00000 1.00000 1.00000 1.00000 octave:9> Q=randstocsym(5),sum(Q),sum(Q) Q = 0.394640 0.151145 0.048322 0.361546 0.044347 0.151145 0.206288 0.026103 0.252914 0.363550 0.048322 0.026103 0.063054 0.450989 0.411532 0.361546 0.252914 0.450989 -0.041929 -0.023520 0.044347 0.363550 0.411532 -0.023520 0.204091 ans = 1.00000 1.00000 1.00000 1.00000 1.00000 ans = 1.00000 1.00000 1.00000 1.00000 1.00000 octave:10> sum(Q(4,1:3)) ans = 1.0654 octave:11> Q Q = 0.394640 0.151145 0.048322 0.361546 0.044347 0.151145 0.206288 0.026103 0.252914 0.363550 0.048322 0.026103 0.063054 0.450989 0.411532 0.361546 0.252914 0.450989 -0.041929 -0.023520 0.044347 0.363550 0.411532 -0.023520 0.204091 octave:12> Q=Q+1 Q = 1.39464 1.15115 1.04832 1.36155 1.04435 1.15115 1.20629 1.02610 1.25291 1.36355 1.04832 1.02610 1.06305 1.45099 1.41153 1.36155 1.25291 1.45099 0.95807 0.97648 1.04435 1.36355 1.41153 0.97648 1.20409 octave:13> sum(Q) ans = 6 6 6 6 6 octave:14> sum(Q,2) ans = 6 6 6 6 6 octave:15> Q=Q/6 Q = 0.23244 0.19186 0.17472 0.22692 0.17406 0.19186 0.20105 0.17102 0.20882 0.22726 0.17472 0.17102 0.17718 0.24183 0.23526 0.22692 0.20882 0.24183 0.15968 0.16275 0.17406 0.22726 0.23526 0.16275 0.20068 octave:16> sum(Q) ans = 1.00000 1.00000 1.00000 1.00000 1.00000 octave:17> Q=randstocsym(5),sum(Q),sum(Q) Q = 0.187442 0.206864 0.194855 0.180698 0.230141 0.206864 0.222902 0.165523 0.222613 0.182099 0.194855 0.165523 0.186178 0.161280 0.292163 0.180698 0.222613 0.161280 0.216520 0.218890 0.230141 0.182099 0.292163 0.218890 0.076707 ans = 1 1 1 1 1 ans = 1 1 1 1 1 octave:18> Q=randstocsym(5),sum(Q),sum(Q) Q = 0.239705 0.130706 0.225122 0.137670 0.266797 0.130706 0.212919 0.264806 0.209132 0.182437 0.225122 0.264806 0.169102 0.156037 0.184932 0.137670 0.209132 0.156037 0.171990 0.325171 0.266797 0.182437 0.184932 0.325171 0.040662 ans = 1.00000 1.00000 1.00000 1.00000 1.00000 ans = 1.00000 1.00000 1.00000 1.00000 1.00000 octave:19> Q=randstocsym(5),sum(Q),sum(Q) Q = 0.259023 0.202303 0.136041 0.230492 0.172140 0.202303 0.280489 0.079757 0.243594 0.193857 0.136041 0.079757 0.161552 0.284560 0.338091 0.230492 0.243594 0.284560 0.128102 0.113252 0.172140 0.193857 0.338091 0.113252 0.182660 ans = 1.0000 1.0000 1.0000 1.0000 1.0000 ans = 1.0000 1.0000 1.0000 1.0000 1.0000 octave:20> Q=randstocsym(5),sum(Q),sum(Q) Q = 0.290236 0.387598 0.060551 0.197767 0.063848 0.387598 0.130745 0.066107 0.196679 0.218871 0.060551 0.066107 0.086855 0.360536 0.425951 0.197767 0.196679 0.360536 0.154547 0.090470 0.063848 0.218871 0.425951 0.090470 0.200858 ans = 1.00000 1.00000 1.00000 1.00000 1.00000 ans = 1.00000 1.00000 1.00000 1.00000 1.00000 octave:21> rand('seed',1) octave:22> rand ans = 0.86680 octave:23> rand ans = 0.91265 octave:24> rand('seed',1) octave:25> rand ans = 0.86680 octave:26> v=rand(1,4) v = 0.912651 0.093661 0.166461 0.740808 octave:27> V=repmat(v,4,1) V = 0.912651 0.093661 0.166461 0.740808 0.912651 0.093661 0.166461 0.740808 0.912651 0.093661 0.166461 0.740808 0.912651 0.093661 0.166461 0.740808 octave:28> W=repmat(v',1,4) W = 0.912651 0.912651 0.912651 0.912651 0.093661 0.093661 0.093661 0.093661 0.166461 0.166461 0.166461 0.166461 0.740808 0.740808 0.740808 0.740808 octave:29> V>=W ans = 1 0 0 0 1 1 1 1 1 0 1 1 1 0 0 1 octave:30> idx=V>=W idx = 1 0 0 0 1 1 1 1 1 0 1 1 1 0 0 1 octave:31> ~idx ans = 0 1 1 1 0 0 0 0 0 1 0 0 0 1 1 0 octave:32> 1-idx ans = 0 1 1 1 0 0 0 0 0 1 0 0 0 1 1 0 octave:33> % NOOOO octave:33> octave:33> v=rand(1,4) v = 0.76157 0.89258 0.14991 0.89530 octave:34> P=metropolis(v) P = 0.268247 0.436243 0.042498 0.253011 0.372211 0.428711 0.011061 0.188016 0.215896 0.065856 0.199959 0.518288 0.215217 0.187444 0.086783 0.510556 octave:35> invariantPotenze(P) error: 'tol' undefined near line 17 column 15 error: called from: error: /home/accounts/personale/clrmrc90/aa1819/sistemi_stocastici/invariantPotenze.m at line 17, column 13 octave:35> invariantPotenze(P,1e-6) error: 'maxit' undefined near line 18 column 39 error: called from: error: /home/accounts/personale/clrmrc90/aa1819/sistemi_stocastici/invariantPotenze.m at line 18, column 1 octave:35> invariantPotenze(P,1e-6,100) ans = 0.282128 0.330663 0.055536 0.331673 octave:36> v v = 0.76157 0.89258 0.14991 0.89530 octave:37> v/sum(v) ans = 0.282128 0.330663 0.055536 0.331673 octave:38> A=randn(4) A = -1.845810 -1.259281 0.156176 0.894340 0.774733 -0.506133 -0.030375 0.042623 1.176704 0.302817 -0.963407 0.446974 0.599976 0.121794 -0.242482 1.048581 octave:39> A7=A^7; octave:40> A7_dummy=(((((A*A)*A)*A)*A)*A)*A; octave:41> A7-A7_dummy ans = 5.3291e-15 -7.1054e-15 1.7764e-15 1.7764e-15 -8.8818e-15 0.0000e+00 0.0000e+00 -8.8818e-16 -2.8422e-14 -2.1316e-14 3.5527e-15 4.4409e-15 -3.5527e-15 -1.7764e-15 0.0000e+00 2.2204e-16 octave:42> A2=A*A; octave:43> A4=A2*A2; octave:44> A6=A4*A2; octave:45> A7-A6*A ans = 1.7764e-15 -3.5527e-15 8.8818e-16 5.3291e-15 -8.8818e-15 0.0000e+00 1.7764e-15 1.7764e-15 -1.7764e-14 0.0000e+00 1.7764e-15 2.6645e-15 -1.7764e-15 -1.7764e-15 0.0000e+00 -2.2204e-16 octave:46> H= [-1,1] H = -1 1 octave:47> espilon=1e-3 espilon = 0.0010000 octave:48> exp(-H/epsilon) error: 'epsilon' undefined near line 1 column 8 error: evaluating argument list element number 1 octave:48> epsilon=1e-3 epsilon = 0.0010000 octave:49> exp(-H/epsilon) ans = Inf 0 octave:50> epsilon=1e-2 epsilon = 0.010000 octave:51> exp(-H/epsilon) ans = 2.6881e+43 3.7201e-44 octave:52> v=exp(-H/epsilon) v = 2.6881e+43 3.7201e-44 octave:53> v=v/sum(v) v = 1.0000e+00 1.3839e-87 octave:54> format long e octave:55> v=v/sum(v) v = 1.00000000000000e+00 1.38389652673674e-87 octave:56> epsilon=1e-3 epsilon = 1.00000000000000e-03 octave:57> v=exp(-H/epsilon) v = Inf 0.00000000000000e+00 octave:58> v=v/sum(v) v = NaN 0.00000000000000e+00 octave:59> v = 1./sum(exp(-(repmat(H’,1,2)-repmat(H,2,1))/epsilon)) error: invalid character 'â' (ASCII 226) near line 1, column 27 parse error: syntax error >>> v = 1./sum(exp(-(repmat(H’,1,2)-repmat(H,2,1))/epsilon)) ^ octave:59> v = 1./sum(exp(-(repmat(H',1,2)-repmat(H,2,1))/epsilon)) v = 1.00000000000000e+00 0.00000000000000e+00 octave:60> H H = -1.00000000000000e+00 1.00000000000000e+00 octave:61> epsilon=1e-1 epsilon = 1.00000000000000e-01 octave:62> v = 1./sum(exp(-(repmat(H',1,2)-repmat(H,2,1))/epsilon)) v = 9.99999997938846e-01 2.06115361819020e-09 octave:63> epsilon=1e-2 epsilon = 1.00000000000000e-02 octave:64> v = 1./sum(exp(-(repmat(H',1,2)-repmat(H,2,1))/epsilon)) v = 1.00000000000000e+00 1.38389652673674e-87 octave:65> epsmet ans = Columns 1 through 6: 2.6616e-01 8.8394e-02 6.2747e-02 8.8394e-02 1.0000e-01 8.8394e-02 1.0000e-01 1.9607e-01 7.0985e-02 1.0000e-01 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 1.0000e-01 7.0985e-02 1.9607e-01 1.0000e-01 1.0000e-01 5.7792e-02 5.1085e-02 3.6263e-02 5.1085e-02 5.3369e-01 5.1085e-02 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 1.9607e-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 1.0000e-01 2.1540e-02 1.9040e-02 1.3515e-02 1.9040e-02 3.7271e-02 1.9040e-02 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 1.0000e-01 Columns 7 through 10: 2.9121e-02 8.8394e-02 1.0000e-01 8.8394e-02 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 4.6411e-02 1.0000e-01 1.0000e-01 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.6830e-02 5.1085e-02 1.0000e-01 5.1085e-02 3.2945e-02 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 6.2726e-03 1.9040e-02 8.2620e-01 1.9040e-02 3.2945e-02 1.0000e-01 1.0000e-01 1.9607e-01 guarda P(5,5) e P(9,9) octave:65> epsmet ans = Columns 1 through 6: 2.6616e-01 8.8394e-02 6.2747e-02 8.8394e-02 1.0000e-01 8.8394e-02 1.0000e-01 1.9607e-01 7.0985e-02 1.0000e-01 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 1.0000e-01 7.0985e-02 1.9607e-01 1.0000e-01 1.0000e-01 5.7792e-02 5.1085e-02 3.6263e-02 5.1085e-02 5.3369e-01 5.1085e-02 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 1.9607e-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 1.0000e-01 2.1540e-02 1.9040e-02 1.3515e-02 1.9040e-02 3.7271e-02 1.9040e-02 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 1.0000e-01 Columns 7 through 10: 2.9121e-02 8.8394e-02 1.0000e-01 8.8394e-02 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 4.6411e-02 1.0000e-01 1.0000e-01 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.6830e-02 5.1085e-02 1.0000e-01 5.1085e-02 3.2945e-02 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 6.2726e-03 1.9040e-02 8.2620e-01 1.9040e-02 3.2945e-02 1.0000e-01 1.0000e-01 1.9607e-01 guarda P(5,5) e P(9,9) octave:65> epsmet ans = 2.6616e-01 8.8394e-02 6.2747e-02 8.8394e-02 1.0000e-01 8.8394e-02 2.9121e-02 8.8394e-02 1.0000e-01 8.8394e-02 1.0000e-01 1.9607e-01 7.0985e-02 1.0000e-01 1.0000e-01 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.5359e-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 1.0000e-01 7.0985e-02 1.9607e-01 1.0000e-01 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 5.7792e-02 5.1085e-02 3.6263e-02 5.1085e-02 5.3369e-01 5.1085e-02 1.6830e-02 5.1085e-02 1.0000e-01 5.1085e-02 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 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 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 1.0000e-01 3.2945e-02 1.9607e-01 1.0000e-01 1.0000e-01 2.1540e-02 1.9040e-02 1.3515e-02 1.9040e-02 3.7271e-02 1.9040e-02 6.2726e-03 1.9040e-02 8.2620e-01 1.9040e-02 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.9607e-01 guarda P(5,5) e P(9,9) octave:65> epsmet epsilon = 2.0000e+01 ans = 2.6616e-01 8.8394e-02 6.2747e-02 8.8394e-02 1.0000e-01 8.8394e-02 2.9121e-02 8.8394e-02 1.0000e-01 8.8394e-02 1.0000e-01 1.9607e-01 7.0985e-02 1.0000e-01 1.0000e-01 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.5359e-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 1.0000e-01 7.0985e-02 1.9607e-01 1.0000e-01 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.0000e-01 5.7792e-02 5.1085e-02 3.6263e-02 5.1085e-02 5.3369e-01 5.1085e-02 1.6830e-02 5.1085e-02 1.0000e-01 5.1085e-02 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 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 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 1.0000e-01 3.2945e-02 1.9607e-01 1.0000e-01 1.0000e-01 2.1540e-02 1.9040e-02 1.3515e-02 1.9040e-02 3.7271e-02 1.9040e-02 6.2726e-03 1.9040e-02 8.2620e-01 1.9040e-02 1.0000e-01 1.0000e-01 7.0985e-02 1.0000e-01 1.0000e-01 1.0000e-01 3.2945e-02 1.0000e-01 1.0000e-01 1.9607e-01 guarda P(5,5) e P(9,9) epsilon = 1.0000e+01 ans = 3.6148e-01 7.8134e-02 3.9371e-02 7.8134e-02 1.0000e-01 7.8134e-02 8.4805e-03 7.8134e-02 1.0000e-01 7.8134e-02 1.0000e-01 2.3876e-01 5.0389e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0854e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.7846e-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.0000e-01 5.0389e-02 2.3876e-01 1.0000e-01 1.0000e-01 1.0854e-02 1.0000e-01 1.0000e-01 1.0000e-01 3.3400e-02 2.6097e-02 1.3150e-02 2.6097e-02 7.2013e-01 2.6097e-02 2.8325e-03 2.6097e-02 1.0000e-01 2.6097e-02 1.0000e-01 1.0000e-01 5.0389e-02 1.0000e-01 1.0000e-01 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.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 1.0000e-01 1.0854e-02 2.3876e-01 1.0000e-01 1.0000e-01 4.6396e-03 3.6251e-03 1.8267e-03 3.6251e-03 1.3891e-02 3.6251e-03 3.9346e-04 3.6251e-03 9.6112e-01 3.6251e-03 1.0000e-01 1.0000e-01 5.0389e-02 1.0000e-01 1.0000e-01 1.0000e-01 1.0854e-02 1.0000e-01 1.0000e-01 2.3876e-01 guarda P(5,5) e P(9,9) epsilon = 1.0000e+00 ans = 7.5759e-01 8.4805e-03 8.9498e-06 8.4805e-03 1.0000e-01 8.4805e-03 1.9240e-12 8.4805e-03 1.0000e-01 8.4805e-03 1.0000e-01 2.9989e-01 1.0553e-04 1.0000e-01 1.0000e-01 1.0000e-01 2.2688e-11 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.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 1.0000e-01 1.0553e-04 2.9989e-01 1.0000e-01 1.0000e-01 2.2688e-11 1.0000e-01 1.0000e-01 1.0000e-01 1.7275e-06 1.4650e-07 1.5461e-10 1.4650e-07 9.0000e-01 1.4650e-07 3.3238e-17 1.4650e-07 1.0000e-01 1.4650e-07 1.0000e-01 1.0000e-01 1.0553e-04 1.0000e-01 1.0000e-01 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 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 1.0000e-01 2.2688e-11 2.9989e-01 1.0000e-01 1.0000e-01 4.6217e-15 3.9194e-16 4.1363e-19 3.9194e-16 2.6753e-10 3.9194e-16 8.8923e-26 3.9194e-16 1.0000e+00 3.9194e-16 1.0000e-01 1.0000e-01 1.0553e-04 1.0000e-01 1.0000e-01 1.0000e-01 2.2688e-11 1.0000e-01 1.0000e-01 2.9989e-01 guarda P(5,5) e P(9,9) epsilon = 1.0000e-01 ans = 8.0000e-01 1.9240e-12 3.2971e-42 1.9240e-12 1.0000e-01 1.9240e-12 6.9524e-109 1.9240e-12 1.0000e-01 1.9240e-12 1.0000e-01 3.0000e-01 1.7137e-31 1.0000e-01 1.0000e-01 1.0000e-01 3.6134e-98 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 2.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 1.0000e-01 1.7137e-31 3.0000e-01 1.0000e-01 1.0000e-01 3.6134e-98 1.0000e-01 1.0000e-01 1.0000e-01 2.3674e-49 4.5551e-60 7.8058e-90 4.5551e-60 9.0000e-01 4.5551e-60 1.6459e-156 4.5551e-60 1.0000e-01 4.5551e-60 1.0000e-01 1.0000e-01 1.7137e-31 1.0000e-01 1.0000e-01 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 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 1.0000e-01 3.6134e-98 3.0000e-01 1.0000e-01 1.0000e-01 4.4462e-135 8.5546e-146 1.4660e-175 8.5546e-146 1.8780e-87 8.5546e-146 3.0911e-242 8.5546e-146 1.0000e+00 8.5546e-146 1.0000e-01 1.0000e-01 1.7137e-31 1.0000e-01 1.0000e-01 1.0000e-01 3.6134e-98 1.0000e-01 1.0000e-01 3.0000e-01 guarda P(5,5) e P(9,9) epsilon = 1.0000e-03 ans = 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 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 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 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 guarda P(5,5) e P(9,9) octave:66> diary off