dP/dt = (-a * P + G) (1 + e ^ -(t-b)/d)^(-1)