thus converges(at least) quadratically [48, 5.4.2]. We use (a, b) as initial values for the iteration.