r = r - 1.6
            z = poly_tail_a[:num].inject(0) {|a, pn| a = a * r + pn}.fdiv(
                poly_tail_a[:den].inject(0) {|a, pn| a = a * r + pn}