function polyomino_equal_test ( )

%*****************************************************************************80
%
%% polyomino_equal_test tests polyomino_equal.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    02 May 2020
%
%  Author:
%
%    John Burkardt
%
  timestamp ( );
  fprintf ( 1, '\n' );
  fprintf ( 1, 'polyomino_equal_test\n' );
  fprintf ( 1, '  MATLAB/Octave version %s\n', version ( ) );
  fprintf ( 1, '  polyomino_equal determines if polyominoes are equal.\n' );

  p1 = [ ...
    1, 0, 1, 1;
    1, 1, 1, 0;
    0, 1, 1, 0 ];
  [ m1, n1 ] = size ( p1 );
  polyomino_print ( m1, n1, p1, '  The polyomino P1:' );

  p2 = [ ...
    1, 0, 1, 1;
    1, 1, 1, 0;
    0, 1, 0, 0 ];
  [ m2, n2 ] = size ( p2 );
  polyomino_print ( m2, n2, p1, '  The polyomino P2:' );

  p3 = [ ...
    1, 0, 1, 1, 0;
    1, 1, 1, 0, 0;
    0, 1, 1, 0, 0 ];
  [ m3, n3 ] = size ( p3 );
  polyomino_print ( m3, n3, p3, '  The polyomino P3:' );

  p4 = [ ...
    1, 0, 1, 1;
    1, 1, 1, 0;
    0, 1, 1, 0 ];
  [ m4, n4 ] = size ( p4 );
  polyomino_print ( m4, n4, p4, '  The polyomino P4:' );

  fprintf ( 1, '\n' );
  value = polyomino_equal ( m1, n1, p1, m1, n1, p1 );
  fprintf ( 1, '  P1 == P1? %s\n', mat2str ( value ) );
  value = polyomino_equal ( m1, n1, p1, m2, n2, p2 );
  fprintf ( 1, '  P1 == P2? %s\n', mat2str ( value ) );
  value = polyomino_equal ( m1, n1, p1, m3, n3, p3 );
  fprintf ( 1, '  P1 == P3? %s\n', mat2str ( value ) );
  value = polyomino_equal ( m1, n1, p1, m4, n4, p4 );
  fprintf ( 1, '  P1 == P4? %s\n', mat2str ( value ) );
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'polyomino_equal_test\n' );
  fprintf ( 1, '  Normal end of execution.\n' );
  fprintf ( 1, '\n' );
  timestamp ( );

  return
end
