function R=hilbert(x);
// x is vector!
x=real(x);
h=zeros(x);
n=length(x);
R=fft(x,-1);
h(1)=1;
if modulo(n,2)==0;
h(n/2+1)=1;
h(2:n/2)=2;
else
h(2:(n+1)/2)=2;
end;
R=fft(R.*h,1);
endfunction;




