load "a.txt" 
load "b.txt" 
load "x.txt" 
fdelta = max(abs(A * x - b)) ;
disp("Solution delta:") ; 
disp(fdelta) ;
