def matvec_p(v, A, a, mu, p):
    for i in range(p):
        v = (a * (A.dot(v) - (mu * v)))
    return v
def matvec_p(v, A, a, mu, p):
    for i in range(p):
        v = (a * (A.dot(v) - (mu * v)))
    return v