clf
for i = 0:pi/12:20*pi,
p1=[cos(i),sin(i)];
p2=[cos(i+pi/2),sin(i+pi/2)];
p3=[cos(i+pi),sin(i+pi)];
p4=[cos(i+pi/2*3),sin(i+pi/2*3)];
rect=[p1;p2;p3;p4;p1];
%cla;
fill(rect(:,1),rect(:,2),'b');
axis([-5 5 -5 5]);
axis equal;
drawnow;
%pause(0.1);
end