function  c=rfft(x)

% c=rfft(x)
% fft via recurrence; calls on rft

n=length(x);
w=0:n-1;
w=exp(-(2/n)*pi*i*w);
c=rft(n,x,w);
