(* #53, Section 3.9 *)
a=2t+1
v=Integrate[a,t] +c
s=Integrate[v,t]+d