function errors_test06 ( )

%*****************************************************************************80
%
%% ERRORS_TEST06
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    22 August 2017
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, '+=============================================+\n' );
  fprintf ( 1, '|                                             |\n' );
  fprintf ( 1, '|  Exercise 6                                 |\n' );
  fprintf ( 1, '|                                             |\n' );
  fprintf ( 1, '|    A = [ 64919121    - 159018721 ]          |\n' );
  fprintf ( 1, '|        [ 41869520.5  - 102558961 ]          |\n' );
  fprintf ( 1, '|                                             |\n' );
  fprintf ( 1, '|    B = [ 1 ]                                |\n' );
  fprintf ( 1, '|        [ 0 ]                                |\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, '|    205,117,922.0                            |\n' );
  fprintf ( 1, '|     83,739,041.0                            |\n' );
  fprintf ( 1, '|                                             |\n' );
  fprintf ( 1, '+=============================================+\n' );
  fprintf ( 1, '\n' );

  n = 2;

  a = [ 64919121.0, -159018721.0; ...
        41869520.5, -102558961.0 ];

  b = [ 1.0; 0.0 ];

  x1 = [ 205117922.0; 83739041.0  ];

  x2 = a \ b;

  r = a * x2 - b;

  r8vec3_print ( n, x1, x2, r, '  Exact x, computed x, residual:' );

  return
end
