function x = wilk04_solution ( )

%*****************************************************************************80
%
%% WILK04_SOLUTION returns the solution of the WILK04 linear system.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    28 November 2007
%
%  Author:
%
%    John Burkardt
%
%  Parameters:
%
%    Output, real X(4,1), the solution of the system.
%
  x(1:4,1) = [ ...
   -9.061709180193406e+15, ...
    9.456494826647572e+11, ...
   -8.311117178175363e+07, ...
    1.102484908044364e+04 ]';

  return
end