To write a program to find the LU Decomposition of a matrix.
- Hardware โ PCs
- Anaconda โ Python 3.7 Installation / Moodle-Code Runner
- Define the package as scipy.linalg import lu.
- Get input from user and print L and U matrix by 'print'.
- Define a package as "from scipy.linalg import lu_factor, lu_solve" and create the variable as 'X' include the package in that variable
- print the variable 'X'
(i) To find the L and U matrix
#Developed by: Kannan.S
#RegisterNumber: 212223230098
import numpy as np
from scipy.linalg import lu
A=np.array(eval(input()))
P,L,U=lu(A)
print(L)
print(U)
(ii) To find the LU Decomposition of a matrix
#Developed by: Kannan.S
#RegisterNumber: 212223230098
# To print X matrix (solution to the equations)
import numpy as np
from scipy.linalg import lu_factor, lu_solve
A = np.array(eval(input()))
b = np.array(eval(input()))
lu, piv = lu_factor(A)
X = lu_solve((lu , piv),b)
print(X)
(i) To find the L and U matrix
(ii) To find the LU Decomposition of a matrix
Thus the program to find the LU Decomposition of a matrix is written and verified using python programming.