diag([1,2,3,4]) ans = 1 0 0 0 0 2 0 0 0 0 3 0 0 0 0 4 spdiags([1;2;3;4],0,4,4) ans = (1,1) 1 (2,2) 2 (3,3) 3 (4,4) 4 spdiags([[1;2;3;4],[5;6;7;8]],[0,1],4,4)[5;6;7;8] ??? spdiags([[1;2;3;4],[5;6;7;8]],[0,1],4,4)[5;6;7;8] | Error: Unbalanced or unexpected parenthesis or bracket. spdiags([[1;2;3;4],[5;6;7;8]],[0,1],4,4) ans = (1,1) 1 (1,2) 6 (2,2) 2 (2,3) 7 (3,3) 3 (3,4) 8 (4,4) 4 full(ans) ans = 1 6 0 0 0 2 7 0 0 0 3 8 0 0 0 4 spdiags([[1;2;3;4],[5;6;7;8]],[0,-1],4,4) ans = (1,1) 1 (2,1) 5 (2,2) 2 (3,2) 6 (3,3) 3 (4,3) 7 (4,4) 4 full(ans) ans = 1 0 0 0 5 2 0 0 0 6 3 0 0 0 7 4 A=rand(4) A = 0.8147 0.6324 0.9575 0.9572 0.9058 0.0975 0.9649 0.4854 0.1270 0.2785 0.1576 0.8003 0.9134 0.5469 0.9706 0.1419 b=rand(4,1) b = 0.4218 0.9157 0.7922 0.9595 x=A\b x = 17.2819 0.8395 -15.9067 1.0883 A*x ans = 0.4218 0.9157 0.7922 0.9595 b b = 0.4218 0.9157 0.7922 0.9595 y=(A'\b)' y = 4.5439 -1.4105 -3.0671 -1.7661 b'/A ans = 4.5439 -1.4105 -3.0671 -1.7661 A/b ??? Error using ==> mrdivide Matrix dimensions must agree. A/b' ans = 1.0136 0.6634 0.4684 0.6984 [D,V]=eig(A') D = -0.5308 -0.6543 -0.3872 -0.0825 -0.3371 0.1107 -0.1442 -0.6424 -0.5902 0.6246 -0.3503 -0.0637 -0.5062 0.4117 0.8406 0.7593 V = 2.4021 0 0 0 0 -0.0346 0 0 0 0 -0.7158 0 0 0 0 -0.4400 [V,D]=eig(A') V = -0.5308 -0.6543 -0.3872 -0.0825 -0.3371 0.1107 -0.1442 -0.6424 -0.5902 0.6246 -0.3503 -0.0637 -0.5062 0.4117 0.8406 0.7593 D = 2.4021 0 0 0 0 -0.0346 0 0 0 0 -0.7158 0 0 0 0 -0.4400 A'*V ans = -1.2751 0.0226 0.2772 0.0363 -0.8097 -0.0038 0.1032 0.2826 -1.4178 -0.0216 0.2507 0.0280 -1.2159 -0.0142 -0.6017 -0.3341 V*D ans = -1.2751 0.0226 0.2772 0.0363 -0.8097 -0.0038 0.1032 0.2826 -1.4178 -0.0216 0.2507 0.0280 -1.2159 -0.0142 -0.6017 -0.3341 pi ans = 3.1416 v=V(:,1) v = -0.5308 -0.3371 -0.5902 -0.5062 v=v' v = -0.5308 -0.3371 -0.5902 -0.5062 v*A ans = -1.2751 -0.8097 -1.4178 -1.2159 v*D(1,1) ans = -1.2751 -0.8097 -1.4178 -1.2159 norm(v) ans = 1 uiopen('/home/accounts/personale/clrmrc90/aa1314/sistemi_stocastici/randstoc.m', true); P=rand(4) P = 0.6557 0.6787 0.6555 0.2769 0.0357 0.7577 0.1712 0.0462 0.8491 0.7431 0.7060 0.0971 0.9340 0.3922 0.0318 0.8235 sum(P) ans = 2.4746 2.5718 1.5645 1.2437 sum(P,2) ans = 2.2669 1.0108 2.3954 2.1815 P=randstoc(4) P = 0.3423 0.2162 0.0921 0.3495 0.1632 0.1964 0.2521 0.3884 0.3899 0.3141 0.1828 0.1132 0.0160 0.3689 0.2998 0.3153 sum(P,2) ans = 1 1 1 1 invariantPotenze(P,1e-6,100) ans = 0.2055 0.2776 0.2183 0.2985 v=ans v = 0.2055 0.2776 0.2183 0.2985 v*P ans = 0.2055 0.2776 0.2183 0.2985 [V,D]=eig(P') V = Columns 1 through 2 -0.4062 -0.7013 -0.5486 0.2490 - 0.1216i -0.4314 0.3570 - 0.3183i -0.5899 0.0953 + 0.4399i Columns 3 through 4 -0.7013 0.1933 0.2490 + 0.1216i -0.8162 0.3570 + 0.3183i 0.0852 0.0953 - 0.4399i 0.5378 D = Columns 1 through 2 1.0000 0 0 0.0837 + 0.1952i 0 0 0 0 Columns 3 through 4 0 0 0 0 0.0837 - 0.1952i 0 0 -0.1306 w=V(:,1) w = -0.4062 -0.5486 -0.4314 -0.5899 w=w/sum(w) w = 0.2055 0.2776 0.2183 0.2985 w=w' w = 0.2055 0.2776 0.2183 0.2985 v v = 0.2055 0.2776 0.2183 0.2985 norm(v-w,1) ans = 6.9278e-08 [0,0,0,0,1]/[P-eye(4),ones(4,1)] ans = 0.2055 0.2776 0.2183 0.2985 sum([1/3,1/3,1/3]) ans = 1 sum([1/3,1/3,1/3])==1 ans = 1 1/3 ans = 0.3333 format long e 1/3 ans = 3.333333333333333e-01 sum([1/3,1/3,1/3])==1 ans = 1 diary off