function F = fourier_transform(x, f, xi)

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