2015-10-13-163157.sagews

maxima_calculus("domain: real")
dnf = solve(abs((x-1)/(x-5))<1/3,x)
qf = apply(qepcad_formula.or_, map(qepcad_formula.and_, dnf)) # reformat the solution
qepcad(qf, vars='(x)') # simplify

real x - 2 < 0 /\ x + 1 > 0

%giac
solve(abs((x-1)/(x-5)) <= 1/3, x) 

list[((x>=-1) and (x<=2))]