% script per sperimentare la rappresentazione tridimensionale

% numero di punti in cui viene diviso l'intervallo
% per la variabile x e la variabile y
nx = 25;
ny = 25

x = linspace(0,1,nx);
y = linspace(0,1,ny);

% costruzione del reticolo
[ X Y ] = meshgrid(x,y);

% calcolo della funzione sul reticolo
Z = X .* (1-X) .* Y .* (1-Y);


% grafico della supeficie 3D
figure(103)
surf(X,Y,Z)
xlabel(' x ')
ylabel(' y ')
zlabel('f(x,y) = x*(1-x)*y*(1-y)')

% grafico delle curve di livello
figure(104)
contour(X,Y,Z)
xlabel(' x ')
ylabel(' y ')



