clear; clc;

N = 1e5;
m = 0;
for i = 1 : N
    
    s = 0;
    while s < 1
        s = s + rand(1);
        m = m + 1;
    end
    
end
e_estimator = m / N