function uvp_gms_test ( )

%*****************************************************************************80
%
%% uvp_gms_test samples the GMS solution.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    19 August 2020
%
%  Author:
%
%    John Burkardt
%
  nu = 1.0;
  rho = 1.0;
  t = 1.0;

  xmin =-1.0;
  xmax = 1.0;
  ymin = -1.0;
  ymax = 1.0;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'uvp_gms_test\n' );
  fprintf ( 1, '  MATLAB/Octave version %s\n', version ( ) );
  fprintf ( 1, '  GMS flow\n' );
  fprintf ( 1, '  Estimate the range of velocity and pressure\n' );
  fprintf ( 1, '  %g <= X <= %g, %g <= Y <= %g\n', xmin, xmax, ymin, ymax );
  fprintf ( 1, '  T = %g\n', t );
  fprintf ( 1, '  Kinematic viscosity NU = %g\n', nu );
  fprintf ( 1, '  Fluid density RHO = %g\n', rho );

  n = 1000;
  seed = 123456789;
  [ x, seed ] = r8vec_uniform_ab ( n, xmin, xmax, seed );
  [ y, seed ] = r8vec_uniform_ab ( n, ymin, ymax, seed );

  [ u, v, p ] = uvp_gms ( nu, rho, n, x, y, t );

  fprintf ( 1, '\n' );
  fprintf ( 1, '           Minimum       Maximum\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  U:  %14.6g  %14.6g\n', min ( u ), max ( u ) );
  fprintf ( 1, '  V:  %14.6g  %14.6g\n', min ( v ), max ( v ) );
  fprintf ( 1, '  P:  %14.6g  %14.6g\n', min ( p ), max ( p ) );

  return
end
