Thank you indeed, CRGreathouse, for signalling this error.
Using Merten's Theorem, I find, that g_alt must be divided by 2*exp(-gamma)=1,122918967...,
with gamma=0,577215664901... being Euler-Mascheroni's constant.
Then, g_alt will be asymptotically equivalent to g_HL
I'm now working...