function fem2d_pack_test16 ( )

%*****************************************************************************80
%
%% TEST16 tests REFERENCE_TO_PHYSICAL_T6.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    14 February 2003
%
%  Author:
%
%    John Burkardt
%
  n = 16;

  ref = [ ...
    0.00, 0.00; ...
    1.00, 0.00; ...
    0.00, 1.00; ...
    0.50, 0.00; ...
    0.50, 0.50; ...
    0.00, 0.50; ...
    0.25, 0.75; ...
    0.75, 0.25; ...
    0.40, 0.10; ...
    0.30, 0.20; ...
    0.20, 0.30; ...
    0.10, 0.40; ...
    0.10, 0.10; ...
    0.20, 0.20; ...
    0.30, 0.30; ...
    0.40, 0.40 ]';
  t = [ ...
    0.0, 0.0; ...
    2.0, 0.0; ...
    0.0, 4.0; ...
    1.0, 0.0; ...
    1.0, 1.0; ...
    0.0, 2.0 ]';

  fprintf ( 1, '\n' );
  fprintf ( 1, 'TEST16\n' );
  fprintf ( 1, '  For an order 6 triangle,\n' );
  fprintf ( 1, '  REFERENCE_TO_PHYSICAL_T6 maps a reference point to\n' );
  fprintf ( 1, '    a physical point.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '      XSI     ETA  ==>  X	Y\n' );
  fprintf ( 1, '\n' );

  phy = reference_to_physical_t6 ( t, n, ref );

  for j = 1 : n
    fprintf ( 1, '  %8f  %8f  %8f  %8f\n', ref(1:2,j), phy(1:2,j) );
  end 

  return
end
