1;

k=1;
epsilon=[.5 1 2];
for i=1:3
  for j=1:3
    axis off square
    subplot (3,3, k)
#    supertoroid([epsilon(i) epsilon(j)], [2 1 1 3])    
    superhyper([epsilon(i) epsilon(j)], [1 1 1])    
#    superellipse([epsilon(i) epsilon(j)], [1 1 1])    
    k=k+1;
  endfor;
axis off square
endfor;

print -dpng "se.png"