z^n=conj(z)

[alb]

je ne pense pas qu'il sera possible de remédier à ce problème ?
simplifier(csolve(z^8=conj(z),z)) renvoie list[0,1,((i)*sqrt(3)-1)/2,((-i)*sqrt(3)-1)/2]
alors qu'on devrait avoir 0 et les racines neuvièmes de l'unité.

[parisse]

Le système passe en x/y et il appelle gsolve, qui ne renvoie pas toutes les solutions actuellement (les polynômes irréductibles de degré >=4 qui apparaissent ne sont pas factorisés). Mais je peux surement améliorer ça, dans le cas de 2 variables, on pourrait appeler directement le résultant en x et y et sortir des intervalles d'isolation des racines. A mettre dans la TODO list!