'Computes resultant of ``f`` and ``g`` via PRS.'

def resultant(f, g, includePRS=False):

    (lev, dom, per, F, G) = f.unify(g)
    if includePRS:
        (res, R) = dmp_resultant(F, G, lev, dom, includePRS=includePRS)
        return (per(res, kill=True), list(map(per, R)))
    return per(dmp_resultant(F, G, lev, dom), kill=True)
