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)