        hh(ii) = patch(isosurface(xeva, yeva, zeva, reshape(fstar,size(xeva)), levels(ii)));
        set(hh(ii), <span class="string">'Facecolor'</span>, cmap(ii,:), <span class="string">'Edgecolor'</span>, <span class="string">'none'</span>, <span class="string">'facealpha'</span>, (1-(5-abs(levels(ii)))/5));
        axis([-3 3 -3 3 -3 3])
        caxis([-4,5]);
        xlabel(<span class="string">'x'</span>);ylabel(<span class="string">'y'</span>);zlabel(<span class="string">'z'</span>)
        colorbar
 