Hola, necesito ayuda para descubrir mi error el el siguiente programa python, que pretende calcular el producto de matrices, suponiendo de antemano que se cumple el requisito de dimension para poderse multiplicar.

def producto_matrices(a,b):
    nfilas1=len(a)
    ncols1=len(a[0])
    nfilas2=len(b)
    ncols2=len(b[0])
    filaprod=[0]*nfilas1
    for i in range(nfilas1):
        for j in range(ncols2):
            prod=[]
            for k in range(nfilas2):
                filaprod[i]+=a[i][k]*b[k][j]
                filaprod.append(filaprod[i])
    return prod.append(filaprod)

La matriz donde vas a guardar el producto, está llena de ceros? debes primero llenar esa matriz de ceros para poder ir agregando los productos respectivos. Puedes utilizar la librería Numpy o lo puedes hacer por iteraciones. Ten en cuenta que la dimensión de la matriz del producto debe ser: si el producto es (M*N) *  (N*I) entonces las dimensiones serán: (M*I)