function y=sigmoid(x)
y=1./(1+exp(-x));
