bulirsch-stoer algorithm