function H(a,b,c,d,e,f)R=M(3,4,f)*M(2,4,e)*M(4,1,d)*M(1,2,c)*M(1,3,b)*M(2,3,a)*(dec2bin(0:15)'-48.5);Z=R(3,:)+2;R=2*R./Z;Q=[1,2,10,12,11,9,10,14,16,12,4,8,16,15,11,3,7,15,13,9,1,5,13,14,6,8,7,5,6,2,4,3,1];plot(R(1,Q),R(2,Q));function A=M(a,b,v)C=@cosd;S=@sind;A=eye(4);A(sub2ind([4,4],[a,a,b,b],[a,b,a,b]))=[C(v),-S(v),S(v),C(v)]; function H(a,b,c,d,e,f) R=M(3,4,f)*M(2,4,e)*M(4,1,d)*M(1,2,c)*M(1,3,b)*M(2,3,a)*(dec2bin(0:15)'-48.5); Z=R(3,:)+2; R=2*R./Z; Q=[1,2,10,12,11,9,10,14,16,12,4,8,16,15,11,3,7,15,13,9,1,5,13,14,6,8,7,5,6,2,4,3,1]; plot(R(1,Q),R(2,Q)); function A=M(a,b,v) C=@cosd;S=@sind; A=eye(4);A(sub2ind([4,4],[a,a,b,b],[a,b,a,b]))=[C(v),-S(v),S(v),C(v)];