from numpy import array
from scipy.linalg import lu
a = array([[2.0, 4.0, 4.0, 4.0], [1.0, 2.0, 3.0, 3.0], [1.0, 2.0, 2.0, 2.0],
    [1.0, 4.0, 3.0, 4.0]])
pl, u = lu(a, permute_l=True)
Is there a standard solution for Gauss elimination in Python?