solve returning infeasible solution

[lukamar]

Hi,
By solving the system in giac 1.7.0-25

Code : Tout sélectionner

sol:=solve([(2*x1^5+4*x1^3+2*x1*x2+2*x1)/(x1^2*x2+x2),x1^4+2*x1^2+x2^2-3],[x1,x2])

I get four solutions:

Code : Tout sélectionner

[[0,sqrt(3)],[0,-(sqrt(3))],[sqrt((sqrt(2*(sqrt(17)-1))-2)/2),(-(sqrt(17))+1)/2],[-(sqrt((sqrt(2*(sqrt(17)-1))-2)/2)),(-(sqrt(17))+1)/2]]

However, the second expression in the system does not vanish for the fourth solution: the command

Code : Tout sélectionner

normal(subs(x1^4+2*x1^2+x2^2-3,[x1,x2],sol[3]))

yields

Code : Tout sélectionner

-(sqrt(17))

[parisse]

Normal is wrong (evalf returns approx 0), it's probably the same bug as your other report, I will investigate.

[lukamar]

Everything is fine in Xcas 1.7.0-35 in Windows.