>> syms theta a k;
>> z=a^2*k*sqrt(a^2+k^2)/2*int(theta*sin(2*theta),0,pi)
z =
-(pi*a^2*k*(a^2 + k^2)^(1/2))/4
>> syms t a;
>> x=a*(t-sin(t));
>> y=a*(1-cos(t));
>> s=int(sqrt(diff(x,t)^2+diff(y,t)^2),0,2*pi)
s =
8*(a^2)^(1/2)