function hw7p1 ( )

  v = [ ...
    1.0, 10.0; ...
    6.0, 10.0; ...
    6.0, 14.0; ...
    1.0, 14.0; ...
    4.0, 11.0; ...
    5.0, 11.0; ...
    5.0, 12.0; ...
    4.0, 12.0 ];

  e = [ ...
    1, 2; ...
    2, 3; ...
    3, 4; ...
    4, 1; ...
    8, 7; ...
    7, 6; ...
    6, 5; ...
    5, 8 ];

  hdata = [];
  hdata.hmax = 0.5;

  [ p, t ] = mesh2d ( v, e, hdata );

  p

  t

  mesh2d_to_xml ( 'rectangle.xml', p, t );
  
  return
end
