integrate sqrt(U^2 - 2Y*g*t + (g*t)^2) dt