k = 1:2:20;

bk = 20/pi * (1./k - cos(k/1000) ./ k + 2*cos(k/2000) ./ k);

x = 1:20000;
y(x) = 0;
y(k*1000) = bk((k-1)/2+1);
plot(x, y, 'LineWidth', 1.5)
title('Négyszög jel spektruma', 'FontSize', 15)
xlabel('frekvencia [ Hz ]',  'FontSize', 12)
ylabel('amplitudó [ V ]',  'FontSize', 12)
set(gca,'XTick',[1000:2000:19000])