   hold on; axis equal; l = light;
    axis ([0 100 10 90  20 80]);
    view([40 20]);
 
    % decorate
    strTitle = sprintf('Ellipsoid, ori=[%d %d %d]', orientations(i, :));
    title(strTitle);

    snapnow;
    
    % add surfacic mesh
    drawEllipsoid(elli, 'FaceColor', 'r');
    
    % decorate
    title(strTitle);
    legend('isosurface', 'shape mesh', 'location', 'northeast');
    snapnow;
end


%% Cuboid

clf;

% iterate on orientations
for i = 1:4
    cubo = [center 90 