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