I found that the solve command doesn't work properly with some systems of equations containing radicals. For xample: solve({b = RootOf(_Z^3+c^2+1),c = RootOf(_Z^2+_Z+1), b = RootOf(_Z^3-RootOf(_Z^2+_Z+1))}); doesn't return any results and no "_SolutionsMayBeLost" warning. But using solve(convert({b = RootOf(_Z^3+c^2+1),c = RootOf(_Z^2+_Z+1), b = RootOf(_Z^3-RootOf(_Z^2+_Z+1))},radical)); solves the system. What can be done about it to avoid such situations?