Numerical solve
When I use a dictionary with solve, the result is product with radicals. I have found no way in the documentation to obtain a decimal expension. But If I call the solution as sol = solve() it's not difficult to obtain directly the decimal expansion with sol[0].rhs().n(). So is there an unseen option in solve to obtain this and if not could it be implemented ?
Could you write down a complete example that actually shows the problem.