clear

% grafici tridimensionali

% vogliamo fare il grafico di
%
% z = x(1-x)y(1-y) 
%
% sul rettangolo [0,1] X [0,1]

n = 101;
m = 101;
x = linspace(0,2*pi,n);
y = linspace(0,2*pi,m);

% creiamo una 'mesh'
% vediamo cosa contiene [X,Y] omettendo il ;
[X,Y] = meshgrid(x,y)

% calcoliamo la funzione
%Z = X .* (1-X) .* Y .* (1-Y);
Z = sin( X .^2 + Y .^2 ) ./ ( X .^2 + Y .^2 );

figure(1)
% ripuliamo la figura da altri grafici
clf
% grafico a superfice (approsimata con tasselli)
surf(X,Y,Z);

figure(2)
clf
% grafico a curve di livello
contour(X,Y, Z);
axis square

figure(3)
clf
% grafico a 'griglia'
mesh(X,Y,Z);

figure(4)
% ripuliamo la figura da altri grafici
clf
% grafico a superfice (approsimata con tasselli)
% con curve di livello
surfc(X,Y,Z);

