n=8; M=zeros(n);
for zei = 1:n
for spa = zei:n
if (zei^2 + spa^2)<1.5*n^2
M(zei,spa) = 2;
end
end
end
T1 = [1 0 0; 0 1 -5; 0 0 1]; Mx = [1 0 0; 0 -1 0; 0 0 1]; T2 = [1 0 0; 0 1 5; 0 0 1]; L = [2 2 3; 2 0 0; 1 1 1] T3 = [1 0 -2; 0 1 -10; 0 0 1]; R = [0 -1 0; 1 0 0; 0 0 1]; T4 = [1 0 2; 0 1 10; 0 0 1] Ttot = T4*R*T3 * T2*Mx*T1 % = [ 0 1 2; 1 0 8 ; 0 0 1]; % Ltr = [4 2 2; 10 10 11; 1 1 1]