d a*(1-exp(-0.25*t)) /dt