I want to substitute the solution back into the original equation. I get caught up in RootOf and have to manually do the substitutions.
F := [x^2+y+z-1, y^2+x+z-1, z^2+x+y-1];
soln1 := solve(F);
for s in soln1 do
subs(s,F)
end do;
The 4th soln has RootOf.
soln2 := solve(_Z^2 + 2*_Z - 1);
for s in soln2 do
evala(subs({x=s,y=s,z=s},F))
end do;
How do I do this all in one step?