def q_continued_fraction_complement(a, x, with_error = false)
          q = q_continued_fraction(a, x, with_error)
          with_error ? [1.0 - q.first, q.last + 2.0 * Float::EPSILON * (1.0 - q.first).abs] : 1.0 - q
        end