def MDL(N, rho, k):
    """Minimum Description Length

    .. math:: MDL(k) = N log \\rho_k + p \\log N

    :validation: results
    """
    from numpy import log
    mdl = N * log(rho) + k * log(N)
    return mdl
