function [ p, l, u ] = a123_plu ( )

%*****************************************************************************80
%
%% A123_PLU returns the PLU factors of the A123 matrix.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    09 March 2015
%
%  Author:
%
%    John Burkardt
%
%  Parameters:
%
%    Output, real P(3,3), L(3,3), U(3,3), the PLU factors.
%
  p = [ ...
     0,     1,     0; ...
     0,     0,     1; ...
     1,     0,     0 ];

  l = [ ...
     1.000000000000000,                   0,                   0; ...
     0.142857142857143,   1.000000000000000,                   0; ...
     0.571428571428571,   0.500000000000000,   1.000000000000000 ];

  u = [ ...
   7.000000000000000,   8.000000000000000,   9.000000000000000; ...
                   0,   0.857142857142857,   1.714285714285714; ...
                   0,                   0,   0.000000000000000 ];

  return
end