import numpy as np
import scipy


if __name__ == '__main__':
    # Aufgabe 2c
    A = np.array([[0.8, 2.2, 3.6],
                   [2.0, 3.0, 4.0],
                   [1.2, 2.0, 5.8]])

    b = np.array([[2.4],
                   [1.0],
                   [4.0]])

    p, l, u = scipy.linalg.lu(A)
    print("\np:\n", p, "\nl:\n", l, "\nu:\n", u)