t=var('t')
a=var(a')
r=a*(1-cos(t))
b=r*r
c=(integrate(b,(t,0,2*pi))/2)