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