% zeichnet eine Funktion z = f(x,y) als eine aus Dreiecken bestehende Oberflaeche

figure 

t=linspace(-1.5*pi,1.5*pi,25);
[x,y]=meshgrid(t,t);

z=cos(x+cos(y));
z(z<0)=0;                                    % negative z - Werte werden 0 gesetzt  

tri = delaunay(x,y);                         % fuehrt mit dem x,y Netz eine Delaunay Triangulation durch

h=trisurf(tri,x,y,z);                        % Die Oberflaeche z=f(x,y) wird mit Dreiecken gezeichnet und ausgefuellt
shading interp
set(h,'edgecolor','k')
xlabel('x')
ylabel('y');
zlabel('z=cos((x+cos(y))')
title('trisurf(tri,x,y,z)','fontsize',16);

view(-25,65)
colormap summer