for n = 1:1e5
   y = real(ifft(fft(y)));
endfor
n_iterations = n
figure(1)
plot(x,y_original,'b', x,y,'r');
figure(2)
plot(x,y - y_original)
