Мне нужна помощь с изменением значений внутри матрицы идентичности, вот код, который у меня есть
def gaussJordan(A,b):
n=len(b)
ide=eye(3)
Ab=concatenate((A,b), axis=1)
c=concatenate((Ab,ide), axis=1)
, и я хочу изменить ide
на другие значения
for e in range(n):
t=c[e][e]
for j in range(e,n+1):
c[e][j] = c[e][j]/t #Normalizar fila e
print('Normalizacion')
print(c)
for i in range(n):
if i!=e:
t=c[i][e]
for j in range(e,n+1):
c[i][j]=c[i][j]-t*c[e][j]
print('Reduccion')
print(c)
x=zeros([n,1])
for i in range(n):
x[i]=c[i][n]
return x
С этой матрицей значений
A=[[4.0,2.0,5.0],[2.0,5.0,8.0],[2.0,4.0,3.0]]
b=[[18.0],[27.3],[16.2]]
Функция
gaussJordan(A,b)