def matrixMultiply(A, B):
  # 獲取A的行數和列數
  A_row, A_col = shape(A)
  # 獲取B的行數和列數
  B_row, B_col = shape(B)

  # 不能運算情況的判斷
  if(A_col != B_row):
    raise ValueError

  # 最終的矩陣
  result = []

  # zip 解包后是轉置后的元組，強轉成list, 存入result中
  BT = [list(row) for row in zip(*B)] 

  # 開始做乘積運算 
  for A_index in range(A_row):
    # 用于記錄新矩陣的每行元素
    rowItem = []
    for B_index in range(len(BT)): 
      # num 用于累加
      num = 0   
      for Br in range(len(BT[B_index])): 
        num += A[A_index][Br] * BT[B_index][Br]
      # 累加完成后，將數據存入新矩陣的行中
      rowItem.append(num) 
    result.append(rowItem) 
  return result