function face_test ( )

%*****************************************************************************80
%
% face_test: Example polygonal geometries for MESHFACES.
%
%  Author:
%
%    Darren Engwirda
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'FACE_TEST\n' ); 
%
%  Case 1:
%
      node = [0.0, 0.0; 1.0,0.0; 1.0,1.0; 0.0,1.0; 1.01,0.0; 1.01,1.0; 3.0,0.0; 3.0,1.0];
      edge = [1,2; 2,3; 3,4; 4,1; 2,5; 5,6; 6,3; 5,7; 7,8; 8,6];
      face{1} = [1,2,3,4];
      face{2} = [5,6,7,2];
      face{3} = [8,9,10,6];

      meshfaces ( node, edge, face );
%
%  Case 2
%

      % Geometry
      dtheta = pi/36;
      theta = (-pi:dtheta:(pi-dtheta))';
      node1 = [cos(theta), sin(theta)];
      node2 = [-2.0,-2.0; 2.0,-2.0; 2.0,2.0; -2.0, 2.0];
      edge1 = [(1:size(node1,1))',[(2:size(node1,1))'; 1]];
      edge2 = [1,2; 2,3; 3,4; 4,1];

      edge = [edge1; edge2+size(node1,1)];
      node = [node1; node2];

      face{1} = 1:size(edge1,1);
      face{2} = 1:size(edge,1);

      meshfaces(node,edge,face);
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'FACE_TEST\n' ); 
  fprintf ( 1, '  Normal end of execution.\n' );

  return
end
