function y=ifft_real(beta,x)
N=length(beta);
y=beta(1); 
for k=1:N/2
   y=y+(beta(k+1)+beta(N-k+1))*cos(k*x)+i*(beta(k+1)-beta(N-k+1))*sin(k*x);
end
y=2*y/N;
