sigma = 1/(1 + np.exp(-a*(t-k-b)))