integrate sqrt((1+.5*cos(N*theta))^2+(-.5*sin(N*theta))^2) d(theta)