i1 : R = QQ[x,y]; S = {x^2 - 1/3, x*y - 1}; roots = phcSolve(S);
using temporary files /tmp/M2-5203-1PHCinput and /tmp/M2-5203-1PHCoutput
|
i4 : r0 = roots#0#Coordinates#1 o4 = 1.73205080756888 o4 : CC (of precision 53) |
i5 : newRoots = refineSolutions(S,roots,64)
using temporary file /tmp/M2-5203-3PHCoutput for storing refined solutions
o5 = {{.57735, 1.73205}, {-.57735, -1.73205}}
o5 : List
|
i6 : newRoots#0 -- recall that solutions are of type Point
o6 = {.57735, 1.73205}
o6 : Point
|
i7 : r1 = newRoots#0#Coordinates#1 o7 = 1.73205080756887729352744634150587236694280525381038062805580698 o7 : CC (of precision 213) |