function y = derivfreefall(cd, v, g, m, t)
	y = sqrt(g*m./cd).*(tanh(sqrt(g*cd./m)*t) - ...
	    t*sqrt(g.*cd/m).*sech(sqrt(g.*cd/m)*t).^2)./(2*cd);
return