function Pi = Pi1(N)
%
s = 0;
for i = 1:N
    x = rand(1,2);
    if (x(1)^2+x(2)^2 <=1)
        s = s+1;
    end
end
Pi = 4*s/N;
