
function idrascript

A=[-0.370 0.050  0.050 0.070;
    0.050 -0.116 0     0.050;
    0.050 0     -0.116 0.050;
    0.070 0.050 0.050 -0.202];

b=[-2; 0; 0; 0];
[L,U,P]=lu(A);

% PA=LU. Risolvo LUx=PAx=Pb
bb=P*b;
y=triinf(L,bb);
xv=trisup(U,y);

% Soluzione esatta.
xv_esatta=A\b;

% Display.
A
b
n=size(A,1)
xv
xv_esatta