clear all;
clc;

A = magic(10);

x_exact = [1:10]';

b = A*x_exact;

x = lu_solver(A,b);

err = norm(x-x_exact);

fprintf('\n \t Errore con la soluzione esatta: %1.3e \n',err)