Plot the real value of the Fourier component of electric field?

I have attached a code below I am trying to plot the real value of the fourier component of electric field dq l, and I am taking contribution of fields from all positions dx j.

When I ListPlot3D for nn=200, the Mathematica quits its kernel by itself without any error or warning.
I wonder is it the memory issue or something else because the code works well for few points such as
nn=20. But I have to have nn=200.

Code-1.nb (10.2 KB)

Those are large symbolic expressions being generated. I’d suggest using approximate numbers in the definition of Eq so as to get a numerical table instead.

Also if you can redefine Eq in a way that uses Fourier rather than explicit sums, it would be much faster.

Thank you for the advise.