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



%% Cube

% iterate on orientations
for i = 1:4
    
    % cylinder representation
    cube = [center 60 orientations(i, :)];
    
    % generation of 3D image
    img = discreteCube(lx, ly, lz, cube);
    
    % display image isosurface
    clf; set(gca, 'fontsize', 14);
    patch(isosurface(img, .