function pi = odhad_pi(n)
s = 0;
for i = 1:n
    x = 2*rand-1;
    y = 2*rand-1;
    r = sqrt(x^2+y^2);
    if(r <= 1)
        s = s + 1;
    end
end

pi = 4*s/n;