function [ l, d ] = bodewig_ldlt ( )

%*****************************************************************************80
%
%% bodewig_ldlt() returns the LDLT factors of the BODEWIG matrix.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    05 December 2019
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer N, the number of rows and columns of A.
%
%  Output:
%
%    real L(4,4), D(4,4), the factors.
%
  l = [ ...
   1.000000000000000,   0.000000000000000,   0.000000000000000,   0.000000000000000;
   0.500000000000000,   1.000000000000000,   0.000000000000000,   0.000000000000000;
   1.500000000000000,   0.142857142857143,   1.000000000000000,   0.000000000000000;
   2.000000000000000,  -0.857142857142857,  -5.363636363636365,   1.000000000000000 ];

  d = diag ( [ 2.0, -3.5, 1.571428571428571, -51.636363636363654 ] );

  return
end

