]=meshacylinder(c1,c2,.2); >>>>>>> [cnode,cface]=meshcheckrepair(cnode,cface,'dup'); >>>>>>> >>>

done via surfboolean('decouple'). Qianqian or opt(1).radbound=4; opt(2).radbound=2; opt(1).side='lower'; opt(2).side='lower'; opt(1).A=diag

mergesurf /surfboolean (if they do intersect) to properly merge the two domains. if you want to create a cavity, you should use the "holes" list in the s2m. here is an example

then use surfboolean > to combine with a cylindrical mesh to create an inclusion. > > in your case, if you want to combine with a sphere, you can simply replace > the meshacylinder

the 'surfboolean' failed in your script is >>> because the bottom of the cylinders are overlapped with that of the >>> bounding box. So I extrude the cylinders

]=meshacylinder([5 15 5],[5 -15 5],2.5,1,2);/ > / > / > /% Create mesh for a 2.5 mm radius cylinder:/ > /[no3,fa3,el3]=meshacylinder([-5 15 -5],[-5 -15 -5],2.5,1,2

problems with surfboolean, I followed your advise to use surf2vol() to convert each component of three surfaces into a volumetric image. Then I combined these volumetric masks into