That hyperbolic-Bessel-function ratio was from

Continued Fraction -- from Wolfram MathWorld
One can use the

Fundamental recurrence formulas (Wikipedia) for them.

For a continued fraction

x = b0 + a1 /(b1 + a2/(b2 + ...))

we set up A's and B's such that

A(0) = b(0)

B(0) = 1

A(1) = b(1)*b(0) + a(1)

B(1) = b(1)

A(n) = b(n)*A(n-1) + a(n)*A(n-2)

B(n) = b(n)*B(n-1) + a(n)*B(n-2)

x(n) = A(n)/B(n) converges to the desired result as n -> infinity

So here we have a(n) = 1 and b(n) = i

^{n}
I evaluated with Mathematica's function RecurrenceTable, and I found no convergence for the first 100 terms. In fact, I find that the values recur with a period of 24. The recurrence period of b(n), 4, evenly divides this, so this recurrence is exact.

So one ought to do the general case of a(n) = 1 and b(n) = c*b[sup]n[/sub], and then take c -> 1 and b -> i.

Let's do A and B together.

p(n) = c*b

^{n}*p(n-1) + p(n-2)

p(0) = q0

p(1) = c*b*q0 + q1

We can use a generating function:

P(t) = sum

_{n from 0 to oo} p(n)*t

^{n}
That gives us

(1 - t

^{2})*P(t) - c*(b*t)*P(b*t) = q0 + q1*t

For b = c = 1, one gets the Fibonacci series:

P(t) = (q0 + q1*t)/(1 - t - t

^{2})

where each power of t is a sum of powers of ((sqrt(5)+1)/2) and ((sqrt(5)-1)/2).

If b is an nth root of unity, one can solve this equation by taking t -> b*t, t -> b

^{2}*t, until the cycle is complete. One can then solve for P(t) by factoring out P(b*t), P(b

^{2}*t), etc.

One gets a rational function of t, implying a Fibonacci-like solution.

For c = 1 and b = -1, one gets a power series in 12th roots of unity, yielding a recurrence period of 12, and for c = 1 and b = i, one gets 24th roots of unity, yielding a recurrence period of 24.

So one ought to treat c = 1 as a limit, to be approached from below.

For b = -1, the roots' absolute values are offset by about sqrt(1-c)/2 from unity, and for b = i, about sqrt(1-c)/4 from unity.

My patience has run out here; it's fairly straightforward algebra the rest of the way.