“…In this paper, we continue the study of efficient algorithms for the computation 8 of the Riemann zeta function over the complex plane, introduced by Borwein [7] 9 and extended by Belovas [1], Belovas and Sabaliauskas [4], Belovas, Sakalauskas and Throughout this paper, we denote by Φ(x) the cumulative distribution function of the standard normal distribution, and by Φ(x) we denote the corresponding tail distribution Φ(x) = 1 − Φ(x). Γ(s), B(x, y) and W(x) denote the gamma function, the beta function and the Lambert W function respectively.…”