v0*t-v1*integral ln(1-q/m*t)dt-1/2*a*t^2