function x = SolveAxEqualsb(A, b)
[U, c] = GaussianElimination(A, b);
x = BackSolve(U, c);

