Add "tolerance" and "debug" as keyword arguments to `newton`

Added "tolerance" and "debug" as keyword arguments to newton,
following the implementation of David (branch improve_root_finding).
They have as default values `eps(one(T))` (with T the type of
the elements of the interval) and `false`.