trabajoMatrix

SolucionMetodoCramerMatrix15x3.py

import pandas as pd
import numpy as np
from numpy.linalg import inv
from scipy import linalg

df = pd.read_table('DataEjercicio7', delimiter='\t',decimal=',')
matrix1 = df.loc[:, ['2', '3', '4']]
#matrix1.astype(float)
#matrix1 = matrix1.values

vectort7 = df.loc[:, '11']
#vectort7 = vectort7.values

matrix1T = matrix1.T

#La multiplicación de dos matrices, Matriz transportada (Mt) * M MMULT
matrix1Mult = matrix1T.dot(matrix1)

#Buscar el q de la mariz q = (Q/n)
matrix1D = matrix1Mult.div(len(matrix1.index.values))

#la inversa del la matriz (M^-1) MINVERSA
#matriz1Minversa = inv(matrix1D)
matriz1Minversa = linalg.inv(matrix1D)

#el deteterminante total de la matriz
Delta = np.linalg.det(matriz1Minversa)

#la multipicación de una matriz transpuesta por lun vector
matrix1vectort7 = matrix1T.dot(vectort7)
print(matrix1vectort7)

#Buscar el q de la mariz q = (Q/n)
matrix1Dvetor7 = matrix1vectort7.div(len(matrix1.index.values))
matrix1Dvetor7 = matrix1Dvetor7
matrix1Dvetor7.T.iloc[0:]

for i in range(len(matrix1vectort7)):
    matriz1Minversa1 = matriz1Minversa.copy()
    #print(matriz1Minversa1)
    matriz1Minversa1[:,i]= matrix1Dvetor7.T.iloc[0:]
    print(matriz1Minversa1)
    Di = np.linalg.det(matriz1Minversa1)
    print(Di)
    out = Di/Delta
    print(out)