picture_paper.m

figure(100)
dx = [0:1:1]; dy = [0:1:1]; dz = [0:1:1];
[x,y,z] = meshgrid(dx,dy,dz);
X1_coarse = [x(:) y(:) z(:)];
T1_coarse = delaunay3(x,y,z,{'QJ'});
[X1_coarse,T1_coarse]=tetrarefine3(X1_coarse,T1_coarse);
    
dx = [0:1:1]; dy = [0:1:1]; dz = [-1:1:0];
[x,y,z] = meshgrid(dx,dy,dz);
X2_coarse = [x(:) y(:) z(:)];
T2_coarse = delaunay3(x,y,z,{'QJ'});

subplot(2,3,1);show_mesh(T1_coarse,X1_coarse),axis off
subplot(2,3,4);show_mesh(T2_coarse,X2_coarse),axis off

[X1_coarse,T1_coarse]=tetrarefine3(X1_coarse,T1_coarse);
[X2_coarse,T2_coarse]=tetrarefine3(X2_coarse,T2_coarse);
subplot(2,3,2);show_mesh(T1_coarse,X1_coarse),axis off
subplot(2,3,5);show_mesh(T2_coarse,X2_coarse),axis off

[X1_coarse,T1_coarse]=tetrarefine3(X1_coarse,T1_coarse);
[X2_coarse,T2_coarse]=tetrarefine3(X2_coarse,T2_coarse);
subplot(2,3,3);show_mesh(T1_coarse,X1_coarse),axis off
subplot(2,3,6);show_mesh(T2_coarse,X2_coarse),axis off