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) |