function [ node_xyz, tet_node ] = tet_mesh_order4_example_set ( node_num, ...
  tet_num )

%*****************************************************************************80
%
%% TET_MESH_ORDER4_EXAMPLE_SET sets an example linear tet mesh.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license. 
%
%  Modified:
%
%    19 August 2009
%
%  Author:
%
%    John Burkardt
%
%  Parameters:
%
%    Input, integer NODE_NUM, the number of nodes.
%
%    Input, integer TET_NUM, the number of tetrahedrons.
%
%    Output, real NODE_XYZ(3,NODE_NUM), the node coordinates.
%
%    Output, integer TET_NODE(4,TET_NUM), the nodes 
%    forming each tet.
%
  tet_node(1:4,1:tet_num) = [ ...
         1,   2,   4,  10; ...
         2,   4,   5,  10; ...
         2,   5,  10,  11; ...
         2,   3,   5,  11; ...
         4,   5,  10,  13; ...
         3,   5,   6,  11; ...
         5,  10,  11,  13; ...
         4,   5,   7,  13; ...
         5,   6,   8,  14; ...
         5,   7,   8,  13; ...
         6,   8,   9,  14; ...
        11,  13,  14,  19; ...
        12,  14,  15,  20; ...
         3,   6,  11,  12; ...
         5,   6,  11,  14; ...
         6,   9,  14,  15; ...
         6,  11,  12,  14; ...
         6,  12,  14,  15; ...
         7,   8,  13,  16; ...
         5,   8,  13,  14; ...
        10,  11,  13,  19; ...
         8,   9,  14,  17; ...
        11,  12,  14,  20; ...
         5,  11,  13,  14; ...
         8,  13,  14,  16; ...
         9,  14,  15,  17; ...
        13,  14,  16,  22; ...
         8,  14,  16,  17; ...
        14,  15,  17,  23; ...
        14,  16,  17,  22; ...
         9,  15,  17,  18; ...
        15,  17,  18,  23; ...
        14,  17,  22,  23; ...
        13,  14,  19,  22; ...
        11,  14,  19,  20; ...
        14,  15,  20,  23; ...
        15,  20,  21,  23; ...
        21,  23,  24,  29; ...
        20,  22,  23,  28; ...
        14,  19,  20,  22; ...
        15,  18,  23,  24; ...
        12,  15,  20,  21; ...
        15,  21,  23,  24; ...
        16,  17,  22,  25; ...
        19,  20,  22,  28; ...
        17,  18,  23,  26; ...
        20,  21,  23,  29; ...
        14,  20,  22,  23; ...
        17,  22,  23,  25; ...
        18,  23,  24,  26; ...
        22,  23,  25,  31; ...
        17,  23,  25,  26; ...
        23,  24,  26,  32; ...
        23,  25,  26,  31; ...
        18,  24,  26,  27; ...
        24,  26,  27,  32; ...
        23,  26,  31,  32; ...
        22,  23,  28,  31; ...
        20,  23,  28,  29; ...
        23,  24,  29,  32; ...
        24,  29,  30,  32; ...
        30,  32,  33,  38; ...
        29,  31,  32,  37; ...
        23,  28,  29,  31; ...
        24,  27,  32,  33; ...
        21,  24,  29,  30; ...
        24,  30,  32,  33; ...
        25,  26,  31,  34; ...
        28,  29,  31,  37; ...
        26,  27,  32,  35; ...
        29,  30,  32,  38; ...
        23,  29,  31,  32; ...
        26,  31,  32,  34; ...
        27,  32,  33,  35; ...
        31,  32,  34,  40; ...
        26,  32,  34,  35; ...
        32,  33,  35,  41; ...
        32,  34,  35,  40; ...
        27,  33,  35,  36; ...
        33,  35,  36,  41; ...
        32,  35,  40,  41; ...
        31,  32,  37,  40; ...
        29,  32,  37,  38; ...
        32,  33,  38,  41; ...
        33,  38,  39,  41; ...
        39,  41,  42,  47; ...
        38,  40,  41,  46; ...
        32,  37,  38,  40; ...
        33,  36,  41,  42; ...
        30,  33,  38,  39; ...
        33,  39,  41,  42; ...
        34,  35,  40,  43; ...
        37,  38,  40,  46; ...
        35,  36,  41,  44; ...
        38,  39,  41,  47; ...
        32,  38,  40,  41; ...
        35,  40,  41,  43; ...
        36,  41,  42,  44; ...
        40,  41,  43,  49; ...
        35,  41,  43,  44; ...
        41,  42,  44,  50; ...
        41,  43,  44,  49; ...
        36,  42,  44,  45; ...
        42,  44,  45,  50; ...
        41,  44,  49,  50; ...
        40,  41,  46,  49; ...
        38,  41,  46,  47; ...
        41,  42,  47,  50; ...
        42,  47,  48,  50; ...
        48,  50,  51,  56; ...
        47,  49,  50,  55; ...
        41,  46,  47,  49; ...
        42,  45,  50,  51; ...
        39,  42,  47,  48; ...
        42,  48,  50,  51; ...
        43,  44,  49,  52; ...
        46,  47,  49,  55; ...
        44,  45,  50,  53; ...
        47,  48,  50,  56; ...
        41,  47,  49,  50; ...
        44,  49,  50,  52; ...
        45,  50,  51,  53; ...
        49,  50,  52,  58; ...
        44,  50,  52,  53; ...
        50,  51,  53,  59; ...
        50,  52,  53,  58; ...
        45,  51,  53,  54; ...
        51,  53,  54,  59; ...
        50,  53,  58,  59; ...
        49,  50,  55,  58; ...
        47,  50,  55,  56; ...
        50,  51,  56,  59; ...
        51,  56,  57,  59; ...
        50,  55,  56,  58; ...
        51,  54,  59,  60; ...
        48,  51,  56,  57; ...
        51,  57,  59,  60; ...
        52,  53,  58,  61; ...
        53,  54,  59,  62; ...
        50,  56,  58,  59; ...
        53,  58,  59,  61; ...
        54,  59,  60,  62; ...
        53,  59,  61,  62; ...
        54,  60,  62,  63 ]';

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

  return
end