Computational Physics, Physics, Linear Equation, Numerical Methods, Linear System.
k in range(q): sum = sum + L[i][k]*U[k][q] L[i][q] = (A[i][q] - sum)/U[q][q] return L, U Untuk memanggil fungsi ini dan mengimplementasikannya pada matriks A. Kita akan gunakan kode berikut. # matrix A A = [[3, 2, -1], [2, -2, 4], [-1, 0.5, -1]] # calculate L and U L, U = decomposition(A) # show L