function y = sigmoid(u)


% function y = sigmoid(u,theta)
% y=1./(1+exp(-u/theta));
y=1./(1+exp(-u));
