zad1.m

Satatistical program for variancy

clear all; clc; clf;

i=1;

for N=10:10:1000
    for K=1:1000
        y = rand(1,N);
        X(K) = 1/N*sum(y);
        
    end
    X2p(i) = 1/K * sum(X);
    Si(i)=1/(K-1)*sum((X-X2p(i)).^2);
    i=i+1;
end
plot(Si, 'b-');
hold on;
plot(Si, 'r*')