function a = bodewig_matrix ( )

%*****************************************************************************80
%
%% bodewig_matrix returns the BODEWIG matrix.
%
%  Example:
%
%    2   1   3   4
%    1  -3   1   5
%    3   1   6  -2
%    4   5  -2  -1
%
%  Properties:
%
%    A is symmetric: A' = A.
%
%    Because it is symmetric, it is normal: A*A'=A'*A.
%
%    Because it is normal, it is diagonalizable.
%
%    A is not positive definite.
%
%    A is not diagonally dominant.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    28 September 2007
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    Joan Westlake,
%    Test Matrix A27,
%    A Handbook of Numerical Matrix Inversion and Solution of Linear Equations,
%    John Wiley, 1968.
%
%  Output:
%
%    real A(4,4), the matrix.
%
  a(1:4,1:4) = [ ...
    2.0,  1.0,  3.0,  4.0; ...
    1.0, -3.0,  1.0,  5.0; ...
    3.0,  1.0,  6.0, -2.0; ...
    4.0,  5.0, -2.0, -1.0 ];

  return
end
