x=[-1,1,1,-1,-1,-1,1,1,-1,-1,1,1,1,1,-1,-1];
y=[-1,-1,1,1,-1,-1,-1,1,1,-1,-1,-1,1,1,1,1];
z=[-1,-1,-1,-1,-1,1,1,1,1,1,1,-1,-1,1,1,-1];
U=[x;y;z];
a=50;
Rx=[1,0,0;0,cos(2*pi/a), -sin(2*pi/a); 0,sin(2*pi/a), cos(2*pi/a)];
Ry=[cos(2*pi/a),0, -sin(2*pi/a); 0,1,0; sin(2*pi/a), 0, cos(2*pi/a)];
Rz=[cos(2*pi/a), -sin(2*pi/a),0; sin(2*pi/a), cos(2*pi/a),0;0,0,1];
for k=0:100
  clf
  axis([-2,2,-2,2,-2,2])
  p=rand;
  if  k < 33
    R=Rx;
  elseif k < 66
    R=Ry;
  else 
    R=Rz;
 end  
 U=R*[U(1,:);U(2,:);U(3,:)];
 line(U(1,:),U(2,:),U(3,:));
 F(k+1) = getframe;
end
movie(F)