it gives wrong numbers in calculating Weierstrass zeta function

[Zou-Bo]

I am doing the calculation that follows the equation (13) and (14) in https://doi.org/10.1063/1.5042618
I use tau = 1+sqrt(3)*im. It turns out that it gives incorrect results when tau varies slightly.