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)