pi.integrate(x^1/2+1)^2 from 0 to 2