function gray_scott_movie_test ( )

%*****************************************************************************80
%
%% gray_scott_movie_test() tests gray_scott_movie().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    07 December 2020
%
%  Author:
%
%    John Burkardt
%
  addpath ( '../gray_scott_movie' );

  timestamp ( );
  fprintf ( 1, '\n' );
  fprintf ( 1, 'gray_scott_movie_test():\n' );
  fprintf ( 1, '  gray_scott_movie() solves the 2D Gray-Scott \n' );
  fprintf ( 1, '  reaction-diffusion equation, and makes a movie.\n' );
%
%  Set the grid size.
%
  n = 200;
%
%  Use the same grid size for x and y.
%
  nx = n;
  ny = n;
%
%  Assume the grid spacing is 1.
%
  dx = 1.0;
  dy = 1.0;
%
%  Set the diffusion coefficients for U and V.
%
  du = 0.3;
  dv = 0.15;
%
%  There are several interesting choices for the feed and kill rates.
%
%  Rings.
%
  if ( false )
    f = 0.01;
    k = 0.1;
%
%  Splits.
%
  elseif ( false )
    f = 0.0367;
    k = 0.0649;
%
%  Coral.
%
  elseif ( true )
    f = 0.0545;
    k = 0.0620;
%
%  Bubbles.
%
  elseif ( false )
    f = 0.012;
    k = 0.05;
%
%  Waves.
%
  elseif ( false )
    f = 0.025;
    k = 0.05;
  end
%
%  Set the time step.
%
  dt = 1.0;
%
%  Choose the initial condition.
%
  if ( true )
    initial_condition = 'wavefront';
  elseif ( false )
    initial_condition = 'bar';
  else
    initial_condition = 'square';
  end

  gray_scott_movie ( nx, ny, f, k, dx, dy, dt, du, dv, initial_condition );
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'gray_scott_movie_test():\n' );
  fprintf ( 1, '  Normal end of execution.\n' );
  timestamp ( );

  rmpath ( '../gray_scott_movie' );

  return
end
function timestamp ( )

%*****************************************************************************80
%
%% timestamp() prints the current YMDHMS date as a timestamp.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    14 February 2003
%
%  Author:
%
%    John Burkardt
%
  t = now;
  c = datevec ( t );
  s = datestr ( c, 0 );
  fprintf ( 1, '%s\n', s );

  return
end

