| LIB "solve.lib";
ring r = 0,(x,y),lp;
// compute the intersection points of two curves
ideal s= x2 + y2 - 10, x2 + xy + 2y2 - 16;
fglm_solve(s,10);
→ // name of new ring: rC
→ // list of roots: rlist
rlist;
→ [1]:
→ [1]:
→ 2.8284271247
→ [2]:
→ 1.4142135624
→ [2]:
→ [1]:
→ -2.8284271247
→ [2]:
→ -1.4142135624
→ [3]:
→ [1]:
→ 1
→ [2]:
→ -3
→ [4]:
→ [1]:
→ -1
→ [2]:
→ 3
|