c0=meshcentroid(node,elem(:,1:4));
plotmesh(node(:,1:3),[elem(:,1:4) c0(:,3).*c0(:,2)],'facealpha','0.7')
plotmesh(node(:,[1 2 3 2]),elem(:,1:4),'facealpha','0.7')
[nobbx,fcbbx]=meshabox([40 0 0], [160, 120, 60], 10);
[nosp,fcsp]=meshasphere(s0, 5, 1);
[no,fc]=mergemesh(nobbx, fcbbx, nosp, fcsp);
[cfg0.node, cfg0.elem]=s2m(no,fc(:,1:3),1,40,'tetgen',[41 1 1;s0]);
cfg0.seg=cfg0.elem(:,5);
plotmesh(cfg0.node,[cfg0.elem,cfg0.seg],'z<20')
view([0 0 1])
--
-Chuanyue
You received this message because you are subscribed to the Google Groups "iso2mesh-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to iso2mesh-user...@googlegroups.com.
To post to this group, send email to iso2mes...@googlegroups.com.
Visit this group at https://groups.google.com/group/iso2mesh-users.
For more options, visit https://groups.google.com/d/optout.