clear();
theta = 0.5;
N=10000;
Esp=0;
Var=0;
alpha=0;
Variance=0;
x=0;
for i=1:N
x=1-rand();
alpha= (log(x))/(-1 * theta) ;
Esp = Esp + (alpha/N);
Var = Var + (alpha * alpha / N);
end;
Variance = Var - (Esp*Esp);
disp(Esp);
disp(Variance);
