FIXED is_integral. Very slow, but at least it seems to be accurate. Using sympy to calculate all roots in the Z domain. Also unit_tests.py uses pyparsing to handle input.