function x=sys_solve(A,b);
%risolve un sis. lineare generico
%funziona solo quando lu(..) ritorna
%una fatt. lu senza pivoting
[L,U]=lu(A);
y=infsolve(L,b); %risolve L*y=b
x=supsolve(U,y); %risolve U*x=y