Subject: error in solving x^4> -2 via solve_rat_ineq
From: Robert Marik
Date: Sun, 24 Jan 2010 19:22:01 +0000 (UTC)
Volker van Nek <volkervannek <at> googlemail.com> writes:
>
> Hi Robert,
> This is one of the examples which show that it might lead to trouble if
> you want to compare a float to zero. But simplification helps. In my
> patch for solve_rat_ineq I am going to use
> is(equal(imagpart(x),0))
> as the test. is(equal(a,b)) internally uses ratsimp(a-b).
Thank you for detailed explanation.
>
> (I also had a look at the underlying algsys bug. At first sight the fix
> seems to be quite easy.)
Great, thanks.
Robert