function [ node_xyz, tetra_node ] = tet_mesh_order10_example_set ( ...
  node_num, tetra_num )

%*****************************************************************************80
%
%% TET_MESH_ORDER10_EXAMPLE_SET sets an example tet mesh.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    03 August 2009
%
%  Author:
%
%    John Burkardt
%
%  Parameters:
%
%    Input, integer NODE_NUM, the number of nodes.
%
%    Input, integer TETRA_NUM, the number of tetrahedrons.
%
%    Output, real NODE_XYZ(3,NODE_NUM), the node coordinates.
%
%    Output, integer TETRA_NODE(10,TETRA_NUM), the nodes forming each tet.
%
  tetra_node(1:10,1:tetra_num) = [ ...
    4,   3,   5,   1,  16,  19,  17,  11,  10,  12; ...
    4,   2,   5,   1,  13,  19,  14,  11,   9,  12; ...
    4,   7,   3,   5,  21,  16,  18,  19,  24,  17; ...
    4,   7,   8,   5,  21,  22,  27,  19,  24,  25; ...
    4,   6,   2,   5,  20,  13,  15,  19,  23,  14; ...
    4,   6,   8,   5,  20,  22,  26,  19,  23,  25 ]';

  node_xyz(1:3,1:node_num) = [ ...
   0.0,  0.0,  0.0; ...   
   0.0,  0.0,  1.0; ...    
   0.0,  1.0,  0.0; ...    
   0.0,  1.0,  1.0; ...    
   1.0,  0.0,  0.0; ...  
   1.0,  0.0,  1.0; ...
   1.0,  1.0,  0.0; ...    
   1.0,  1.0,  1.0; ...    
   0.0,  0.0,  0.5; ...    
   0.0,  0.5,  0.0; ...   
   0.0,  0.5,  0.5; ...   
   0.5,  0.0,  0.0; ...   
   0.0,  0.5,  1.0; ...    
   0.5,  0.0,  0.5; ...    
   0.5,  0.0,  1.0; ...   
   0.0,  1.0,  0.5; ...    
   0.5,  0.5,  0.0; ...   
   0.5,  1.0,  0.0; ...   
   0.5,  0.5,  0.5; ...    
   0.5,  0.5,  1.0; ...   
   0.5,  1.0,  0.5; ...    
   0.5,  1.0,  1.0; ...    
   1.0,  0.0,  0.5; ...    
   1.0,  0.5,  0.0; ...    
   1.0,  0.5,  0.5; ...    
   1.0,  0.5,  1.0; ...   
   1.0,  1.0,  0.5 ]';

  return
end
