monte carlo pi

format long;
len = 6;
rzuty = [5 10 100 1000 10000 100000];
bledy = zeros(1,len);
for i=1:len
    trafienia = 0;
    for j=1:rzuty(i)
        pos = rand(1,2);
        if pos(1)^2 + pos(2)^2 <= 1
            trafienia = trafienia + 1;
        end
    end
    pole = trafienia / rzuty(i) * 4;
    bledy(i) = abs(pole - pi);
end
srednia = mean(bledy)
srednia = ones(1, len) * srednia
semilogx(rzuty,bledy,rzuty,srednia)