function polyomino_embed_number_test ( )

%*****************************************************************************80
%
%% POLYOMINO_EMBED_NUMBER_TEST tests POLYOMINO_EMBED_NUMBER.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    25 April 2018
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'POLYOMINO_EMBED_NUMBER_TEST:\n' );
  fprintf ( 1, '  POLYOMINO_EMBED_NUMBER reports the number of ways a\n' );
  fprintf ( 1, '  fixed polyomino can be embedded in a region.\n' );

  mr = 4;
  nr = 4;
  r = [ ...
    0, 1, 1, 1; 
    1, 1, 0, 1;
    1, 1, 1, 1; 
    1, 0, 1, 1 ];

  polyomino_print ( mr, nr, r, '  The given region R:' );

  mp = 3;
  np = 2;
  p = [ ...
    0, 1;
    0, 1;
    1, 1 ];

  polyomino_print ( mp, np, p, '  The given polyomino P:' );

  number = polyomino_embed_number ( mr, nr, r, mp, np, p );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  As a fixed polyomino, P can be embedded in R in %d ways\n', number );
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'POLYOMINO_EMBED_NUMBER_TEST:\n' );
  fprintf ( 1, '  Normal end of execution.\n' );

  return
end
