%% Flaeche im Raum, kartesisches Grundgitter
clc; clf; clear;
% 1D-Koordinaten des Gitters
range=4;
x = linspace(-range , range, 30);
y = linspace(-range,  range, 30);

% 2D-Koordinaten des gesamten Gitters
[ xx, yy ] = meshgrid(x, y);

% z-werte der Funktion ausrechnen

zz = sin(pi*sqrt(xx.^2+yy.^2))./sqrt(xx.^2+yy.^2);

% Darstellen
subplot(1,2,1);
surf(xx, yy, zz);

title('via kartesische Koord.')
xlim([-range, range]); ylim([-range, range]); zlim([-range, range]);
%axis equal;


%% Flaeche im Raum, Polarkoordinatengitter

% 1D-Koordinaten des Gitters
phi = linspace(0 , 2*pi, 30);     % Winkel
rad = linspace(0,  sqrt(range^2+range^2), 20);        % Radius

% 2D-Polarkoordinaten des gesamten Gitters
[ pp, rad ] = meshgrid(phi, rad);

% 2D kart. Koordinaten des gesamten Gitters
[xx, yy] = pol2cart(pp,rad);

% z-werte der Funktion ausrechnen

zz = sin(pi*rad)./rad;

% Darstellen
subplot(1,2,2); 
surf(xx, yy, zz);
title('via Polarkoord.')
xlim([-range, range]); ylim([-range, range]); zlim([-range, range]);
%axis equal;