Finding roots polynomial 2D
I have polynomial 2D degree 2 $(1,x,y,xy,x^2,y^2) $ defined on triangle. How to find all roots?
add a comment
2 Answers
You can solve exactly in algebraic numbers (AA) and then get results in reals with any precision by using .n() method:
R.<x,y> = AA[]
eq1 = 2*x^2 + x*y +x - 6
eq2 = 3*x^2 + 2*y^2- y - 6
sols = R.ideal([eq1,eq2]).variety()
[tuple(map(lambda t: s[t].n(digits=20), R.gens())) for s in sols]
I found sequence:
var('x y')
eq1 = 2*x^2 + x*y +x == 6
eq2 = 3*x^2 + 2*y^2- y == 6
solutions = solve([eq1, eq2], x, y)
print(solutions)
now I want increase precision: R = RealField(100) how use with solve?
Your Answer
Stats
Asked: 2022-09-18 10:43:54 +0200
Seen: 91 times
Last updated: Sep 18 '22
