Solve ndarray Ax=b via LU factorization with pivoting - View it on GitHub