from NumPy import *
def Mult(str1,str2): # строки в виде 0;2;6
   arr1=str1.split(';')
   arr2=str2.split(';')
   if len(arr1)==len(arr2):
      for x in arr:
         res=res+arr1[x]*arr2[x]
      return res
   else:
   return "Не могу перемножить."
a=[[input() for y in range(4)] for x in range(4)]
b=[[input() for y in range(4)] for x in range(4)]
res=[[input() for y in range(4)] for x in range(4)]
for y in range(4):
   for x in range(4):
      res[y][x]=Mult(';'.join(a[y],';'.join(b[:,x])

>>> import numpy as np
>>> a = np.array(range(16))
>>> a = a.reshape(4,4)
>>> a
array([[ 0,  1,  2,  3],
       [ 4,  5,  6,  7],
       [ 8,  9, 10, 11],
       [12, 13, 14, 15]])
>>> b=np.array(range(15,-1,-1))
>>> b = b.reshape(4,4)
>>> b
array([[15, 14, 13, 12],
       [11, 10,  9,  8],
       [ 7,  6,  5,  4],
       [ 3,  2,  1,  0]])
>>> a*b
array([[ 0, 14, 26, 36],
       [44, 50, 54, 56],
       [56, 54, 50, 44],
       [36, 26, 14,  0]])
>>> np.dot(a, b)
array([[ 34,  28,  22,  16],
       [178, 156, 134, 112],
       [322, 284, 246, 208],
       [466, 412, 358, 304]])

def Mult(arr1,arr2):
   if len(arr1)==len(arr2):
      for x in arr:
         res=res+arr1[x]*arr2[x]
      return res
   else:
   return "Не могу перемножить."
a=[[input() for y in range(4)] for x in range(4)]
b=[[input() for y in range(4)] for x in range(4)]
res=[[input() for y in range(4)] for x in range(4)]
for y in range(4):
   for x in range(4):
      res[y][x]=Mult(a[y],b[:,x])