(*Zmodyfikowany wariant poprzedniego programu. Nie ma już kryterium testowania "trafień". Każdy punkt jest włączany, ale z losową wagą.*)

L=10000;

V=1;

n=0;
err=0;

Do[
x=RandomReal[];
n=n+Exp[-x];
err=err+Exp[-2x];
,{i,1,L}]

n

(*Błąd jest mniejszy, niż w poprzednim przypadku.*)

V*N[Sqrt[(err/L-(n/L)^2)/L]]

V*N[n/L]

NIntegrate[Exp[-x],{x,0,1}]