Status: Accepted
Owner: plaes800
Labels: Type-Defect Priority-Medium
New issue 1637 by plaes800: solve([Eq(x+a*y,1), Eq(x+y+z, 1), Eq(z+3*x)],
x, y, z, a) troubles
http://code.google.com/p/sympy/issues/detail?id=1637
In [38]: x, y, z, a = symbols('xyza')
In [39]: solve([Eq(x+a*y,1), Eq(x+y+z, 1), Eq(z+3*x)], x, y, z, a)
---------------------------------------------------------------------------
PolynomialError Traceback (most recent call last)
/home/plaes/code/sympy/<ipython console> in <module>()
/home/plaes/code/sympy/sympy/solvers/solvers.pyc in solve(f, *symbols,
**flags)
327 soln = solve_linear_system(matrix, *symbols,
**flags)
328 else:
--> 329 soln = solve_poly_system(polys)
330
331 # Use swap_dict to ensure we return the same type as
what was
/home/plaes/code/sympy/sympy/solvers/polysys.pyc in
solve_poly_system(system,
*symbols)
113 head, tail = f.unify_with(system)
114
--> 115 solutions = solve_reduced_system([head] + tail, True)
116
117 if solutions is None:
/home/plaes/code/sympy/sympy/solvers/polysys.pyc in
solve_reduced_system(system,
entry)
70 f = univariate.pop()
71 else:
---> 72 raise PolynomialError("Not a zero-dimensional system")
73
74 zeros = roots(Poly(f, f.symbols[-1])).keys()
PolynomialError: Not a zero-dimensional system
Maple gives following solutions:
x = x; a = -(x-1)/(2x+1) ; z = - 3x and y = 2x +1
Also (maybe I'm doing it wrong):
In [43]: print solve([Eq(x+a*y,1), Eq(x+y+z, 1), Eq(z+3*x)], x)
-------> print(solve([Eq(x+a*y,1), Eq(x+y+z, 1), Eq(z+3*x)], x))
None
--
You received this message because you are listed in the owner
or CC fields of this issue, or because you starred this issue.
You may adjust your issue notification preferences at:
http://code.google.com/hosting/settings
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en
-~----------~----~----~----~------~----~------~--~---