function x = wilson_solution ( )

%*****************************************************************************80
%
%% wilson_solution() returns the WILSON solution.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    29 November 2019
%
%  Author:
%
%    John Burkardt
%
%  Output:
%
%    real X(4,1), the solution vector.
%
  x(1:4,1) = [ 1.0; ...
               2.0; ...
               3.0; ...
               4.0 ];

  return
end
