Implement decomposition based solve functions for solving Ax = b

I had written this wrong. `x_{i} = b_{i} / A_{i, i}`, so this is the correct formula since A is not a column vector.