Weierstrass zeta function