function output = CalcFallingSpeed2(a,b,c,d)

output = sqrt(a*b/c)*tanh(sqrt(b*c/a)*d);