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

%*****************************************************************************80
%
%% BODEWIG_PLU returns the PLU factors of the BODEWIG matrix.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    23 March 2015
%
%  Author:
%
%    John Burkardt
%
%  Parameters:
%
%    Output, real P(4,4), L(4,4), U(4,4), the PLU factors.
%
  p = [ ...
    0.0, 0.0, 0.0, 1.0; ...
    0.0, 1.0, 0.0, 0.0; ...
    0.0, 0.0, 1.0, 0.0; ...
    1.0, 0.0, 0.0, 0.0 ];

  l = [ ...
    1.00, 0.00,              0.00,              0.00; ...
    0.25, 1.00,              0.00,              0.00; ...
    0.75, 0.647058823529412, 1.00,              0.00; ...
    0.50, 0.352941176470588, 0.531531531531532, 1.00 ];

  u = [ ...
    4.00,  5.00, -2.00,              -1.00; ...
    0.00, -4.25,  1.50,               5.25; ...
    0.00,  0.00,  6.529411764705882, -4.647058823529412; ...
    0.00,  0.00,  0.00,               5.117117117117118 ];

  return
end
