I am probably missing something about how to use sage:
If I run this:
x = var('x')
a = var('a')
solve([a*x>0,a>0],[x])
I get this results:
[[0 < x, a > 0], [x < 0, -a > 0, a > 0]]
However, I would expect something like:
[0 < x, a > 0]
What am I missing?
Thanks
Formally, there is nothing wrong since the second list of inequalities leads to the empty set (because of -a > 0, a > 0). That said, one could indeed expect that that list would have been automatically removed. In general, the solve function is not reliable since it might fail silently (e.g. give less solution than expected), for such systems of inequelities, an alternative is to use qepcad, whose installation is unfortunately currently broken (it used to be an experimental package).
Asked: 2015-05-26 11:38:05 +0100
Seen: 762 times
Last updated: May 27 '15
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
Reported as http://trac.sagemath.org/ticket/18902