box% sage
┌────────────────────────────────────────────────────────────────────┐
│ SageMath version 10.0, Release Date: 2023-05-20                    │
│ Using Python 3.11.3. Type "help()" for help.                       │
└────────────────────────────────────────────────────────────────────┘
sage: R.<x,r> = PolynomialRing(QQ, order='lex')
sage: ideal(((x+r)^3 - 7)^2 - 50, (x^3 + 7)^2 - 50).groebner_basis()[-1]
r^21 - 140*r^18 + 3098*r^15 + 4788*r^12 - 702391*r^9 + 3996104*r^6 + 2190656*r^3 - 175616
sage: