n=8; M=zeros(n); for zei = 1:n for spa = zei:n if (zei^2 + spa^2)<1.2*n^2 M(zei,spa) = 4; 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 = [4 4 5; 2 0 0; 1 1 1] T3 = [1 0 -4; 0 1 -10; 0 0 1]; R = [0 -1 0; 1 0 0; 0 0 1]; T4 = [1 0 4; 0 1 10; 0 0 1] Ttot = T4*R*T3 * T2*Mx*T1 % = [ 0 1 4; 1 0 6 ; 0 0 1]; % Ltr = [6 4 4; 10 10 11; 1 1 1]