function make_flow_example(n)
[x,y,tri]=make_rand_tri(n);
M=make_adj_matrix(tri);
B=find_boundary(tri,M);
% plot triangulation alone
figure; hold on; axis equal; axis off;
title(['A random trianglulation, n = ', num2str(n)]);
plot_tri(x,y,tri)
% plot triangulation and boundary
figure; hold on; axis equal; axis off;
title('The boundary of the triangulation');
plot_tri(x,y,tri)
M=make_adj_matrix(tri);
B=find_boundary(tri,M);
plot_tri_boundary(x,y,B);
% plot triangulation, boundary and central region
figure; hold on; axis equal; axis off;
title('The thin central regions');
plot_tri(x,y,tri)
plot_tri_boundary(x,y,B);
radii=find_radii(x,y,tri);
plot_central_regions2(x,y,tri,radii);
% plot triangulation, boundary, thin centrals and flows
figure; hold on; axis equal; axis off;
title('The thin region flow');
plot_tri(x,y,tri)
plot_tri_boundary(x,y,B);
radii=find_radii(x,y,tri);
plot_central_regions2(x,y,tri,radii);
propogate2(x,y,tri,radii,B)
% plot triangulation, boundary, thick centrals and flows
figure; hold on; axis equal; axis off;
title('The thick region flow');
plot_tri(x,y,tri)
plot_tri_boundary(x,y,B);
radii=find_radii2(x,y,tri);
plot_central_regions2(x,y,tri,radii);
propogate4(x,y,tri,radii,B)
return