terate on orientations
for i = 1:3
    torus = [center 30 10 angles(i, :)];
    
    % generation of 3D image
    img = discreteTorus(lx, ly, lz, torus);
    
    % display image isosurface
    clf; set(gca, 'fontsize', 14);
    patch(isosurface(img, .5), 'FaceColor', 'g', 'LineStyle', 'none');
        
    % setup display
    hold on; axis equal; l = light;
    axis ([0 100 10 90  20 80]);
    view([40 20]);
 
    % decorate
    strTitle = sprintf('Torus, ori=[%d %d]', angles(i, :));
    