function ff(xmin,xmax,fx)
    clf;
    xx=linspace(xmin,xmax);     yy=xx;
    [x,y]=meshgrid(xx,yy);
    z=eval(fx);     z=(z<=0);
    plot3(x,y,z);     view(-45,3);
    title(fx);     xlabel('x');    ylabel('y');    zlabel('z');    grid on;    axis tight;
end

f3(-6,6,'(x-1).^2-(y-1).^2-1')