Computational Physics, Physics, Linear Equation, Numerical Methods, Linear System.
dan U(k,q). # lower matrix q = p for i in range (q,n): if (i==q): # diagonal of L L[i][q]=1 else: # summation of L(i,k)*U(k,q) sum = 0 for k in range(q): sum = sum + L[i][k]*U[k][q] L[i][q] = (A[i][q] - sum)/U[q][q] Dengan menambahkan implementasi rumus umum matriks U dan L tersebut ke kode