 Number of function calls:	8906
  Total final loss:        	39.027
</pre><h2>Visualise GP Outputs<a name="6"></a></h2><pre class="codeinput"><span class="comment">% Build a colormap</span>
cmap = jet(5);
levels = linspace(-2.5,3.5,5);

<span class="comment">% Generate isosurfaces</span>
figure;
<span class="keyword">for</span> ii = 1:5
    camlight
    lighting <span class="string">gouraud</span>
    hh(ii) = patch(isosurface(xeva, yeva, zeva, reshape(mf,size(xeva)), levels(ii)));
    set(hh(ii), <span clas