function y=fx(x) y=5/sqrt(2*9.81*x)-tanh(3/(2*5)*sqrt(2*9.81*x)); endfunction root=fzero(@fx,3)