% Daten werden in Form eines 3D Contour Plots dargestellt

x=linspace(-2*pi,2*pi,50);
y=linspace(-pi,pi,50);

[xx,yy]=meshgrid(x,y);

zz=cos(xx)+sin(yy);

figure

subplot(2,2,1)
contour3(zz,20);                            % erstellt aus den Daten der Matrix zz einen 3D Contourplot mit
                                            % 20, in Richtung z aequidistanten Hoehenlinien. Auf der x- und y-
                                            % Achse sind die 1. und 2. Indizes von zz aufgetragen
view(-15,30)                                % guenstige Waghl des Blickwinkels

xlabel('1. Index von z');
ylabel('2. Index von z');
zlabel('z');
title('contour3(z,n)','fontsize',16)
axis square

subplot(2,2,2)

contour3(xx,yy,zz,30);                      % 3D Contourplot mit 30 Hoehenlinien der Funktion zz=f(xx,yy)
xlabel('x');
ylabel('y');
zlabel('z');
title('contour3(x,y,z,n)','fontsize',16)
axis square

subplot(2,2,4)

v=[-1.5,-0.5,0.5,1.5];                      % Vektor der darzustellenden Hoehenlinien
[c,h]=contour3(xx,yy,zz,v);                 % plottet length(v) Contourlinien fuer die Werte von v und 
                                            % gibt Handle c (beinhaltet Koordinaten der Linien) und
                                            % Handle h (beinhaltet Informationen wie z. B. Farbe aller Linien)
set(h,'linewidth',3)                        % Allen Linien wird die Dicke 3 zugewiesen
clabel(c,h,'fontsize',12);                  % Beschriftung der Contourlinien
view(-30,30)
axis square

xlabel('x');
ylabel('y');
zlabel('z');
title('contour3(x,y,z,v)','fontsize',16)

set(gcf,'color','w')