Sage cannot solve non-linear polynomial systems symbolically!
x, y, z = var('x y z')
eq1 = x + y + z == 2
eq2 = x^2 + 2*y + 3*z == 4
eq3 = x + y + z^2 == 10
show(solve([eq1,eq2,eq3],x,y,z))
This gives me numerical solutions.
in Mathematica:
{eq1, eq2, eq3} = {x + y + z == 2, x^2 + 2*y + 3*z ==4, x + y + z^2 == 10};
Solve[{eq1, eq2, eq3}, {x, y, z}]
gives me symbolic solutions.
Any way to get symbolic solutions in sagemath ?
You mean solutions involving square roots etc. ( like 12(2−i√2(√33−1)) ) rather than floating point results, right? You should force the integers to be symbols. The following is a "impractical" workaround:
Tip: in the Ask Sage search box, type: polynomial system.