diff --git a/src/klausurvorbereitung/lr_zerlegung.py b/src/klausurvorbereitung/lr_zerlegung.py new file mode 100644 index 0000000..9eac7b6 --- /dev/null +++ b/src/klausurvorbereitung/lr_zerlegung.py @@ -0,0 +1,37 @@ +""" +LR-Zerlegung + +Matrix A: + +[3, 2] +[1, 4] + +""" + +import numpy as np +import scipy as sc + +# Parameter +A = np.array([[3, 2], [1, 4]]) +pr = 3 + +[P, L, R] = sc.linalg.lu(A) + +with np.printoptions(precision=pr): + print(f"P = \n{P}\n\nL = \n{L}\n\nR = \n{R}") + +""" +Ausgabe: + +P = +[[1. 0.] + [0. 1.]] + +L = +[[1. 0. ] + [0.333 1. ]] + +R = +[[3. 2. ] + [0. 3.333]] +"""