function [ u, v, w, p ] = uvwp_poiseuille ( nu, n, x, y, z, t )

%*****************************************************************************80
%
%% uvwp_poiseuille(): solution of the Poisseuille exact Navier Stokes solution.
%
%  Discussion:
%
%    The flow is in the x direction.  The velocity is zero along the
%    cylinder around the x axis with radius 1, which can be used to define
%    a finite flow region.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    12 January 2020
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    real NU, the viscosity.
%
%    integer N, the number of evaluation points.
%
%    real X(N), Y(N), Z(N), the coordinates of the points.
%
%    real T(N), the time coordinates.
%
%  Output:
%
%    real U(N), V(N), W(N), P(N), the solution.
%
  u = 1.0 - y.^2 - z.^2;
  v = zeros(n,1);
  w = zeros(n,1);
  p = -4.0 * x;

  return
end

