def DeadEndMetabolites(model, nodeadendmodel=None):
    """ iteratively identify peripheral metabolites """
    if nodeadendmodel == None:
        nodeadendmodel = model.NoDeadEndModel()
    return list(set(model.Metabolites()).difference(nodeadendmodel.
        Metabolites()))
