Use groebner basis reduction.
Fix docstring for `join_systems` function.