factorización LU
Publicado por Gerardo (7 intervenciones) el 12/06/2019 06:55:30
Buenas, tengo este algoritmo,debe recibir como parámetro una matriz A y retornar dos matrices L Y U que serán su factorización LU, pero no sé como acomodarlo para que compile, gracias.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
import numpy as np
def LU(A):
(m,n) = A.shape
(m,n) = L.shape
(m,n) = U.shape
for k in range (1,n):
L[k,k] = 1
for i in range (k+1,m+1) :
L[i,k] = A[i,k]/ A[k,k]
U[i,k] = 0
for j in range (k+1, n+1):
U[i,j] = A[i,j] - L[i,k]*A[i-1,j]
return L
return U
A = [[10,2,4],[2,8,1],[4,1,20]]
K = LU(A)
print (K)
Valora esta pregunta
0