%% Para gerar Pi a partir de 10.000 pontos
clear
N=1e4;

4*mean( sum( rand(2,N).^2 )<1 )

%%
% analise dos resultados

PIs = [ 3.1444  3.1532 3.1608 3.1264 3.1208 3.1248 3.1280 ...
    3.1524 3.1460 3.1356 3.1604 3.1616 3.1344 3.1512 3.1416 ...
    3.1456 3.1444 3.1224 3.1348 3.1628 ]

%%
clear
N=1e4;
PIexp = 4*mean( sum( rand(2,N,100).^2 )<1 )
