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)