Most symbolic algebra systems provide a function, typically called "prem" or something like that, to compute the pseudo-remainder of two polynomials. Sage does provide other elementary functions for manipulating polynomials, e.g. "factor", but does not seem to provide this one. Perhaps something could be done by using the sage.rings.polynomial package, but that would be overkill for my immediate problem. Is there ANY way to do this simply?