)


figure;
for i = 1:5
    fstar = GP.sampleprior(xstar);
    clf
    for ii = 1:5
        camlight
        lighting gouraud
        hh(ii) = patch(isosurface(xeva, yeva, zeva, reshape(fstar,size(xeva)), levels(ii)));
        set(hh(ii), 'Facecolor', cmap(ii,:), 'Edgecolor', 'none', 'facealpha', (1-(5-abs(levels(ii)))/5));
        axis([-3 3 -3 3 -3 3])
        caxis([-4,5]);
        xlabel('x');ylabel('y');zlabel('z')
        colorbar
    end
    title('Prior Samples')
    pause(0.5)
end


 
 figure;
for 