clear;clc;
% generacja sygnalu
% estymacja g阺to渃i

N = 100000;
przedz = 600;

NN = zeros(2*przedz,1);

x = randn(N,1);



a = min(x);
b = max(x);

W = (b-a)/przedz;

prog=(a:W:b);

delta = (max(abs(x)))/W;
delta = fix(delta/2);
delta = 300;

for i=1:N-1
    NN(fix(x(i)/W) + 1 + delta) = NN(fix((x(i)/W)) + 1 + delta) + 1;
end

for i=1:2*przedz-1
    p(i)=(1/W)*(NN(i)/przedz);
end

plot(p(1:(b-a)/W),'r+')