integrate (u*(ln (m/(m-q*t))) - g * t) dt