newton-solver now also parametrized with polynomial