% script PivotMatTest
clc
for delta = logspace(-2,-18,9)
    A=mat1(delta);
    [L,U,P]=lu(A);   
    b=[1+delta;2];
    y=LTriSol(L,P*b);
    x=UTriSol(U,y);
    disp(sprintf(' %5.0e  %20.15f  %20.15f',delta,x(1),x(2)))
end
