link approx
link genrfncs
link periodic
link rational

procedure sig2cf(a, b, c, limit)
   local n, m, j, rat, seq, results, line

   /limit := 5 * (log(b) * sqrt(b))

   every rat := sqrtapprox(b) do {
      seq := []
      rat := divrat(addrat(rational(a, 1, 1), rat),
         rational(c, 1, 1))
      every j := cfseq(rat.numer, rat.denom) \ limit do
         put(seq, j)
      results := repeater(seq, 1.5)
      if *results[2] ~= 0 then break
      }

   return results

end