%Draw the landscape
function drawlandscape(NX,NY,Z)

x = 1:NX;
y = 1:NY;
[X,Y] = meshgrid(y,x); % strange the y goes first !!!
axis manual
mesh(X,Y,Z,'EdgeColor','black');
surf(X,Y,Z);
xlabel('x');
ylabel('y');
view([+34.5 0]);
%axis([-100 100 -100 100])