Thank you for these hints. I wonder what Sage uses to get the current
result. It must be none of these.
As I am unsure whether it is better to focus on Maxima or Sage w.r.t. my
test case, I have also posted it to the Sage development group:
https://groups.google.com/forum/#!topic/sage-devel/QLEFLgkuXRg
Robert
Am 28.06.2013 20:09, schrieb Volker van Nek:
> Another one is fourier_elim:
>
> (%i1) load(fourier_elim);
> (%o1)
> /usr/local/share/maxima/branch_5_30_base_98_g29f9239_dirty/share/fourier_elim/fourier_elim.lisp
>
> (%i4) fourier_elim([abs(x-1)/abs(x-5)<=1/3],[x]);
> (%o4) [x = - 1] or [x = 2] or [x = 1] or [1 < x, x < 2] or [- 1 < x, x < 1]
>
> Volker van Nek
>
>
>
> 2013/6/28 Stavros Macrakis <macrakis at alum.mit.edu
> <mailto:macrakis at alum.mit.edu>>
>
> Fortunately, Maxima /does/ have a package that solves your equation.
>
> load(to_poly_solve)$
>
> %solve(abs((x-1)/(x-5))<=1/3,x)
> =>
> %union([?1<x,x<1],[1<x,x<2],[x=?1],[x=1],[x=2])
>
>
> This is still not ideal -- the various cases come down do -1 <= x, x
> <=2 -- but still, better than what you had.