function uvwp_poiseuille_test ( )

%*****************************************************************************80
%
%% uvwp_poiseuille_test tests uvwp_poiseuille.
%
%  Location:
%
%    http://people.sc.fsu.edu/~jburkardt/m_src/navier_stokes_3d_exact/uvwp_poiseuille_test.m
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    12 January 2020
%
%  Author:
%
%    John Burkardt
%
  nu = 0.25;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'uvwp_poiseuille_test\n' );
  fprintf ( 1, '  uvwp_poiseuille evaluates the poiseuille solution.\n' );
  fprintf ( 1, '  Estimate the range of velocity and pressure\n' );
  fprintf ( 1, '  at the initial time T = 0, in a region that is an \n' );
  fprintf ( 1, '  infinite cylinder along the x-axis with "radius" 1.0.\n' );
  fprintf ( 1, '  Density nu = %g\n', nu );

  n = 1000;
  x = 10.0 * rand ( n, 1 );
  r = sqrt ( rand ( n, 1 ) );
  t = 2.0 * pi * rand ( n, 1 );
  y = r .* cos ( t );
  z = r .* sin ( t );
  t = 0.0;

  [ u, v, w, p ] = uvwp_poiseuille ( nu, n, x, y, z, 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, '  W:  %14.6g  %14.6g\n', min ( w ), max ( w ) );
  fprintf ( 1, '  P:  %14.6g  %14.6g\n', min ( p ), max ( p ) );

  return
end
