step2;

x=[-3:0.1:3];
y=0*x+1;
plot3(x,y,f(x,y),'linewidth',5,'color','k');

