function y=sigmm(x)
   y=1./(1+efun(-x));



function y= efun(x)
   y= exp(x)
