%% TEST05 tests DTRIS2.
%
  ndim = 2;
  g_num = 9;

  g_xy = [ ...
       0.0, 0.0; ...
       0.0, 1.0; ...
       0.2, 0.5; ...
       0.3, 0.6; ...
       0.4, 0.5; ...
       0.6, 0.4; ...
       0.6, 0.5; ...
       1.0, 0.0; ...
       1.0, 1.0 ]';

  fprintf ( 1, '\n' );
  fprintf ( 1, 'TEST05\n' );
  fprintf ( 1, '  DTRIS2 computes the Delaunay triangulation of a\n' );
  fprintf ( 1, '  pointset in 2D.\n' );
%
%  Set up the Delaunay triangulation.
%
  [ tri_num, nod_tri, tnbr ] = dtris2 ( g_num, g_xy );

  triangulation_print ( g_num, tri_num, g_xy, nod_tri, tnbr );