% esercizio 3
%


NP = 101; %numero di punti in cui dividere intervalli

x=linspace(-10,10, NP);
y=linspace(-10,10, NP);

% creo rete
[ X , Y ] = meshgrid(x,y);

% per il calcolo di Z, uso una funzione definita a parte
Z = fz_es3( X, Y );

figure(100)
clf
surf(X, Y, Z)
xlabel('x');
ylabel('y');
zlabel('z(x,y)');

figure(101)
clf
contour(X, Y, Z)
xlabel('x');
ylabel('y');
axis square

