function cube(x1,y1,z1,x2,y2,z2,c)
if nargin < 7
  c = 'r'
end
patch([x1 x2 x2 x1], [y1 y1 y2 y2], [z2 z2 z2 z2],c);
hold on;
patch([x1 x2 x2 x1], [y2 y2 y1 y1], [z1 z1 z1 z1],c);
patch([x1 x2 x2 x1], [y1 y1 y1 y1], [z1 z1 z2 z2],c);
patch([x1 x2 x2 x1], [y2 y2 y2 y2], [z2 z2 z1 z1],c);
patch([x1 x1 x1 x1], [y1 y1 y2 y2], [z1 z2 z2 z1],c);
patch([x2 x2 x2 x2], [y2 y2 y1 y1], [z1 z2 z2 z1],c);
