function x = solution ( )

%*****************************************************************************80
%
%% solution() returns the solution.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    03 January 2020
%
%  Author:
%
%    John Burkardt
%
%  Output:
%
%    real x(m,k), the solution vectors.
%
  global Ab
  global k
  global m
  global n

  I = eye ( m, n );

  norm_fro = norm ( Ab(1:m,1:n) - I, 'fro' );

  if ( eps < norm_fro )
    fprintf ( 1, 'solution - Fatal error!\n' );
    fprintf ( 1, '  Augmented matrix Ab does not start with the identity.\n' );
    return
  end

  x = Ab(1:m,n+1:n+k);

  return
end

