alpha=1.375;
a=alpha-1;
[s (1)] = 0;
for i = 2:10000
m (i) = (-log(rand))^a;
s (i) = m (i) + s (i-1);
d (i) = s (i)/(i-1); 
endfor

plot (d);

