from sage.all import *
p = 1125899906842679
R = PolynomialRing(GF(p), 'x,y')
R.inject_variables()
I = Ideal(3*x**4 + 740736*x**2 + 4854540*x + 1125884665458743, -x**3 + y**2 + 1125899906719223*x + 1125899906438134)
J = I.elimination_ideal(y)
Here, elimination_ideal leads to an error. I have no clue why. For other examples, like other finite fields, it works. Has someone an idea, why?