| 1 | initial version |
Instead of using solve, you could use solve_ineq which is more appropriate. So, you can replace:
sol=[solve(ieq[0],x[0]) for ieq in Ineq]
with
sol=[solve_ineq(ieq,[x[0]]) for ieq in Ineq]
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.