pkg load image ax = axes('XLim',[-1.5 1.5],'YLim',[-1.5 1.5],'ZLim',[-1.5 1.5]); view(3) grid on [x,y,z] = cylinder([.2 0]); h(1) = surface(x,y,z,'FaceColor','red'); h(2) = surface(x,y,-z,'FaceColor','green'); h(3) = surface(z,x,y,'FaceColor','blue'); h(4) = surface(-z,x,y,'FaceColor','cyan'); h(5) = surface(y,z,x,'FaceColor','magenta'); h(6) = surface(y,-z,x,'FaceColor','yellow'); %get(h, "") t = hgtransform(ax); set(h,'Parent',t); view(3) Rz = eye(4); Sxy = Rz; pause(1) for r = 1:.1:2*pi % Z-axis rotation matrix Rz2 = makehgtform('zrotate', r); % Scaling matrix Sxy = makehgtform('scale',(r/4)); % Concatenate the transforms and % set the transform Matrix property set(t,'Matrix',Rz2); drawnow pause(0.2) end pause(1)