dP/dt = (-k1*P+G)*(1/(1+exp(-(t-k2)/d)))