There is an infinite series, I want to transform it into a function, with digamma functions or something else. I hope someone can provide some guidance and suggestions. $$ E(x,y)=\sum_{n=-\infty}^{\infty} (\frac{1}{\sqrt{(2n-2x)^2+y^2}}-\frac{1}{\sqrt{(2n)^2+y^2}}), \qquad 0<x<1 $$
I have plotted this function with Matplotlib by recursion, shown below (with y = 0.025).