integral u^2*(a+bu)^(1/2) du