disp(' '), disp(' question 1) b)')
disp(' ')
A=toeplitz([2 -1 zeros(1,2)])+diag([2 3 4 5])
n=4; a=-ones(3,1), b=2+[2 3 4 5]', c=a
X0=[0 1 10 1000]';
Y=A*X0

disp(' '), disp(' La solution exacte de l''équation AX=Y est')
X=X0

disp(' ')
disp(' Les erreurs relatives obtenues en résolvant cette équation AX=Y')
disp('directement par Matlab ou à l''aide de la fonction solLU valent')
fprintf('respectivement : %.3g et %.3g.\n\n'...
,norm(A\Y-X0,'inf'),norm(solLU(n,a,b,c,Y)-X0,'inf'))