import numberdb.sage as numberdb
from sage.arith.misc import bernoulli
from sage.rings.complex_arb import ComplexBallField
from sage.rings.integer_ring import ZZ
from sage.rings.rational_field import QQ
CBF = ComplexBallField(numberdb.bits(100, losing=96))
k = ZZ(1)
H = sum(QQ(1) / QQ(j) for j in range(1, int(k) + 1))
log_A = CBF(QQ(bernoulli(k + 1)) * H / QQ(k + 1)) - CBF(-k).zetaderiv(1)
log_A.exp(), log_ABoth values are computed as Sage complex balls with arb at numberdb.bits(digits, losing=96) bits. The imaginary part is required to contain zero before the real part is returned, and $A$ is checked against (2).