integrate from 0 to x P/2* (a*cos(theta)+h*sin(theta)) d theta