function [prob]=sigmoid(x,sigma,lapse)
     
    prob=lapse+((1-2*lapse)./(1+exp(-x./sigma)));
    
end