I ran your code exactly as it's posted here and got [{s: 1.51111111111111, q: 0.0, z: 0.0, p: 0.0}, {s: 1.51111111111111, q: 0.0 , z: 0.0, p: 0.980392156862745}]
I'm using SymPy 1.1.1. Perhaps your version of SymPy is old? By the way, it's better to use the assumption sm.symbols('p, q, z, s', nonnegative=True) although this doesn't change the result. There is a difference: negative=False just says it cannot be negative, but it might be a complex number like 2+3*I. While nonnegative=True means it's a real number >=0. It's rare that setting an property to False is needed, with exception of commutative=False (and that is because commutative is True by default.) -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To post to this group, send email to sympy@googlegroups.com. Visit this group at https://groups.google.com/group/sympy. To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/0271c2dc-d87f-4844-825f-1c4cef8da9db%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.