#### Finding the continued fraction of ep/q

The continued fraction expansions for e1/k and e2/k (k odd) are well known, namely:

However there is no known formula for the partial quotients of the continued fraction expansion of e3, or more generally ep/q, with p distinct from 1,2 and gcd(p,q)=1.
Here are the first 250,000 partial quotients of e3 and here is the table size-sorted. The largest partial quotient is a196440 = 546341.

We perform the algorithm of J.L. Davison, An algorithm for the continued fraction of ep/m, Proceedings of the Eighth Manitoba Conference on Numerical Mathematics and Computing (Univ. Manitoba, Winnipeg, 1978), 169-179, Congress. Numer., XXII, Utilitas Math.

The starting point is a result of R.F.C. Walters in Alternate derivation of some regular continued fractions, J. Austr. Math. Soc 8 (1968), 205-212): If

then pn/rn and qn/snel/m.

We first find the least n=n* such that pn, qn, rn, sn are non-negative and repeatedly apply Raney's factorisation for n*≤ k ≤ n*+N, as in Davison's example in §3.
The number (count) of partial quotients of ep/q found is returned.
We cannot predict the value of count, but it becomes positive for sufficiently large N.
We restrict taking p/q no greater than 14, because of the time taken to compute a0 for larger values of p/q.

This is a BCMATH translation of a BC program.
A faster C version is available.

Inside the program, p and q are replaced by p/gcd(p,q) and q/gcd(p,q).

Enter p (> 0):
Enter q (> 0):
Enter N (10000 ≥ N ≥ 0):