function x=sys_solve(A, b)
    % utilizzando la fattorizzazione LU di A calcolo
    % la soluzione di L*U*x=b
    % L(U*x)=b
    
    [L, U]=my_lu(A);
    y=inf_solve(L, b) % Ly=b, y=Ux
    x=sup_solve(U, y);