\frac{\mathrm{dPain}}{\mathrm{d}t} = -\left(k_1 \text{Pain} + \ \ \right) \left(\frac{1}{1+e^{-(t-k_2)/d}}\right)