%%
% Para avaliar o tempo para geração de Pi por Monte Carlo
% usando loops
disp( 'Usando Loops' )
clear
tic
N=1e5;
M=400;
PIexp = zeros(M,1);
for qPI = 1 : M
    n = 0;
    for i = 1 : N
        x = rand;
        y = rand;
        if (x^2 + y^2)<1
            n = n + 1;
        end
    end
    PIexp(qPI) = 4*n/N;
end
toc
% std( PIexp )


%%
% Para avaliar o tempo para geração de Pi por Monte Carlo
% usando um codigo matLab
disp( 'Usando um codigo MatLab' )
clear
tic
N=1e5;
M=400;
PIexp = 4*mean( sum( rand(2,N,M).^2 )<1 );
toc
% std( PIexp )
