% LMex030514.m
clear all
y=rand(1000,1);
stick=[y./(1-y)];
i=find(stick>1);
stick(i)=1/stick(i);
mean(stick)
median(stick)
% Not at all close to book
syms y
int((1-y)/y,0.5,1)
% matches book
hist(stick);figure(gcf)

