function y=rfft(x)
n=length(x);
if floor(log2(n))~=log2(n)
    error('nnon potenza di 2');
end

if n==1
    y=x;
else
    xodd=x(1:2:n-1);
    xeven=x(2:2:n);
    todd=rfft(xodd);
    teven=rfft(xeven);
    
    w=exp((2*pi*i/n)*(0:n/2-1))'; %vettore che contiene le radici dell'unita'
    
    y0=todd+w.*teven;
    y1=todd-w.*teven;
    y=[y0;y1];
end
