function errors_test14 ( )

%*****************************************************************************80
%
%% errors_test14()
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    22 August 2017
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, '+=============================================+\n' );
  fprintf ( 1, '|                                             |\n' );
  fprintf ( 1, '|  Exercise 14                                |\n' );
  fprintf ( 1, '|                                             |\n' );
  fprintf ( 1, '|    A = [ 64079  57314 ]                     |\n' );
  fprintf ( 1, '|        [ 51860  46385 ]                     |\n' );
  fprintf ( 1, '|                                             |\n' );
  fprintf ( 1, '|    B = [   2 ]                              |\n' );
  fprintf ( 1, '|        [ 305 ]                              |\n' );
  fprintf ( 1, '|                                             |\n' );
  fprintf ( 1, '|  Solve:                                     |\n' );
  fprintf ( 1, '|                                             |\n' );
  fprintf ( 1, '|    A * X = B                                |\n' );
  fprintf ( 1, '|                                             |\n' );
  fprintf ( 1, '|  Correct X:                                 |\n' );
  fprintf ( 1, '|                                             |\n' );
  fprintf ( 1, '|    - 46368.0                                |\n' );
  fprintf ( 1, '|      51841.0                                |\n' );
  fprintf ( 1, '|                                             |\n' );
  fprintf ( 1, '+=============================================+\n' );

  n = 2;

  a = [ 64079.0, 57314.0; ...
        51860.0, 46385.0 ];

  b = [ 2.0; 305.0 ];

  x1 = [ - 46368.0; 51841.0 ];

  x2 = a \ b;

  r = a * x2 - b;

  r8vec3_print ( n, x1, x2, r, '  Exact X, computed X, Residual' );

  return
end
