function [L,U,x]=Solver(A,b)
[L,U]=LU(A);
x=LT_Solver(L,b);
x=UT_Solver(U,x);