import linalg as la
print(
la.Solution(
la.System(
la.TridiagonalMatrix(
[4, 9],
[23, 5],
[5, 16, 2]
),
la.Vector(
[5, 8, 1]
)
),
la.ThomasAlgorithm() # optional
)
)