
lz = 1:100;

<span class="comment">% choose a center not aligned with the grid</span>
center = [50+sqrt(2)-1 50+sqrt(3)-1 50+sqrt(5)-2];

angles = [<span class="keyword">...</span>
    0 0; <span class="keyword">...</span>
    30 0; <span class="keyword">...</span>
    30 30];

<span class="comment">% set of ellipsoid orientations, in format YAW PITCH ROLL</span>
orientations = [<span class="keyword">...</span>
     0   0   0; <span class="keyword">...</span>
    30   0   0; <span class="keyword">...</span