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 ) )