function v=logitpdf(x,mu,b)
v=1./(1+exp(-(x-mu)/b));
v=v.*(1-v)/b;
