使用for循环的Python矩阵乘法

from numpy import linalg
from numpy import matmul



M = [[2,1],[1,3]]
b = [[3],[-1]
Mb = [] 

nrows = 2
ncols = 2
for i in range(0,nrows):
    sum = 0
    for j in range(0,ncols):
        sum = sum + M[i,j] * b[j]
     Mb[?] = ?   

print(Mb)

#其目的是使用循环乘法矩阵。当我试图运行这个程序时，我一直会遇到错误（python列表索引必须是整数而不是元组）。我不确定我需要用什么来代替插入的问号。