integrate u^n sqrt(a + bu) du