[X Y]=meshgrid(-3:0.5:3);
Z=f(X,Y);
mesh(X,Y,Z);
xlabel('X','fontweight','bold','fontsize',12);
ylabel('Y','fontweight','bold','fontsize',12);
zlabel('z','fontweight','bold','fontsize',12);
hold on
axis([-3 3 -3 3 -4 4])
plot(1,1,'r+','linewidth',2)
