theta = 0.25;
Esp = 0;
Esp2 = 0;
Var = 0;
for i=1:10000;
    u=rand(u);
    //disp(u);
    x=-log(abs(1-u))/theta;
    //disp(x);
    Esp = x/10000 + Esp;
    Esp2 = (x*x)/10000 + Esp2;
end;
disp(Esp);
Var=Esp2-Esp^2;
disp(Var);
