clear
N = 100; % numero de pontos a serem gerados
Nd = 0;
for i=1:N
    x = rand;
    y = rand;
    if (x^2 + y^2 )<=1
        Nd = Nd + 1;
    end
end
pi_exp = 4 * Nd / N