X = [2,3,4 ; 3,3,4 ; 4,4,4 ; 5,4,4 ; 5,4,5 ; 4,4,5 ; 3,3,5 ; 2,3,5 ; 2,3,4];
Xprim = [1 0 0 0; 0 1 0 0; 0 0 1 0; 3 2 1 1];
hold on

theta = 30;
OX = [1 0 0 0; 0 cos(theta) sin(theta) 0; 0 -sin(theta) cos(theta) 0; 0 0 0 1];

for theta = 30:30:360
   rad = (pi * theta)/180;
   OY = [cos(rad) sin(rad) 0 0; -sin(rad) cos(rad) 0 0; 0 0 1 0; 0 0 0 1];
   res = xj * OY;
   plot3(res(:,1), res(:,2), res(:,3));
end