function [ u, v, p ] = uvp_poiseuille ( nu, rho, n, x, y, t )

%*****************************************************************************80
%
%% uvp_poiseuille() returns the Poiseuille solution.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    22 July 2015
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    real NU, the kinematic viscosity.
%
%    real RHO, the fluid density.
%
%    integer N, the number of nodes.
%
%    real X(N), Y(N), the coordinates of nodes.
%
%    real T, the current time.
%
%  Output:
%
%    real U(N), V(N), the X and Y velocity.
%
%    real P(N), the pressure.
%
  u(1:n) = 1.0 - y(1:n) .^ 2;

  v(1:n) = 0.0;

  p(1:n) = -2.0 * rho * nu * x(1:n);

  return
end

