Я в своё время делал так:

;

.
Т.о. задача свелась к вычислению арктангенса.
Далее, пользуясь тем, что арктангенс --- нечётная функция и формулой

,
выражаем арктангенс любого аргумента через арктангенс числа от 0 до 1.
Затем, пользуясь формулой (возможно, несколько раз)

,
сводим аргумент к отрезку
![$[0,1/16]$ $[0,1/16]$](https://dxdy-03.korotkov.co.uk/f/2/a/9/2a90328517d3f057910ab0fe073bd82782.png)
.
После этого применяем разложение в ряд Тейлора:
с оценкой погрешности (для
![$x\in[0,1/16]$ $x\in[0,1/16]$](https://dxdy-01.korotkov.co.uk/f/0/2/4/024ca07ce6a9af89f8107a06a30fa30b82.png)
):
(выбираем N таким, чтобы эта оценка нас удовлетворила).
Не гарантирую, что это самый быстрый из известных методов вычисления, но работает неплохо.