l=101
x=rand(1:l);
g=1:1:length(x);
y=pi.*x-.5*pi.*g;
z=rand(1,l);
zz=-log(z);
plot(g,zz)