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 callthe 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 ?