% Script File: no_pivot.m
% Examines solution to
%
%        [ delta 1 ; 1 1][x1;x2] = [1+delta;2]
%
% for a sequence of diminishing delta values.

disp(' Delta            x(1)                   x(2)  ' )
disp('-----------------------------------------------------')
for delta = logspace(-2,-18,9)
   A = [delta 1; 1 1];
   b = [1+delta; 2];
   x = genp(A,b);
   disp(sprintf(' %5.0e   %20.15f  %20.15f',delta,x(1),x(2)))
end