Fourier.m

fourier_function = @(x, k) (2./(pi .^2)) * (1./(k.^2)) * (((-1).^k)-1);
test_range = 0:0.1:10;
k_range = 1:10;

fourier_values = zeros(size(test_range));

for k = k_range
    values = fourier_function(test_range, k) .* e.^(i * k * test_range);
    fourier_values = fourier_values + values;
end

plot(test_range, fourier_values);
input("")