Status: New
Owner: ----
Labels: Type-Defect Priority-Medium
New issue 2228 by [email protected]: solve_poly_system: too complicated
result
http://code.google.com/p/sympy/issues/detail?id=2228
I found two equations, that solve_poly_system handles in strange manner.
eq1 = 15/2 - 5*y/2 - 3*x/4 + y**2/4 + x**2/16 # expand(Ellipse(Point(6, 5),
4, 2).equation())
eq2 = 261/4 - 16*x - 3*y/4 + x**2 + y**2/16 # expand(Ellipse(Point(8, 6),
1, 4).equation())
If i call "result = solve_poly_system([eq1, eq2], x, y)", i'll get very big
result, that simplify() can't cope.
[(x[0].evalf(), x[1].evalf()) for x in result] returns
[(7.02757595283879 - 4.5938932547515e-30⋅ⅈ, 6.93287917762256 +
1.53485614032609e-29⋅ⅈ), ...]
That is very close to result from maxima [[x = 7.027576197387518, y =
6.932879377431907], ...], but has unwanted imaginary part.
--
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.
