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