function y = freefall(cd, v, g, m, t)
	y = v - sqrt(g*m./cd).*tanh(sqrt(g.*cd/m)*t);
end