int_0^1 ln((x+1)/(2x^2))/sqrt(x^2+2x) dx, with 20 digs