d/dt(A*(1+g*t)e^(-q*t)+1/2*a*t^2*e^(-q*t))