function f = inverse_fourier_transform(xi, F, x)

f = zeros(1, length(x));
ind = 0;
for m = 1:length(x)
    x0 = x(m);
    ind = ind + 1;
    integrand = (1/2/pi)* F.*exp(1j*xi*x0);
    f(ind) = trapz( xi, integrand );
end