TODO: use
org.apache.commons.math3.analysis.solvers.NewtonRaphsonSolver
to solve polynomial functions