Solution method for tri-diagonal system