function x=inversefouriertransform(X,omega);
%This function computes an the inverse Fourier transform, generated by the
%function fouriertransform.

step=1/(length(omega)/2/pi*(omega(2)-omega(1)));
x=ifft(fftshift(X))/step;