function out = hilbert(in,lin)
% HILBERT computes real hilbert transform of real input
% Has symbol -i sgn(k)
in = fft(in);
in(1:(lin/2+1)) = -1*in(1:(lin/2+1)); 
in = i*in;
out = real(ifft(in));