% function to calculate the Fourier Transform of the time series St function fourtran = ft(St,dt) fourtran = length(St)*fftshift(ifft(ifftshift(St)))*dt/sqrt(2*pi); return