Very often it happens that using
solve alone, gives huge expressions that can't be used. The simplest thing to do is to wrap it inside an
evalf. But then sometimes, even your system only has real solutions, you may get some complex numbers. When this
evalf(solve()) being used inside an algorithm, then disastrous consequences may arise! If the system consists of a single equation of a single variable, then you may have some more tools. But if you have a system of several equations in several variables, you have less options. I am mostly interested in polynomials, I know several approaches to use and code to solve and get only the real solutions, but my codes might be not very optimized. In Maple 2022, one predefined command which is nice is
RootFinding:-Isolate but it has one issue and it is that this command only likes numeric coefficients which means integers, fraction of integers and float numbers, so no square root or other real numbers of this shape in the coefficients. I thought it might be a good idea to have a list of all solving commands in Maple that only return the real solutions or have the options to restrict to only real solutions.
fsolve is not very ideal, because it only returns one solution.