some changes in modgcd() (bivariate and multivariate case)

Some bivariate helper functions are replaced by multivariate ones.
add some (very short) docstrings and fix a minor error