Comment #2 on issue 2075 by smichr: solve fails for eq1=1 + 2*y/b - (e + x)**2/a; eq3=1 + 2*y/d - (x - e)**2/c
http://code.google.com/p/sympy/issues/detail?id=2075

Here is a single equation example that also fails:

h[1] >>> solve((x**2 - 1)**2 - a, x)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "sympy\solvers\solvers.py", line 241, in solve
    result = roots(poly, cubics=True, quartics=True).keys()
  File "sympy\polys\polyroots.py", line 449, in roots
    for root in _try_decompose(f):
  File "sympy\polys\polyroots.py", line 377, in _try_decompose
    g = factor - Poly(root, x)
  File "sympy\core\decorators.py", line 49, in __sympifyit_wrapper
    return func(a, sympify(b, strict=True))
  File "sympy\polys\polytools.py", line 2837, in __sub__
    return f.sub(g)
  File "sympy\polys\polytools.py", line 1001, in sub
    _, per, F, G = f._unify(g)
  File "sympy\polys\polytools.py", line 320, in _unify
    F = f.rep.convert(dom)
  File "sympy\polys\polyclasses.py", line 264, in convert
    return DMP(dmp_convert(f.rep, f.lev, f.dom, dom), dom, f.lev)
  File "sympy\polys\densebasic.py", line 508, in dmp_convert
    return dup_convert(f, K0, K1)
  File "sympy\polys\densebasic.py", line 483, in dup_convert
    return dup_strip([ K1.convert(c, K0) for c in f ])
  File "sympy\polys\domains\domain.py", line 85, in convert
    raise CoercionFailed("can't convert %s of type %s to %s" % (a, K0, K1))
sympy.polys.polyerrors.CoercionFailed: can't convert DMF(([1], [1]), ZZ) of type
 ZZ(a) to ZZ(a,a**(1/2))


--
You received this message because you are subscribed to the Google Groups 
"sympy-issues" group.
To post to this group, send email to sympy-iss...@googlegroups.com.
To unsubscribe from this group, send email to 
sympy-issues+unsubscr...@googlegroups.com.
For more options, visit this group at 
http://groups.google.com/group/sympy-issues?hl=en.

Reply via email to