used another method to solve Ax = b, much faster than BiCG for 1D