#==== intersection, réunion, différence de deux surfaces implicites ==========
with(plots):
S1:=x^2+y^2+z^2-4:  S2:=(x-1)^2+y^2-1: #--- sphère et cylindre
a,b,c:=3,2,3;  ng:=100:
preci:=grid=[ng,ng,ng]:
display(array([implicitplot3d(max(S1,S2),x=-a..a,y=-b..b,z=-c..c,preci),
               implicitplot3d(min(S1,S2),x=-a..a,y=-b..b,z=-c..c,preci),
               implicitplot3d(min(-S1,S2),x=-a..a,y=-b..b,z=-c..c,preci)
                    ]),style=patchnogrid,color=cyan,scaling=constrained,axes=none);