fN = 2^8
k = 0:fN;
ps = 1./(1+k.^2);
signal = randsig(ps);
plot(signal)