function v = velff(vt, g, t0, t)
%velff Summary of this function goes here
%   Detailed explanation goes here
v = vt*tanh(-g/vt*(t-t0));
end

