xxxxxxxxxx
def multiply_matrix(a, b):
if len(a[0]) != len(b):
return 'matrix multiplication not possible'
f = []
for i in range(len(a)):
f.append([])
for j in range(len(b[0])):
f[i].append(0)
for index in range(len(b)):
f[i][j] += a[i][index] * b[index][j]
return f
TEST:
a = [
[23, 41, 56],
[85, 26, 63],
]
b = [
[17, 34],
[18, 38],
[97, 69],
]
print(multiply_matrix(a, b))
OUTPUT:
[[6561, 6204]
[8024, 8225]]
xxxxxxxxxx
import numpy as np
matrix1 = np.array([[1, 2], [3, 4]])
matrix2 = np.array([[5, 6], [7, 8]])
result = np.dot(matrix1, matrix2)
print(result)
xxxxxxxxxx
A = np.array( [ [1,1,0],[1, -2, 1], [2,1,2]]) #C^3*3
X = np.array( [[1,2,-2]]) #C^1*3
C=A@X