Impact Probability Computation for Neo Resonant Returns Through a Polynomial Representation of the Line of Variations