[X,Y] = meshgrid(-8:.5:0); 
R = (X.^4 + Y.^4).^(1/3)+sin(X)+cos(Y)+eps;
Z = -(R.^2);
figure('OuterPosition',[100, 100, 650, 500])
%surf(X,Y,Z,'EdgeColor','flat')
waterfall(X,Y,Z)
%colormap(winter)
hold on
contour3(X,Y,Z)
%camlight
%lighting phong

%[X1,Y1] = meshgrid(-8:2:0);
%R1 = (X1.^4 + Y1.^4).^(1/3) + eps;
%Z1 = -(R1.^2);
%plot3(X1,Y1,Z1,'o','MarkerEdgeColor','k','MarkerFaceColor','k','MarkerSize',5)