((m^2 + 1) x^2 + 2 m (sqrt(r^2 - u^2) - m u) x + m^2 u^2 - 2 m u sqrt(r^2 - u^2) - u^2)/(m^2+1)