f = @(x,y) sqrt( x.^2 + x.*y + y.^2 ) ;
n = 500 ;
vx = linspace(-1, 1, n) ;
[x, y] = meshgrid(vx, vx) ;
z = f(x,y) ;
livelli = 15 ;
contour(x, y, z, livelli) ;
axis equal ;
xlabel('x') ;
ylabel('y') ;