clear;clc;[x,y,z]=meshgrid(linspace(-1.1,1.1));v=max(cat(4,x.^2+y.^2-1,z-1,-z-1),[],4);isosurface(x,y,z,v,0,z);axis equal
