y=[-3:0.1:3];
x=0*y+1;
plot3(x,y,f(x,y),'linewidth',5,'color','k');
axis([-3 3 -3 3 -4 4]);
hold on
xlabel('X','fontweight','bold','fontsize',12);
ylabel('Y','fontweight','bold','fontsize',12);
zlabel('z','fontweight','bold','fontsize',12);
plot3(1,1,0,'r+','linewidth',2)