integrate 1/(sqrt(1+x^2) sqrt((a-b)+(a+b)*x^2)) assuming x>0