sqrt(u/r2)(1-sqrt(2r1/(r1+r2)))