integrate (x*y-x^2)^(-1/2) from 0 to y