integrate [1/sqrt(u*(1-B*u))] du