Solve[ z == L - (Sqrt[ 1 + (b z)^2 ] - 1) / b, z ]