e="6"></a></h2><pre class="codeinput"><span class="comment">% iterate on orientations</span>
<span class="keyword">for</span> i = 1:4

    <span class="comment">% cylinder representation</span>
    cube = [center 60 orientations(i, :)];

    <span class="comment">% generation of 3D image</span>
    img = discreteCube(lx, ly, lz, cube);

    <span class="comment">% display image isosurface</span>
    clf; set(gca, <span class="string">'fontsize'</span>, 14);
    patch(isosurface(img, .5), <span class="string