>> [x,y]=meshgrid(-10:0.1:10);
z=2*sqrt((x.^2+y.^2)/9-1);
z((x.^2+y.^2)/9-1