function x=logistrnd(a,b,n)
% logistic ranfom number generator
u=rand(1,n);
x=a+b*log(u./(1-u));