And hence we have a model of the volume of the capsule as a function of time t.

Input := 


vol[t_] = vsol[t]/.csolve[[1]]

Output =


            1/3

 -(k (36 Pi)    t)     1/3 3

(----------------- + v0   )

         3