function [ W, f ] = new_fft( x, Fs )
% [ W, f ] = new_fft( x, Fs )

if nargin<2
    Fs = 1;
end

W = fft( x );
N = length(x);
if mod(N,2)==0
    f0 = 0;
    fpositivos = (1:N/2)*(Fs/N);
    fnegativos = -( (N/2-1):-1:1 )*(Fs/N);
    f = [ f0; fpositivos(:); fnegativos(:) ];
else
    f0 = 0;
    fpositivos = (1:(N-1)/2)*(Fs/N);
    fnegativos = -( ((N-1)/2):-1:1 )*(Fs/N);
    f = [ f0; fpositivos(:); fnegativos(:) ];    
end
