Substitute nx = y, dx=dy/n. Now we have this integral: . Since in the interval a-1<y<a for some positive integer a , we can consider the equivalent integrals

where erf is the error function. As n goes to infinity this converges to some constant around 0.683938, but I don't know how to tackle that sum analytically. Hopefully this can help you out though =).