d=t(r+a*e^(-0.013*t))