def jenkins_traub(coefficients, epsilon, max_iterations):
    """Find the roots of a polynomial using the Jenkins-Traub algorithm.
    Wrapper method for the Jenkins-Traub algorithm. Loops through the
    coefficients, finds a root, divides the coefficients by
    that smallest root to find a new set of coefficients, and then starts over.
    Returns a list of the roots.
    """
    # Remove leading coefficients equal to 0
    while len(coefficients) > 0 and coefficients[0] == 0:
        coefficients.pop(0)