function xy = xy2 ( )

%*****************************************************************************80
%
%% XY2 returns XY coordinates for the nodes in example 2.
%
%  Discussion:
%
%    This matrix appears on page 35 of the reference.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    04 September 2015
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    John MacCormick,
%    Nine Algorithms That Changed the Future: The Ingenious Ideas that 
%    Drive Today's Computers,
%    Princeton University Press,
%    ISBN-13: 978-0691158198.
%
%  Parameters:
%
%    Output, integer XY(16,2), the XY coordinates for nodes.
%
  xy = zeros ( 16, 2 );

  xy(1,1) = 1.0;
  xy(1,2) = 1.0;

  xy(2,1) = 1.0;
  xy(2,2) = 3.0;

  xy(3,1) = 2.0;
  xy(3,2) = 3.0;

  xy(4,1) = 1.0;
  xy(4,2) = 5.0;

  xy(5,1) = 3.0;
  xy(5,2) = 5.0;

  xy(6,1) = 3.0;
  xy(6,2) = 3.0;

  xy(7,1) = 4.0;
  xy(7,2) = 2.0;

  xy(8,1) = 4.0;
  xy(8,2) = 5.0;

  xy(9,1) = 5.0;
  xy(9,2) = 5.0;

  xy(10,1) = 5.0;
  xy(10,2) = 4.0;

  xy(11,1) = 8.0;
  xy(11,2) = 5.0;

  xy(12,1) = 5.0;
  xy(12,2) = 3.0;

  xy(13,1) = 6.0;
  xy(13,2) = 3.0;

  xy(14,1) = 7.0;
  xy(14,2) = 3.0;

  xy(15,1) = 8.0;
  xy(15,2) = 1.0;

  xy(16,1) = 3.0;
  xy(16,2) = 1.0;

  return
end
