This function takes the n time domain samples (real or complex) and finds the psd by taking (fft/n)^2. The two sided spectrum is produced by shifting the psd The array freq provides the appropriate frequency values for plotting purposes.