where l = h^2 + (r1+r2)^2−−√