m = 68.1; g = 9.81; c = 0.25; t = (0:0.1:20)'; v = sqrt(m*g/c)*tanh(sqrt(g*c/m)*t);