%  anonymous functions
%
freefall = @(cd, v, g, m, t) [v - sqrt(g*m./cd).*tanh(sqrt(g.*cd/m)*t)]
derivfreefall = @(cd, v, g, m, t) [sqrt(g*m./cd).*(tanh(sqrt(g.*cd/m)*t) - t*sqrt(g.*cd/m).*sech(sqrt(g.*cd/m)*t).^2)./(2*cd)]

%  inline functions
%
freefall = inline('v - sqrt(g*m./cd).*tanh(sqrt(g.*cd/m)*t)', 'cd', 'v', 'g', 'm', 't')
derivfreefall = inline('sqrt(g*m./cd).*(tanh(sqrt(g.*cd/m)*t) - t*sqrt(g.*cd/m).*sech(sqrt(g.*cd/m)*t).^2)./(2*cd)', 'cd', 'v', 'g', 'm', 't')