restarteqs := [u + v + w = 1,
u*x + v*y + w*z = 1/2,
u*x^2 + v*y^2 + w*z^2 = 1/3,
u*x^3 + v*y^3 + w*z^3 = 1/4,
u*x^4 + v*y^4 + w*z^4 = 1/5,
u*x^5 + v*y^5 + w*z^5 = 1/6]soln := solve(eqs, explicit)# Show that each member of soln is a solution of eqsfor s in soln do simplify(subs(s, eqs)) end do# Show that the solutions are unique, if the union of each solution and each remaining solution is not equal to nops(soln[1])for idx to 6 do for jdx from idx+1 to 6 do print(verify(nops(soln[idx]),nops(soln[idx] union soln[jdx]), equal)) end do; end donops(soln[1] union soln[2])