% script per esercizio 4
%


% numero di punti su asse x
nx = 50;
NX = 2*nx ;
% vettore x
x = linspace( -0.01, 0.01, NX );

% numero di punti su asse y
ny = 50;
NY = 2*ny ;
% vettore x
y = linspace( -0.01, 0.01, NY );


% creazione automatica della griglia

[ X, Y ] = meshgrid( x, y )

% calcolo della funzione f

Z = exp( -((X.^2)/4 + Y.^2) ) ./ ((X.^2)/4 + Y.^2) ;

% grafico di superficie
figure(1)
clf
surf(X,Y,Z)

% grafico a curve di libello
figure(2)
clf
contour(X,Y,Z)


