(integrate exp(-(x-mu)^2/(2sigma^2)+I*n*pi*x/L) for x=-L to L)