R.<a,b,c,d,x> = PolynomialRing(ZZ)
y = a*x^3 + b*x^2 + c*x + d

x0 = -b/(3*a)
y0 = y(x=x0)

assert y0-y(x=x0-x+x0)+y0 == y