%estymacja 渞ednio kwadratowa
clear all; clc; clf;

i=1;

for N=10:10000
    y = randn(1,N);
    if (i == 1)
        Psi2(i) = 1/N*sum(y.^2);
    else
        Psi2(i) = 1/N*(  (N-1)*Psi2(i-1)  + y(N).^2);
    end
    i=i+1;
end

XX(1:N)=1;
plot(Psi2, 'b*');
hold on;
plot(XX, 'g-')