restart;
L:=int(sqrt(1+cos(t)^2),t=0..Pi);
evalf(L);
f:=t;
g:=(1/3)*sqrt(t)*(3-t);
plot([f,g,t=0..3]);
df:=diff(f,t);
dg:=diff(g,t);
expand(dg);