integral of (sqrt(1 + x^2 -2x cos y)-x)/(2pi) for x=0..1,y=0..2pi