randstocsym(4) ans = 0.2772 0.2971 0.1269 0.2988 0.2971 0.3534 0.1384 0.2111 0.1269 0.1384 0.2941 0.4405 0.2988 0.2111 0.4405 0.0496 A=ans A = 0.2772 0.2971 0.1269 0.2988 0.2971 0.3534 0.1384 0.2111 0.1269 0.1384 0.2941 0.4405 0.2988 0.2111 0.4405 0.0496 sum(A,2) ans = 1.0000 1.0000 1.0000 1.0000 sum(A,1) ans = 1.0000 1.0000 1.0000 1.0000 A=ones(4)/4 A = 0.2500 0.2500 0.2500 0.2500 0.2500 0.2500 0.2500 0.2500 0.2500 0.2500 0.2500 0.2500 0.2500 0.2500 0.2500 0.2500 uiopen('/home/accounts/personale/clrmrc90/aa1213/sistemi_stocastici/metropolis.m', true); v=[1,2,3,4] v = 1 2 3 4 repmat(v,4,1) ans = 1 2 3 4 1 2 3 4 1 2 3 4 1 2 3 4 repmat(v,4,2) ans = 1 2 3 4 1 2 3 4 1 2 3 4 1 2 3 4 1 2 3 4 1 2 3 4 1 2 3 4 1 2 3 4 V=repmat(v,4,1) V = 1 2 3 4 1 2 3 4 1 2 3 4 1 2 3 4 W=repmat(v',1,4) W = 1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 V V = 1 2 3 4 1 2 3 4 1 2 3 4 1 2 3 4 W W = 1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 1==2 ans = 0 [1,2]==[1,3] ans = 1 0 V V = 1 2 3 4 1 2 3 4 1 2 3 4 1 2 3 4 W W = 1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 V>=W ans = 1 1 1 1 0 1 1 1 0 0 1 1 0 0 0 1 v=rand(1,4) v = 0.1576 0.9706 0.9572 0.4854 V=repmat(v,4,1) V = 0.1576 0.9706 0.9572 0.4854 0.1576 0.9706 0.9572 0.4854 0.1576 0.9706 0.9572 0.4854 0.1576 0.9706 0.9572 0.4854 W=repmat(v',1,4) W = 0.1576 0.1576 0.1576 0.1576 0.9706 0.9706 0.9706 0.9706 0.9572 0.9572 0.9572 0.9572 0.4854 0.4854 0.4854 0.4854 V>=W ans = 1 1 1 1 0 1 0 0 0 1 1 0 0 1 1 1 P=zeros(4) P = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 find(V>=W) ans = 1 5 6 7 8 9 11 12 13 16 P(V>=W)=2 P = 2 2 2 2 0 2 0 0 0 2 2 0 0 2 2 2 Q=randstocsym(4) Q = 0.3114 0.1358 0.2105 0.3422 0.1358 0.2856 0.3252 0.2533 0.2105 0.3252 0.1088 0.3555 0.3422 0.2533 0.3555 0.0490 P(V>=W)=Q(V>=W) P = 0.3114 0.1358 0.2105 0.3422 0 0.2856 0 0 0 0.3252 0.1088 0 0 0.2533 0.3555 0.0490 P(VFile: invariantPotenze.m Line: 1 Column: 3 Function definitions are not permitted at the prompt or in scripts. invariantPotenze(P,1e-6,100) ans = 0.0613 0.3776 0.3723 0.1888 v/sum(v) ans = 0.0613 0.3776 0.3723 0.1888 P^10 ans = 0.0613 0.3775 0.3723 0.1888 0.0613 0.3776 0.3723 0.1888 0.0613 0.3776 0.3723 0.1888 0.0613 0.3775 0.3723 0.1888 P P = 0.3114 0.1358 0.2105 0.3422 0.0221 0.5305 0.3208 0.1267 0.0347 0.3252 0.4598 0.1803 0.1111 0.2533 0.3555 0.2801 P^2 ans = 0.1453 0.2695 0.3276 0.2576 0.0438 0.4209 0.3673 0.1680 0.0539 0.3725 0.3871 0.1864 0.0836 0.3360 0.3677 0.2127 P^5 ans = 0.0653 0.3714 0.3709 0.1924 0.0603 0.3792 0.3726 0.1879 0.0611 0.3778 0.3725 0.1886 0.0625 0.3757 0.3719 0.1899 P^10 ans = 0.0613 0.3775 0.3723 0.1888 0.0613 0.3776 0.3723 0.1888 0.0613 0.3776 0.3723 0.1888 0.0613 0.3775 0.3723 0.1888 w=rand(1,4) w = 0.6555 0.1712 0.7060 0.0318 w=w/sum(w) w = 0.4190 0.1094 0.4513 0.0203 w=w*P w = 0.1508 0.2669 0.3380 0.2443 w=w*P w = 0.0917 0.3339 0.3596 0.2148 w=w*P w = 0.0723 0.3610 0.3681 0.1987 w=w*P w = 0.0653 0.3714 0.3709 0.1925 w=w*P w = 0.0628 0.3753 0.3718 0.1902 w=w*P w = 0.0618 0.3767 0.3721 0.1893 w=w*P w = 0.0615 0.3772 0.3723 0.1890 w=w*P w = 0.0614 0.3774 0.3723 0.1889 w=w*P w = 0.0613 0.3775 0.3723 0.1888 w=w*P w = 0.0613 0.3775 0.3723 0.1888 v v = 0.1576 0.9706 0.9572 0.4854 v/sum(v) ans = 0.0613 0.3776 0.3723 0.1888 H(1)=-1; H(2)=1; epsilon=1e-3; v=exp(-H/epsilon) v = Inf 0 v=v/sum(v) v = NaN 0 v=[1e309,2e309] v = Inf Inf v=[1e308,2e308] v = 1.0e+308 * 1.0000 Inf v=[3e307,4e307] v = 1.0e+307 * 3.0000 4.0000 norm(v) ans = 5.0000e+307 sqrt(v(1)^2+v(2)^2) ans = Inf v = sum(exp(-(repmat(H’,1,2)-repmat(H,2,1))/epsilon)) ??? v = sum(exp(-(repmat(H’,1,2)-repmat(H,2,1))/epsilon)) | Error: The input character is not valid in MATLAB statements or expressions. v = 1/sum(exp(-(repmat(H',1,2)-repmat(H,2,1))/epsilon)) ??? Error using ==> mrdivide Matrix dimensions must agree. v = 1./sum(exp(-(repmat(H',1,2)-repmat(H,2,1))/epsilon)) v = 1 0 epsilon=1e-1; v = 1./sum(exp(-(repmat(H',1,2)-repmat(H,2,1))/epsilon)) v = 1.0000 0.0000 epsilon=0.5; v = 1./sum(exp(-(repmat(H',1,2)-repmat(H,2,1))/epsilon)) v = 0.9820 0.0180 H H = -1 1 epsilon=0.1; v = 1./sum(exp(-(repmat(H',1,2)-repmat(H,2,1))/epsilon)) v = 1.0000 0.0000 format long e v = 1./sum(exp(-(repmat(H',1,2)-repmat(H,2,1))/epsilon)) v = Column 1 9.999999979388463e-01 Column 2 2.061153618190204e-09 epsilon=0.01; v = 1./sum(exp(-(repmat(H',1,2)-repmat(H,2,1))/epsilon)) v = Column 1 1.000000000000000e+00 Column 2 1.383896526736738e-87 % Modify expression to add input arguments. % Example: % a = [1 2 3; 4 5 6]; % foo(a); epsmet epsilon = 10 ans = Columns 1 through 3 3.8317e-01 8.8844e-03 9.4179e-03 1.1371e-02 3.3859e-01 7.8145e-02 2.3921e-02 1.5508e-01 2.1976e-01 2.0279e-01 3.6444e-02 6.7750e-02 5.7152e-02 2.4926e-02 4.3887e-03 7.8092e-02 8.6899e-02 1.2307e-02 2.3401e-01 1.2605e-01 1.0228e-01 8.4830e-03 1.3834e-01 9.9256e-02 5.0131e-03 5.3355e-03 1.2761e-03 9.3966e-02 5.3831e-02 6.0528e-02 Columns 4 through 6 1.5845e-01 1.7112e-01 6.1017e-02 3.6444e-02 9.5515e-02 8.6899e-02 1.3445e-01 3.3375e-02 2.4423e-02 1.6533e-01 1.0980e-01 3.7282e-02 2.8653e-02 6.1113e-01 9.7599e-03 3.7282e-02 3.7399e-02 3.1168e-01 7.3945e-02 4.0278e-02 7.0732e-02 1.0217e-01 6.9472e-02 2.6991e-01 4.7201e-03 3.1507e-02 3.6851e-03 1.4020e-01 6.8602e-02 5.7101e-02 Columns 7 through 9 1.9845e-02 6.6282e-03 1.0805e-01 1.3681e-02 1.3834e-01 1.4718e-01 2.2032e-02 1.9698e-01 6.9861e-02 8.0258e-03 1.0217e-01 1.3021e-01 1.1408e-03 1.8130e-02 2.2681e-01 7.6771e-03 2.6991e-01 1.0165e-01 5.2337e-02 1.2841e-01 9.8654e-02 1.3937e-02 2.4864e-01 2.5675e-02 3.8817e-04 9.3073e-04 9.4635e-01 7.9559e-03 2.4115e-02 2.1837e-02 Column 10 7.3420e-02 5.3831e-02 1.2012e-01 1.4020e-01 1.7903e-02 5.7101e-02 7.3301e-02 2.4115e-02 7.9161e-04 4.7186e-01 guarda P(5,5) e P(9,9) epsilon = 1.0000e-01 ans = Columns 1 through 3 7.2083e-01 2.1878e-13 7.8870e-43 1.1371e-02 4.3042e-01 2.6576e-31 2.3921e-02 1.5508e-01 2.4179e-01 2.0279e-01 3.6444e-02 2.3041e-31 4.0511e-49 4.3508e-60 2.6052e-90 7.8092e-02 8.6899e-02 4.1853e-32 2.3401e-01 1.2605e-01 1.0228e-01 8.4830e-03 1.3834e-01 3.3755e-31 4.8041e-135 1.2591e-145 1.0241e-175 9.3966e-02 5.3831e-02 2.0585e-31 Columns 4 through 6 3.9018e-12 1.7112e-01 1.5025e-12 3.6444e-02 9.5515e-02 8.6899e-02 1.3445e-01 3.3375e-02 2.4423e-02 2.4111e-01 1.0980e-01 3.7282e-02 5.0013e-60 7.7319e-01 1.7036e-60 3.7282e-02 3.7399e-02 3.3167e-01 7.3945e-02 4.0278e-02 7.0732e-02 1.0217e-01 6.9472e-02 2.6991e-01 1.1139e-145 4.2597e-87 8.6960e-146 1.4020e-01 6.8602e-02 5.7101e-02 Columns 7 through 9 1.6269e-108 1.6322e-13 1.0805e-01 4.5545e-98 1.3834e-01 1.4718e-01 2.1568e-68 1.9698e-01 6.9861e-02 2.6720e-98 1.0217e-01 1.3021e-01 6.6294e-157 3.1645e-60 2.2681e-01 2.5559e-98 2.6991e-01 1.0165e-01 5.2337e-02 1.2841e-01 9.8654e-02 4.6400e-98 3.6184e-01 2.5675e-02 3.0495e-242 2.1964e-146 1.0000e+00 2.6487e-98 2.4115e-02 2.1837e-02 Column 10 1.8079e-12 5.3831e-02 1.2012e-01 1.4020e-01 3.1249e-60 5.7101e-02 7.3301e-02 2.4115e-02 1.8681e-146 5.4035e-01 guarda P(5,5) e P(9,9) epsilon = 1.0000e-02 ans = Columns 1 through 3 7.2083e-01 7.9053e-110 0 1.1371e-02 4.3042e-01 3.3870e-299 2.3921e-02 1.5508e-01 2.4179e-01 2.0279e-01 3.6444e-02 2.9365e-299 0 0 0 7.8092e-02 8.6899e-02 5.3340e-300 2.3401e-01 1.2605e-01 1.0228e-01 8.4830e-03 1.3834e-01 4.3020e-299 0 0 0 9.3966e-02 5.3831e-02 2.6235e-299 Columns 4 through 6 1.4099e-108 1.7112e-01 5.4292e-109 3.6444e-02 9.5515e-02 8.6899e-02 1.3445e-01 3.3375e-02 2.4423e-02 2.4111e-01 1.0980e-01 3.7282e-02 0 7.7319e-01 0 3.7282e-02 3.7399e-02 3.3167e-01 7.3945e-02 4.0278e-02 7.0732e-02 1.0217e-01 6.9472e-02 2.6991e-01 0 0 0 1.4020e-01 6.8602e-02 5.7101e-02 Columns 7 through 9 0 5.8977e-110 1.0805e-01 0 1.3834e-01 1.4718e-01 0 1.9698e-01 6.9861e-02 0 1.0217e-01 1.3021e-01 0 0 2.2681e-01 0 2.6991e-01 1.0165e-01 5.2337e-02 1.2841e-01 9.8654e-02 0 3.6184e-01 2.5675e-02 0 0 1.0000e+00 0 2.4115e-02 2.1837e-02 Column 10 6.5329e-109 5.3831e-02 1.2012e-01 1.4020e-01 0 5.7101e-02 7.3301e-02 2.4115e-02 0 5.4035e-01 guarda P(5,5) e P(9,9) epsmet ??? Error: File: epsmet.m Line: 13 Column: 27 Unexpected MATLAB expression. epsmet epsilon = 10 ??? Error using ==> minus Matrix dimensions must agree. Error in ==> epsmet at 13 v = 1./sum(exp(-(repmat(H',1,2)-repmat(H,2,1))/epsilon)); epsmet epsilon = 10 ans = Columns 1 through 3 3.9605e-01 7.3972e-02 4.1325e-02 9.4673e-02 2.2986e-01 3.7606e-02 1.0496e-01 7.4631e-02 1.6687e-01 7.7947e-02 9.7832e-02 5.5996e-02 2.5632e-02 2.2274e-02 1.4877e-02 1.0293e-01 1.3045e-01 5.8702e-02 1.1657e-01 9.5984e-02 9.9505e-02 1.2507e-01 1.1149e-01 3.9590e-02 3.7541e-03 3.1031e-03 1.5866e-03 1.0502e-01 1.1673e-01 6.3454e-02 Columns 4 through 6 6.0904e-02 7.6745e-02 8.0423e-02 9.7832e-02 8.5352e-02 1.3045e-01 1.1113e-01 1.1313e-01 1.1650e-01 2.5983e-01 1.2549e-01 1.0749e-01 3.2749e-02 7.2057e-01 1.9309e-02 1.0749e-01 7.3989e-02 2.1983e-01 1.3615e-01 1.0833e-01 8.3033e-02 7.8677e-02 1.1022e-01 1.0226e-01 3.6795e-03 1.5630e-02 3.2314e-03 8.0460e-02 7.7547e-02 1.0620e-01 Columns 7 through 9 9.8855e-03 9.7723e-02 8.0915e-02 1.0418e-02 1.1149e-01 8.5599e-02 2.1433e-02 7.8568e-02 8.6855e-02 1.4778e-02 7.8677e-02 1.0150e-01 3.0685e-03 2.8765e-02 1.1252e-01 9.0122e-03 1.0226e-01 8.9140e-02 9.0108e-02 9.7466e-02 8.6849e-02 1.0579e-02 2.3222e-01 9.5504e-02 3.4171e-04 3.4621e-03 9.5902e-01 9.3340e-03 9.4385e-02 1.7085e-01 Column 10 8.2055e-02 1.1673e-01 1.2593e-01 8.0460e-02 2.0237e-02 1.0620e-01 8.5998e-02 9.4385e-02 6.1934e-03 1.7603e-01 guarda P(5,5) e P(9,9) ??? Operation terminated by user during ==> epsmet at 16 epsmet epsilon = 10 ans = Columns 1 through 3 3.7827e-01 1.5448e-01 6.0939e-02 1.9771e-01 3.2745e-01 7.1713e-02 1.5478e-01 1.4232e-01 9.2779e-02 2.9263e-02 7.3941e-02 2.2906e-02 3.0867e-02 1.9233e-02 9.9900e-03 4.7504e-02 5.0872e-02 9.1414e-03 7.0843e-02 1.3558e-01 2.4822e-02 9.2432e-02 5.5621e-02 8.9283e-02 7.4575e-03 6.0724e-04 3.2850e-03 1.9127e-02 1.1753e-01 5.4515e-02 Columns 4 through 6 2.2865e-02 9.2418e-02 3.7117e-02 7.3941e-02 7.3698e-02 5.0872e-02 4.5458e-02 7.5970e-02 1.8142e-02 4.1996e-01 1.1783e-01 1.7724e-01 3.0749e-02 7.2591e-01 3.6784e-02 1.7724e-01 1.4095e-01 2.4661e-01 4.1213e-01 1.2478e-01 4.9168e-02 7.7305e-03 8.6862e-02 1.7871e-01 7.8275e-04 1.4637e-02 1.7311e-03 8.4811e-02 5.7081e-02 9.5884e-02 Columns 7 through 9 6.0078e-03 7.2221e-02 1.6074e-01 1.4715e-02 5.5621e-02 1.6751e-02 5.3466e-03 1.7719e-01 1.7983e-01 4.4732e-02 7.7305e-03 2.1592e-02 3.5343e-03 2.2668e-02 1.0537e-01 5.3366e-03 1.7871e-01 4.7752e-02 4.7297e-02 5.8986e-02 6.1332e-03 6.4021e-03 2.9691e-01 9.8162e-02 2.4132e-05 3.5585e-03 9.6255e-01 7.6261e-03 8.7891e-02 1.4793e-01 Column 10 1.4945e-02 1.1753e-01 1.0819e-01 8.4811e-02 1.4896e-02 9.5884e-02 7.0262e-02 8.7891e-02 5.3626e-03 3.2760e-01 guarda P(5,5) e P(9,9) epsilon = 1.0000e-01 ??? Undefined function or variable 'Q'. Error in ==> epsmet at 19 epsmet ans = Columns 1 through 3 2.3899e-01 1.0858e-01 7.4927e-02 1.2284e-01 2.2897e-01 9.6676e-02 1.1941e-01 1.3619e-01 1.0046e-01 9.6799e-02 1.2844e-01 3.1003e-02 3.2873e-02 2.3219e-02 2.1122e-02 1.2193e-01 6.5910e-02 7.9130e-02 8.0049e-02 1.0303e-01 7.9596e-02 5.2692e-02 5.0429e-02 1.6770e-01 3.0386e-02 3.5171e-02 1.4876e-02 1.3160e-01 4.2613e-02 3.3566e-02 Columns 4 through 6 8.5564e-02 5.6880e-02 1.0778e-01 1.2844e-01 4.5452e-02 6.5910e-02 4.3676e-02 5.8248e-02 1.1147e-01 2.0638e-01 1.9378e-01 8.6803e-02 9.8994e-02 5.8741e-01 8.9883e-02 8.6803e-02 1.7595e-01 1.9757e-01 2.1978e-02 2.0758e-01 1.1619e-01 5.1042e-02 6.4720e-02 1.3911e-01 1.5392e-02 1.2773e-02 4.4599e-03 1.1766e-01 8.6589e-02 7.1898e-02 Columns 7 through 9 2.3311e-02 4.6576e-02 1.4107e-01 3.3944e-02 5.0429e-02 1.8473e-01 3.6941e-02 2.3625e-01 1.1007e-01 7.2407e-03 5.1042e-02 8.0840e-02 3.4935e-02 3.3062e-02 3.4270e-02 3.8277e-02 1.3911e-01 2.3424e-02 1.2600e-01 4.5963e-03 1.4611e-01 1.5142e-03 2.2190e-01 1.7792e-01 9.1647e-03 3.3876e-02 8.3080e-01 3.7847e-02 7.2975e-02 6.8795e-02 Column 10 1.1632e-01 4.2613e-02 4.7285e-02 1.1766e-01 4.4234e-02 7.1898e-02 1.1488e-01 7.2975e-02 1.3098e-02 3.3646e-01 guarda P(5,5) e P(9,9) ans = Columns 1 through 3 3.5200e-01 3.4675e-02 3.6884e-02 4.4378e-02 4.5171e-01 9.4671e-02 9.3683e-02 1.8788e-01 1.7901e-01 1.1958e-01 2.3294e-02 3.2708e-02 4.3465e-02 1.6126e-02 1.3456e-02 7.5789e-02 7.9136e-02 4.7326e-02 7.0408e-02 1.6040e-01 1.1871e-01 1.8932e-01 3.2223e-02 5.8968e-02 3.3551e-04 6.1696e-03 8.0178e-04 1.6962e-01 2.5193e-02 4.6248e-02 Columns 4 through 6 9.3435e-02 1.3014e-01 5.9217e-02 2.3294e-02 6.1795e-02 7.9136e-02 6.4911e-02 1.0233e-01 9.3921e-02 1.7926e-01 6.9787e-03 1.7002e-01 1.8212e-03 7.8821e-01 4.7749e-02 1.7002e-01 1.8297e-01 1.8891e-01 1.1426e-01 1.7149e-01 4.4742e-03 1.0966e-01 1.2278e-01 7.6443e-02 7.4943e-03 5.4097e-03 9.4988e-04 1.3936e-01 5.1080e-02 1.5271e-01 Columns 7 through 9 5.9710e-03 1.4792e-01 7.2315e-03 1.7409e-02 3.2223e-02 1.7019e-01 2.5569e-02 1.1702e-01 4.3893e-02 1.2401e-02 1.0966e-01 2.0673e-01 4.8573e-03 3.2042e-02 3.8943e-02 4.8562e-04 7.6443e-02 2.6203e-02 1.4724e-01 1.0339e-01 9.7323e-02 1.1222e-02 3.6373e-01 1.3284e-02 3.8293e-04 4.8155e-04 9.7573e-01 1.3368e-03 2.2367e-02 6.1975e-02 Column 10 1.3253e-01 2.5193e-02 9.1782e-02 1.3936e-01 1.3330e-02 1.5271e-01 1.2317e-02 2.2367e-02 2.2467e-03 3.3012e-01 guarda P(5,5) e P(9,9) ans = Columns 1 through 3 7.6458e-01 8.7623e-03 7.6224e-06 1.0332e-01 3.2248e-01 8.8983e-05 8.5168e-02 8.4316e-02 1.9962e-01 1.0133e-01 8.6639e-02 1.0570e-04 1.6733e-06 1.2968e-07 1.7244e-10 1.1127e-01 9.8684e-02 1.1626e-04 1.0099e-01 1.1927e-01 9.4247e-02 1.0579e-01 1.1796e-01 1.0644e-04 4.3816e-15 3.2626e-16 3.3566e-19 9.4126e-02 9.9073e-02 1.3406e-04 Columns 4 through 6 8.5933e-03 9.6859e-02 9.4359e-03 8.6639e-02 8.8518e-02 9.8684e-02 1.0016e-01 1.1153e-01 1.1016e-01 2.8162e-01 8.6444e-02 1.0155e-01 1.2664e-07 9.0036e-01 1.2132e-07 1.0155e-01 8.2811e-02 3.0887e-01 9.1660e-02 1.1396e-01 9.2082e-02 1.0821e-01 9.2651e-02 1.0787e-01 3.9070e-16 2.6657e-10 3.4818e-16 1.3442e-01 1.0532e-01 9.9995e-02 Columns 7 through 9 1.9432e-12 8.9717e-03 9.4805e-02 2.7059e-11 1.1796e-01 8.3242e-02 2.0261e-08 1.0086e-01 8.1150e-02 2.0796e-11 1.0821e-01 9.9683e-02 3.7879e-17 1.3574e-07 9.9640e-02 2.0891e-11 1.0787e-01 8.8835e-02 9.7560e-02 9.4287e-02 9.8377e-02 2.1392e-11 2.7971e-01 8.3469e-02 8.7480e-26 3.2715e-16 1.0000e+00 2.2135e-11 1.0424e-01 9.8014e-02 Column 10 7.9823e-03 9.9073e-02 1.2703e-01 1.3442e-01 1.5430e-07 9.9995e-02 9.7565e-02 1.0424e-01 3.8416e-16 2.6467e-01 guarda P(5,5) e P(9,9) ans = Columns 1 through 3 8.0209e-01 2.3050e-12 1.1241e-42 1.1980e-01 2.2846e-01 1.9834e-31 3.4093e-02 1.1574e-01 2.3336e-01 6.6046e-02 1.3027e-01 1.9384e-31 1.9527e-49 1.4848e-60 3.7957e-90 1.7574e-01 1.5932e-01 2.5895e-31 2.8001e-02 6.1128e-02 1.8013e-01 1.5316e-01 1.0395e-01 2.2939e-31 5.1323e-135 1.3478e-145 9.2406e-176 6.7363e-02 6.8052e-02 1.8347e-31 Columns 4 through 6 1.2708e-12 8.2483e-02 3.3813e-12 1.3027e-01 3.2596e-02 1.5932e-01 1.1311e-01 4.8627e-02 1.5111e-01 2.9119e-01 1.7438e-01 1.6925e-01 7.9432e-60 8.2454e-01 5.9101e-60 1.6925e-01 1.2975e-01 2.2513e-01 1.5734e-01 6.7001e-02 3.2595e-02 5.0166e-02 7.4171e-02 4.8170e-02 9.7829e-146 3.2951e-87 4.1441e-146 4.3316e-03 3.9096e-02 4.4196e-02 Columns 7 through 9 1.9467e-109 2.9469e-12 1.1543e-01 2.2088e-98 1.0395e-01 1.5756e-01 3.7981e-68 1.3386e-01 6.3034e-02 5.6855e-98 5.0166e-02 1.1436e-01 1.1028e-156 3.3785e-60 1.7546e-01 1.1778e-98 4.8170e-02 4.8443e-02 1.6709e-02 3.4213e-02 1.5892e-01 1.2363e-98 3.0797e-01 8.7135e-02 4.9124e-242 7.4540e-146 1.0000e+00 9.5382e-98 1.7528e-01 4.5745e-02 Column 10 1.2961e-12 6.8052e-02 1.0706e-01 4.3316e-03 1.7808e-60 4.4196e-02 2.6397e-01 1.7528e-01 3.9133e-146 5.5593e-01 guarda P(5,5) e P(9,9) ans = Columns 1 through 3 9.8156e-02 8.0132e-02 1.2162e-01 8.0132e-02 9.6465e-02 9.3158e-02 1.2162e-01 9.3158e-02 8.8075e-02 1.3540e-01 1.1298e-01 1.1008e-01 1.0623e-01 9.6034e-02 1.1142e-01 7.8917e-02 1.1622e-01 7.5368e-02 9.3339e-02 1.0171e-01 8.6105e-02 8.4076e-02 1.1689e-01 8.4337e-02 0 0 0 1.0759e-01 1.0353e-01 1.2000e-01 Columns 4 through 6 1.3540e-01 1.0623e-01 7.8917e-02 1.1298e-01 9.6034e-02 1.1622e-01 1.1008e-01 1.1142e-01 7.5368e-02 8.4574e-02 1.0118e-01 9.7171e-02 1.0118e-01 7.1533e-02 1.0108e-01 9.7171e-02 1.0108e-01 1.3380e-01 7.1715e-02 9.8641e-02 1.0891e-01 9.4385e-02 1.0098e-01 7.4301e-02 0 0 0 1.0633e-01 9.2115e-02 1.2880e-01 Columns 7 through 9 9.3339e-02 8.4076e-02 9.4532e-02 1.0171e-01 1.1689e-01 8.2877e-02 8.6105e-02 8.4337e-02 1.0985e-01 7.1715e-02 9.4385e-02 8.6184e-02 9.8641e-02 1.0098e-01 1.2079e-01 1.0891e-01 7.4301e-02 8.5432e-02 1.1769e-01 8.6370e-02 1.0476e-01 8.6370e-02 9.5044e-02 1.6229e-01 0 0 1.0000e+00 1.3076e-01 1.0133e-01 7.3810e-02 Column 10 1.0759e-01 1.0353e-01 1.2000e-01 1.0633e-01 9.2115e-02 1.2880e-01 1.3076e-01 1.0133e-01 0 3.5730e-02 guarda P(5,5) e P(9,9) diary off