function []=prob_9_1_3()
%
% Written by:
% -- 
% John L. Weatherwax                2006-12-14
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----
  
clc; close all; 

disp( 'the (3x3) Hilbert matrix' ); 
A = hilb(3)

x_0 = [ 1 1 1 ].'; 
A*x_0 

x_1 = [ 0, 6, -3.6 ].'; 
A*x_1

disp( 'the difference between the two' ); 
A*x_0 - A*x_1