added a computation of the exact breaking threshold by solving a transcendental equation using Newton's method. The old breaking_thr routine was moved to pessimistic_breaking_thr()