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

figure 

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

z=cos(x).*cos(y);
z(z<0)=nan;                                  % jene z mit z<0 werden mit einem nicht 
                                             % numerischen Wert (not a number) besetzt und im Plot weggelassen

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

trimesh(tri,x,y,z)                           % Die Oberflaeche z=f(x,y) wird mit Dreiecken gezeichnet
xlabel('x')
ylabel('y');
zlabel('z=cos(x).*cos(y)')
title('trimesh(tri,x,y,z)','fontsize',16);