If you use "solve(equat,dict=True", where the right side of equat must be zero, solve gives a list of a dictionary of results for all variables appearing in equat. The result looks like this: [{....}].