Eu tenho essa função aleatória para calcular piEstilo Monte Carlo:

max=10000000;
format long;

in = 0;
tic
for k=1:max
    x = rand();
    y = rand();
    if sqrt(x^2 + y^2) < 1
        in = in + 1;
    end
end

toc
calc_pi = 4*in/max 
epsilon = abs(pi - calc_pi)

calcPi(100000000);

Se eu pudesse repetir isso 10 a 100 vezes, esse algoritmo poderia competir como recorde mundial? Em caso afirmativo, como posso encontrar o número de iterações que fornecerá o dígito Nth?