import numpy as np
from scipy import linalg as la
A=np.array([[5,6,-4],
[2,1,-2],
[7,9,2]], dtype=np.float64)
b = np.array([10,14,-11], dtype=np.float64)
x = la.solve(A,b)
print(x)
print(A @ x)
L,U=la.lu(A,permute_l=True)
d=la.solve(L,b)
x=la.solve(U,d)
print(x)
print(A @ x)
aW1wb3J0IG51bXB5IGFzIG5wCmZyb20gc2NpcHkgaW1wb3J0IGxpbmFsZyBhcyBsYQoKQT1ucC5hcnJheShbWzUsNiwtNF0sCiAgICAgICAgICAgIFsyLDEsLTJdLAogICAgICAgICAgICBbNyw5LDJdXSwgZHR5cGU9bnAuZmxvYXQ2NCkKYiA9IG5wLmFycmF5KFsxMCwxNCwtMTFdLCBkdHlwZT1ucC5mbG9hdDY0KQp4ID0gbGEuc29sdmUoQSxiKQpwcmludCh4KQpwcmludChBIEAgeCkKTCxVPWxhLmx1KEEscGVybXV0ZV9sPVRydWUpCmQ9bGEuc29sdmUoTCxiKQp4PWxhLnNvbHZlKFUsZCkKcHJpbnQoeCkKcHJpbnQoQSBAIHgp