Simplify definition of mid
Added find_roots function. Defaults to Newton; Krawczyk selectable by default.
Playing with Krawczyk to see if it can be made to work. I think the bisection is where the check for a root should be done, not the K operator