Operaciones con matrices

Para poder realizar operaciones que estén bien hechas con matrices, se debe de importar de numpy la librería matrix.

Ya solo antes de un arreglo solo pones matrix([elementos]), y así se pueden realizar hasta divisiones de matrices. Tambien se pueden sacar determinantes sin necesidad de un ciclo for.

Solo se pone np.linalg.det(matriz)  y ya te calcula el determinante de cualquier matriz cuadrada.
Ya si se quiere sacar la dimensión y el tamaño se usa .ndim y .shape respectivamente, después del nombre de la variable donde se tenga el arreglo, y también se debe de importar de numpy el np.

• Código:

#operaciones con matrices

import numpy as np

from numpy import matrix

a=np.array([[1,2], [3,4]])

print(np.linalg.det(a))

a=np.array([[[1,2],[3,4]],[[1,2],[2,1]],[[1,3],[3,1]]])

print(a.shape)

print(a.ndim)

print(np.linalg.det(a))

b=np.array([[3,2,1],[6,2,-5],[-2,1,2]])

print(b.shape)

print(b.ndim)

print(np.linalg.det(b))

#ejercicio

x=matrix([[2,4,1],[1,-2,3],[5,0,-1]])

y=matrix([[3,-1,-2],[0,5,6],[0,0,9]])

z=matrix([[2,0,-1],[0,-1,2],[1,-2,5]])

c=np.array([[3,2,1],[0,2,-5],[-2,1,4]])

p1=-x-y+z

p2=x+y-z

p3=(3*x) + (z/2)

print("\nParte 1")

print(p1)

print("\nParte 2")

print(p2)

print("\nParte 3")

print(p3)

print("\nDeterminate de matriz: %0.2f"%np.linalg.det(c))

• Ejecución: